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Abstract 



Results are presented from analyses of jet data produced in pp collisions at y/s 
= 630 and 1800 GeV collected with the D0 detector during the 1994-95 Fer- 
milab Tevatron Collider run. We discuss details of detector calibration, and 
jet selection criteria in measurements of various jet production cross sections 
at y/s = 630 and 1800 GeV. The inclusive jet cross sections, the dijet mass 
spectrum, the dijet angular distributions, and the ratio of inclusive jet cross 
sections at y/s = 630 and 1800 GeV are compared to next-to-leading-order 
QCD predictions. The order calculations are in good agreement with the 
data. We also use the data at yfs = 1800 GeV to rule out models of quark 
compositeness with a contact interaction scale less than 2.2 TeV at the 95% 
confidence level. 
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I. INTRODUCTION 



The quark model which suggested that hadrons are composite particles was first proposed 
in the early 1960's and was confirmed as the quark-parton model in a series of exper- 
iments at the Stanford Linear Accelerator Center in the late 1960's and early 1970's |§. 
The model developed during the 1970's into the theory of strong interactions, Quantum 
Chromodynamics (QCD) ||, which describes the interactions of quarks and gluons (called 
partons). Together with the theory of electroweak interactions, QCD forms the foundation 
of the standard model of particle physics (SM), which, thus far, describes accurately the 
interactions of all known elementary particles. 

Perturbative QCD (pQCD) [||] predicts the production cross sections at large transverse 
momentum (p?) for parton-parton scattering in proton-antiproton (pp) collisions Pj-B|. The 
outgoing partons from the parton-parton scattering hadronize to form jets of particles. High- 
Pt jets were observed clearly during experimentation at the CERN ISR || and the CERN 
pp collider jnj. Significant deviations from predictions of pQCD can only be observed if 



the uncertainties in both experimental measurements and theoretical predictions are small. 
Calculations of high-pr jet production involve the folding of parton scattering cross sections 
with experimentally determined parton distribution functions (PDFs). These predictions 
have recently improved with next-to-leading-order (NLO) QCD calculations [Tl -13] and im- 
proved PDFs Hl4| , |l5| . These (D(a s 3 ) predictions reduce theoretical uncertainties to ~30% fl6| 
(where a s is the strong coupling parameter). 

In this paper we describe the production of hadronic jets as observed with the D0 
detector at the Fermilab Tevatron pp collider at center-of-mass (CM) energies of 630 and 
1800 GeV. High p T jet production at ^Js = 1800 GeV probes the structure of the proton 
where the interacting partons carry a fraction of the proton momentum, 0.1 ^ x ^ 0.66, 
for momentum transfers (Q) of 2.5 x 10 3 < Q 2 < 2.3 x 10 5 GeV 2 , where Q 2 = E\ and is 
equivalent to a distance scale of 10 -4 fm (see Fig. [I]). Measurements of the inclusive jet cross 
section, the dijet angular distribution, and the dijet mass spectrum, can be used to test the 
predictions of pQCD. Additionally, new phenomena such as quark compositeness [0 would 
reveal themselves as an excess of jet production at large transverse energy (E T ) and dijet 
mass (M) relative to the predictions of QCD. 

Previous measurements by the CDF Collaboration of the inclusive jet cross section 



and the inclusive dijet mass spectrum |L8| have reported an excess of jet production relative 
to a specific QCD prediction. More recent analysis of the dijet angular distribution by the 
CDF Collaboration |19] has excluded models of quark compositeness in which the contact 
interaction scale is less than 1.6 TeV at the 95% confidence level. 

This paper presents a detailed description of five measurements previously published 
by the D0 Collaboration: the inclusive jet cross sections at y/s = 1800 GeV |8|j20|| and 
630 GeV [2l]J22] |, the ratio of inclusive jet cross sections at y/s = 630 and 1800 GeV |2~T| , |2"2" 



the dijet angular distribution p3[^4|, and the dijet mass spectrum [^5[ at ^fs = 1800 GeV. 
In addition to the analyses presented in this paper, D0 has recently measured the inclusive 
jet cross section as a function of Et and pseudorapidity, \rj\, at y/s = 1800 GeV |2f| where 
7] = — ln[tan(#/2)] and 6 is the polar angle. 

For jet measurements, the most critical component of the D0 detector is the calorime- 



ter [28|. A thorough understanding of the jet energy scale, jet resolutions, and knowledge of 
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FIG. 1. The x and Q 2 range of the data set analyzed by the D0 experiment at \/i = 1.8 TeV 
(D0 [8] and CDF [7] Inclusive Jets with | rj\ < 0.7) compared with the data used to produce PDFs 
[27]. Also shown is the extended x and Q 2 reach of the D0 measurement of the inclusive jet cross 
section with 1 77 1 < 3 as presented in Ref . [26] . 
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biases caused by jet triggering and reconstruction are necessary. After detector calibration, 
small experimental uncertainties can be achieved and precise statements can be made about 
the validity of QCD predictions. These results can then be used as the basis for searches for 
physics beyond the standard model. 

In this paper we discuss the theoretical predictions for the inclusive jet cross section, 
the inclusive dijet mass cross section, and the dijet angular distribution. We describe the 
various measurements undertaken to understand and calibrate the D0 detector for jet mea- 
surements. Finally, four different physics measurements performed at D0 are presented: 
the inclusive jet cross section, the ratio of inclusive jet cross sections, the dijet angular 
distribution, and the dijet mass spectrum. The measurements presented here constitute a 
stringent test of QCD, with a total uncertainty substantially reduced relative to previous 
measurements. Further, the results represent improved limits on the existence of phenomena 
not predicted by the standard model. 
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II. CALORIMETER 

The D0 detector is described in detail elsewhere [P^ , p9|1 . The D0 uranium-liquid argon 
sampling calorimeters are uniform in structure and provide coverage for a pseudorapidity 
range 1 77 1 < 4.5. They are nearly compensating with an e/7r ratio of less than 1.05 above 30 
GeV. The central and end calorimeters are approximately 7 and 9 interaction lengths deep 
respectively, ensuring containment of most particles except high-p^ muons and neutrinos. 
The calorimeters are segmented into cells of size A7/ x A0 = 0.1 x 0.1, where is the 
azimuthal angle. These characteristics along with excellent energy resolution for electrons 
(~ 15% / J E[GeV\) and pions (~ 50%/ yj E [GeV]) make the D0 calorimeters especially well 
suited for jet measurements. 

The calorimeter is divided into three sections (see Fig. ^): a central calorimeter (CC) 
covering low values of pseudorapidity, two end calorimeters (EC) covering forward and back- 
ward pseudorapidities, and the Intercryostat Detector (ICD) covering the gaps between the 
CC and EC (0.8 < \ f] \ < 1.6). The CC and EC calorimeters each consist of an inner electro- 
magnetic (EM) section, a fine hadronic (FH) section, and a coarse hadronic (CH) section. 
Each EM section is 21 radiation lengths deep and is divided into four longitudinal samples. 
The hadronic sections are divided into four (CC) and five (EC) layers. The ICD consists 
of scintillator tiles inserted in the space between the CC and EC cryostats. The Tevatron 
accelerator's Main Ring, which is used for preaccelerating protons, passes through the CH 
calorimeters. 



0.0 0.5 1.0 




FIG. 2. A schematic view of one quarter of the D0 calorimeters. The shading pattern indicates 
the distinct readout cells. The rays indicate the pseudorapidity intervals defined from the center 
of the detector. 
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III. JET DEFINITIONS 

A jet is a collection of collimated particles produced by the hadronization of a high- .Bp 
quark or gluon. In the measurements presented in this paper, we measure the energy and 
direction of the jets produced in pp interactions and compare the measurements to various 
theoretical predictions. 

In addition to the jets produced by the high-p-^ parton-parton scattering, there are many 
particles produced by the hadronization of the partons in the proton and antiproton that 
were not involved in the hard scattering process. Because of this there is no unequivocal 
method for experimentally selecting the particles that belong to a jet produced in high-p-r 
scattering. It is preferable to use a standard definition of a jet to facilitate comparisons 
of measurements from different experiments, and with theoretical predictions. In 1990 the 



so-called Snowmass Jet Algorithm [30] was adopted as a standard definition. 

A jet algorithm can be run on several different input variables, calorimeter cells, and 
particles or partons produced by a Monte Carlo Simulation. To differentiate the results 
of the same algorithm being run on these different input we describe the resulting jets as 
follows: A jet (or calorimeter jet) is the result of the jet algorithm being run on calorimeter 
cells; A particle jet is created from particles produced by a Monte Carlo simulation after the 
hadronization step; Finally, a parton jet is formed from partons before hadronization takes 
place. 

A. The Snowmass Accord 

The Snowmass Jet Algorithm defines a jet as a collection of partons, particles, or 
calorimeter cells contained within a cone of opening angle 1Z. All objects in an event have 

a distance from the jet center, TZi = \J (rji — Jftet) 2 + {4>i ~ 0jet) 2 > where rj- jet and je t define 
the direction of the jet and (rji, fa), are the coordinates of the parton, particle, or center of 
the calorimeter cell. If TZi < TZ then the object is part of the jet. The Snowmass suggested 
value of 1Z = 0.7 was used in these measurements. The E T of the jet is given by 

Et= E E t, (3-1) 

i£lli<1l 

where i is an index for the i th parton or cell. The direction of the jet is then given by 

Viet = ^r E E W, (3.2) 

T ieTl^ll 

The Snowmass algorithm gives a procedure for finding the jets: 

• Determine a list of jet "seeds," each with a location r]j et , fa 3et . 

• Form a jet cone with direction r/j et , </>j e t- 
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• Recalculate the E^ and direction of the jet. 

• Repeat steps 2 and 3 until the jet direction is stable. 

The definition of the jet seed is not given by the algorithm. At the parton level these seeds 
could be the partons, points lying between pairs of partons, or even a set of points randomly 
located in rj-cfi space. Experimentally, the seed could be defined as any cell above a given 
i?T threshold, all cells in the calorimeter, or clusters of calorimeter cells. It is up to each 
experimentalist or theorist to define a seed. 



B. The D0 Experiment's Jet Algorithm 

At the calorimeter level in the D0 experiment, jets are defined in two sequential proce- 
dures. In the first, or clustering, procedure all the energy depositions that belong to a jet 
are accumulated, and in the second the r], <p, and of the jet are defined. The clustering 
consists of the following steps: 

1) Calorimeter towers (a set of four calorimeter cells of size A77 x A<ft = 0.2 x 0.2) with 
-Sr > 1 GeV are ordered in E^. Starting with the highest-.&r tower, preclusters are 
formed from contiguous towers around these seed towers. 

2) The jet direction (77j et ,^>j et ) is calculated using Eq. |3]2| from the energy deposit pattern 
in a fixed cone of size TZ around the precluster center. 

3) The energy deposited in a cone of size 1Z around the jet axis is summed and the jet 
direction (?7jet,0jet) is recalculated using the Snowmass algorithm (Eq. [3.2|) . 

4) Step 3 is iterated until the jet direction is stable. This is typically achieved in two or 
three iterations. 

5) Only jets with E^ > 8 GeV are retained. 

6) Jets are merged or split according to the following criteria: two jets are merged into 
one jet if more than 50% of the E^ of the jet with the smaller E^ is contained in 
the overlap region. If less than 50% of the E^ is contained in the overlap region, 
the jets are split into two distinct jets and the energy of each calorimeter cell in the 
overlap region is assigned to the nearest jet. The jet directions are recalculated using 



an alternative definition as given in Eq. 3.3 



The D0 jet algorithm and the Snowmass algorithm calculate the final direction of the jet 
differently. In the D0 jet algorithm the final 77 and <fi of the jet are defined as: 



jot = tan 1 



YE i 



Viet = - In [tan(0 jot /2)] 



, t = tan- (3.3) 
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where i corresponds to all towers whose centers are within the jet radius 71, 
E l x = Ei sin 8i cos <j> i7 E l y = Ei sin 6i sin 0j, and E l z = E^ cos Q{. 

Applying the 8 GeV Et threshold to jets before merging and splitting has two important 
consequences. The first is that two jets of < 8 GeV cannot be merged into a single jet 
to create a jet with E^ > 8 GeV. The second is that jets may have E^ < 8 GeV if they were 
involved in splitting. 



C. Corrected Jets 

In this paper a "true" or corrected jet is the jet that would be found by the D0 jet 
algorithm if it was applied to the particles produced by the high-p T parton-parton scat- 
tering before they hit the calorimeter. The jet does not include the particles produced by 
hadronization of the partons not involved in the hard scattering (the underlying event). The 
differences between jets observed in the calorimeter and the "true" jets are determined using 
Monte Carlo (MC) simulations of pp interactions. The direction and E^ of the "true" jets 
are calculated using the Snowmass definition (see Eqs. |3.1| and ^.2| ) and are denoted by E^ cl , 
77 ptcl , and ptcl (where ptcl denotes particle). 



D. Differences between the D0 and Snowmass Algorithms in Data 

Because the D0 and Snowmass algorithms calculate the location, and hence angle, of 
the jet differently, a study to measure the differences was performed. The same data events 
were reconstructed using the two different algorithms and the differences in location were 
compared. There were no differences in the or <p of the jets. However, there were small 
differences in the jet \ r]\, which increase as a function of the 1 77 1 of the jets and decrease as 
the transverse energy of the jets increases. Figure |3| shows the average difference between 
the 1 77 1 of jets with E<? > 40 GeV reconstructed using the two different algorithms. As can 
be seen, the difference is small, even at a large \ rj\. 



E. Jet Algorithms at NLO 

In pQCD calculations of parton-parton scattering at leading order (LO, 0(a s 2 )) there 
can only be two partons in the final state. These partons are well separated and always 
form two jets when the Snowmass algorithm is applied. At next-to-leading order (NLO, 
0(a s 3 )), three partons can be produced in the final state. The Snowmass algorithm at the 
NLO parton level is illustrated in Fig. §|(a). F° r an y two partons in the final state, the seeds 
direction is given by applying Eq. [O] to the two partons. If the partons lie within a distance 
TZ from the seed's the two partons are combined to form a jet. 

In the Snowmass algorithm the partons contributing to a single jet can have a maxi- 
mum separation of 27Z. Consider a two parton final state with the partons separated by 
2TZ. The experimentally observed energy pattern will be determined by the parton show- 
ering, hadronization, and calorimeter response. Application of the D0 jet algorithm to 
the calorimeter energy deposition that results from the hadronization of the two partons 
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FIG. 3. The average difference between the 1 77 1 of jets reconstructed using the D0 algorithm 
and the Snowmass algorithm for D0 data. 



will produce one or two jets depending on the splitting and merging criteria. The Snow- 
mass algorithm is only capable of finding one jet, and hence cannot match the experimental 
measurement. 

This example illustrates the different treatment of jets at the parton and calorimeter level. 
To accommodate the differences between the jet definitions at the parton and calorimeter 
levels, an additional, purely phenomenological parameter has been suggested in Ref. [£31 



The variable is called !Z sep and is the maximum allowed distance (A7?.) between two partons 
in a parton jet, divided by the cone size used: lZ sep = /S.1Z/1Z. This algorithm is illustrated 
in Fig. |](b) and is referred to as the modified Snowmass algorithm. 

The value of lZ sep depends on details of the jet algorithm used in each experiment. At the 
parton level 7Z sep is the distance between two partons when the clustering algorithm switches 
from a one-jet to a two-jet final state, even though both partons are contained within the 
jet defining cone. The value of 7Z sep depends on the experimental splitting/merging scheme. 
After several studies an 7Z sep value of 1.3 was found to best simulate the D0 merging and 



splitting criteria [32 



F. Jet Reconstruction Efficiency 

The jet algorithm does not reconstruct all jets with the same efficiency Primarily this 
is due to calorimeter energy clusters not containing a seed tower of greater than 1 GeV. 
Since the jet algorithm explicitly depends on the Et of a seed tower used to begin searching 
for a jet, the seed tower distributions are studied to determine if jets are likely to have seed 
towers below threshold. Figure [5] shows the seed-tower- E^ distribution of jets for an E^ 
range of 18 to 20 GeV (other E^ ranges have similar distributions). The distribution is 
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FIG. 4. Illustration and description of the jet definitions at NLO parton level as used by the 
D0 experiment, a) The jet definition in NLO according to Snowmas. Parton -1- and -2- are 
combined into jet -J-, if the parton distance to the jet axis is less than R. The jet axis is defined by 
partons 1 and 2, according to the Snowmass definition, b) The jet definition in NLO according to 
the modified Snowmass with lZ sep - Use the standard Snowmass clustering, but in addition require 
the distance between the two partons to be less than lZ sep . 
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fitted with 



A ■ exp 



-0.5- (-W + e 

y/2-X 



(3.4) 



where W = (E^ — x)/\ an d A V an d ^ are f ree parameters in the fit. Assuming that the 
seed towers are smoothly distributed in E T , the fraction of jets not containing a seed tower 
exceeding 1 GeV is determined from the fit and used to calculate the jet reconstruction 
efficiency. Figure ^| shows the reconstruction efficiency for jets as a function of jet E^. For 
jets with an Ey of 20 GeV and 72. = 0.7, the reconstruction efficiency is 99%. 
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FIG. 5. Seed tower distributions for 1Z = 0.7 cone jets with an £/p range of 18-20 GeV. The 

data is represented by the solid histogram and the fit is given by the dashed curve. 



G. Biases in the Jet Algorithm 

The rj dependence of the calorimeter energy response together with algorithm related 
effects may cause a bias in the reconstructed jet rj. The r\ of the jet, assuming perfect 
position resolution, is: 

^P tcl = v + p(E, rj) (3.5) 



where p(E, rj) is the possible bias. To measure the bias, the herwig [33| Monte Carlo 



event generator and the D0 detector simulation, D0GEANT |34]| , are used. Jets are recon- 



structed at both the particle and calorimeter level. Statistically, (p(E,r))) can be obtained 
as (rj ptcl — rj) where a matching criterion is used to associate the particle jet to the recon- 
structed calorimeter jet. Figure |7| shows the r] bias for all jet energies as a function of rj. 
The bias in 77 is less than 0.02 for all rj. The magnitude of the bias is greatest when part of 
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FIG. 7. herwig Monte Carlo simulation of the r\ bias for all jet energies as a function of the 
reconstructed jet rj. 
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the jet falls into the Intercryostat Region (0.8 < r\ < 1.6), which is the least instrumented 
region of the calorimeter. 

A similar study was performed to measure a possible bias in (azimuth). Since the 
calorimeter has a symmetric tower structure in 0, no bias is expected. The bias in 
(azimuth) was measured to be small — much less than 0.01 radians. Any bias introduced 
by this effect will be small for most physics analyses since the ATI between jets is typically 
used rather than the absolute or r] position. The analyses presented in this publication 
are not corrected for these effects. 
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IV. THEORETICAL PREDICTIONS 

A. NLO QCD Predictions 

Within the framework of pQCD, high jet production can be described as an elastic 
collision between a single parton from the proton and a single parton from the antiproton [|3| . 
After the collision, the outgoing partons form localized streams of particles called "jets." 
Predictions for the inclusive jet cross section, the dijet angular distribution, and the dijet 
mass spectrum are in general given by H]: 



o = Yl J dx!dx 2 fi fj (^2,^1) 



'j 
a 



( i\ 

X 1 Pi,X 2 P2,Otg [H R ) , — , — 

/Af ^r. 



(4.1) 



where (x\(2)i I^f) represent the PDFs of the proton (antiproton) defined at factorization 
scale up, & is the parton scattering cross section, P\(i) is the momentum of the proton (an- 
tiproton), X\(2) is the fraction of the proton (antiproton) momentum carried by the scattered 
parton, Q is the hard scale that characterizes the parton scattering (which could be the 
of a jet, the dijet mass of the event, etc.), and fiR is the renormalization scale. 

The parton scattering cross sections have been calculated to next-to-leading order (NLO, 
0(a s 3 )). The perturbation series requires renormalization to remove ultraviolet divergences. 
This introduces a second scale to the problem, In addition, an arbitrary factorization 
scale, /ip, is introduced to remove the infrared divergences. Qualitatively, it represents the 
scale that separates the short- and long-range processes. A parton emitted with transverse 
momentum relative to the proton less than the scale fip will be included in the PDF, while 
a parton emitted at large transverse momentum will be included in a. The scales hr and 
should be chosen to be of the same order as the hard scale, Q, of the interaction. The larger 
the number of terms included in the perturbative expansion, the smaller the dependence on 
the values of [Ir and /zp. If all orders of the expansion could be included, the calculation 
should have no dependence on the choice of scales. In this article the renormalization 
scale is written as the product of a constant, D, and the hard scale of the interaction, 
/i = DQ. Typically, the renormalization and factorization scales are set to the same value, 
— Hr — 1*>f- 

Several pQCD NLO calculations have been performed |TT|--[T3|. In this paper we use 
the event generator jetrad and a version of the analytic calculation eks [^j that 
integrates the cross section over bins. Both programs require the selection of a renormal- 
ization and factorization scale, a set of parton distribution functions, and a jet clustering 
algorithm. Two partons are combined if they are contained within a cone of opening angle 
1Z = V At]' 2 + Acf) 2 = 0.7, and are also within lZ sep = 1.3 (see Section |TJ). The authors of 
jetrad have provided several choices for the renormalization scale. We have chosen a scale 
proportional to the E^ of the leading jet: \x = DE^ &X , where D is constant in the range 
0.25 < D < 2.00. An alternative scheme sets the scale to be proportional to the center-of- 
mass energy of the two outgoing partons: \i = C\fl = C ^JX\XiS where C is constant in the 
range 0.25 < C < 1.00, x\ = Yl E^e 11 * />/s, x 2 = Enz~ m I yfs, and i = 1, . . . ,n where n 
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is the number of jets in the event. The authors of eks prefer a different definition of the 
renormalization scale: the of each jet in the event, /i = DE 3 ^ 1 (a version of eks that uses 
the renormalization scale \i = DE™* is also available). 

Several choices of PDF are considered. The CTEQ3M H and MRS (A') g| PDFs are 
fits to collider and fixed target data sets published before 1994. CTEQ4M ]TJJ updates 
these fits using data published before 1996, and CTEQ4A repeats the fits with values of 
a s (M z ) fixed in the range 0.110 to 0.122 (CTEQ4M corresponds to an a s (M z ) of 0.116). 
CTEQ4HJ |14]| adjusts the gluon distributions to fit a CDF inclusive jet cross section mea- 
surement by increasing the effective weighting of the CDF data. MRST incorporates 
all data published before 1998. In addition to the standard MRST PDF, two alternative 
PDFs are provided that vary the gluon distributions within the range allowed by experi- 
mental observations. The resulting distributions are called MRST(gf)and MRST(gj). 

Since the publication of the MRST and CTEQ4 PDFs, problems were found in the 
implementations of the QCD evolution of the parton distributions in Q 2 f37|. This was 



caused by approximations to NLO QCD to reduce the time required for computation. The 
removal of these approximations could lead to changes of approximately 5% in the theoretical 
predictions presented in this paper. Currently, PDFs calculated without the approximations 
are not available for use with jetrad and eks. 



1. Inclusive Jet Cross Section 

The inclusive jet cross section may be expressed in several ways. Theoretical calculations 
are normally expressed in terms of the invariant cross section 

In the D0 experiment the measured variables are the transverse energy (Et) and pseudo- 
rapidity (77). In terms of these variables, the cross section is expressed as 

(4.3) 



dE T dr] 



where the two are related by 



d 3 cx d 3 a 1 d 2 a , A A . 

E = > , (4.4) 

dp 3 d 2 pTdy 2ttE t dE^dr/ ' 

where y is the rapidity of the jet. The final expression follows if the jets are assumed 
to be massless. For most measurements, the cross section is averaged over some range of 
pseudorapidity: in this paper \r) \ < 0.5 and 0.1 < \r) \ < 0.7. 

The inclusive jet cross section measures the probability of observing a hadronic jet with 
a given and r\ in a pp collision. The term "inclusive" indicates that the presence or 
absence of additional objects in an event does not affect the selection of the data sample. 
An event which contains three jets that pass the selection criteria, for instance, will be 
entered into the cross section three times. The inclusive measurement is sometimes denoted 
cr(pp-+ jet + X). 
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Theoretical predictions for the inclusive jet cross section are generated using the jetrad 
and eks programs. Our reference prediction is the jetrad calculation with \i = 0.5i?™ ax , 
TZ sep = 1.3, and the CTEQ3M PDF. The predictions are smoothed by fitting to the function 



ae; 



1 - 



2E^ 



(4.5) 



where Vq (Et) is a sixth order polynomial. The resulting uncertainty due to smoothing is 
less than 2% for a given E^. The uncertainty in the calculations resulting from the choices of 
different renormalization scales and PDFs is approximately 30% and varies as a function of 
Et [16|. Figure H shows the variations in the predictions for the inclusive jet cross section at 
s/s = 1800 GeV for jetrad. The uncertainties in the inclusive jet cross section at y/s = 630 
GeV are of a similar size. 
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FIG. 8. The difference between alternative predictions and the reference prediction 
(fi = 0.5£^ ax , CTEQ3M) for the inclusive jet cross section for | rf ct \ < 0.5 at y/s = 1.8 TeV. Shown 
are the alternative predictions for the choices (a) fi = (0.25, 0.75, 1.0, 2.0)£'™ ax , (b) fi = (0.25, 0.5, 
l.Ojv 7 ! and 0.5^ ct , (c) CTEQ4M, CTEQ4HJ, MRS(A'), and MRST, and (d) for a s = 0.110-0.122 
using the CTEQ4A PDFs compared with the calculation using CTEQ4M, for which a s = 0.116. 



2. Ratio of Inclusive Jet Cross Sections at yfs = 1800 and 630 GeV 

While it is possible to compare the inclusive jet cross sections as a function of E? for 
both center-of-mass energies, the data will differ greatly in both magnitude and E^ range 
(see Fig. |9|(a)). If we express the cross section as a dimensionless quantity 

4 d 3 a _ El d 2 a 
ETE df = 2^dE^> (4 ' 6) 
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and calculate it as a function of xt = 2Et/^, the "scaling" hypothesis, which is motivated 
by the Quark-Parton Model, predicts that it will be independent of the center-of-mass energy. 
However, QCD leads to scaling violation through the running of a s and the evolution of the 
PDFs. 
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FIG. 9. Different presentations of the NLO inclusive jet cross sections at (dashed line) 
= 1800 and (solid line) 630 GeV. Theoretically, the scaled dimensionless cross sections (b) 
should be nearly exponential and close to one another. Panel (c) shows the ratio of dimensionless 
inclusive jet cross sections at ^/s = 630 and 1800 GeV for \ rf et \< 0.5. 



By taking the ratio of the cross sections at = 1800 and 630 GeV, many of the 
theoretical and experimental uncertainties are reduced. Variations in the prediction resulting 
from the choice of renormalization scale, factorization scale, and PDF are approximately 10% 
and vary as a function of xt- This is a significant reduction in the theoretical uncertainty 
compared to the uncertainties in the inclusive jet cross sections. The theoretical predictions 
for the ratio of the inclusive jet cross sections are calculated using the jetrad and eks 
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programs. Figure [H] shows the variations in the ratio between inclusive jet cross sections at 
fa = 630 and 1800 GeV for jetrad. 
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FIG. 10. The difference between alternative predictions and the reference prediction 
= 0.5£f? ax , CTEQ3M) of the ratio of inclusive jet cross sections at *fs = 630 and 1800 GeV 
for \rf ct \ < 0.5. Shown are the alternative predictions for the choices (a) fi = (0.25, 0.75, 1.0, 
2.0)£^p ax , (b) ix = (0.25, 0.5, 1.0)Vs and 0.5£$*, (c) CTEQ4M, CTEQ4HJ, MRS(A'), and MRST, 
and (d) for a s = 0.110 — 0.122 using the CTEQ4A PDFs compared with the calculation using 
CTEQ4M, for which a s = 0.116. 



3. Dijet Angular Distributions 

At leading order two jets are produced. The invariant mass of the jets is given by 

M 2 = s = 4p 2 T cosh 2 (Ay/2) (4.7) 

where s = X1X2S, is the CM energy squared of the interacting partons, and Ay is the 
separation in rapidity of the two jets. If we assume that the jets are massless we can write 
the dijet invariant mass as 

M 2 = 2E$ tl E% t2 [cosh (Arj) - cos (A0)] , (4.8) 

where is the azimuthal angle with respect to the beam. Since the dijet mass represents 
the center-of-mass energy of the parton-parton interaction, it directly probes the parton 
scattering cross section. The presence of higher-order processes can result in the production 
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of additional jets. In this case the mass is calculated using the two highest-.&r jets in the 
event. 

If only two partons are produced in a parton-parton interaction, and we neglect the 
intrinsic transverse momentum of the scattering partons, the two jets will be back-to-back 
in azimuth and balance in transverse momentum. The resulting two-jet inclusive cross 
section at LO can be written as a function of the pr and rapidity (y 3 , y 4 ) of the jets || 



a 



dy z dydp 



(4.9) 



This can be rewritten in terms of the dijet invariant mass and the center of mass scattering 
angle, 6*, using the transformation || 



resulting in 



dp^dy^dy^ = Adx\dx 2 d cos 9* (4. 10) 



>l ' a Y f 1 d Xl dx 2 5 ( Xl x 2 s - M 2 ) 4^7T- (4-11) 

~? Jo v ' d cos u* 



dMd cos 6* 



The dijet angular distribution as measured in the dijet center-of-mass frame is sensitive 
to the QCD matrix elements. Angular distributions for the qg — > qg, qq — > qq, and gg — > gg 
processes are similar. The properties of parton-parton scattering are almost independent of 



the partons involved (see Fig. [11]). The dominant process in QCD parton-parton scattering 
is t-channel exchange, which results in angular distributions peaked at small center-of-mass 
scattering angles. Many theoretical predictions for phenomenology beyond the SM have 
an isotropic angular distribution and could be detected using the measurement of the dijet 
angular distribution. 

At small center-of mass-scattering angles, 8*, the dijet angular distribution predicted by 
leading order QCD is proportional to the Rutherford cross section: 

da 1 , 

(4.12) 



dcosfl* sin 4 (#*/2) 

It is conventional to measure the angular distribution in the variable \i rather than cos#*, 
where: 

X= t— P^U expflAql). (4.13) 
1 — | costr | 

Plotting the dijet angular distribution in the variable x flattens out the distribution and 
facilitates comparison to theory [[| (dajdx is uniform for Rutherford scattering). The dif- 
ferential angular cross section measured in this analysis is: 



d" 



dM dx drj\), 



" , (4.14) 



oost 



where r]^ oost = 0.5(r/i + 772)- The predictions are calculated using jetrad. 
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FIG. 11. Quark-antiquark and quark-gluon angular distributions, normalized to the angular 
distribution for gluon-gluon scattering. 



4- Inclusive Differential Dijet Mass Cross Section 

The inclusive triple differential dijet mass cross section is obtained by integrating over 
cos 8* and is given by 

dMdr]i etl d^ ct2 ^ 4 ' 15 ' ) 

where rf et1 ' 2 are the pseudorapidities of the jets. We integrate the cross section over a range 
of pseudorapidity such that both jets satisfy | rf et \ < 1.0. The NLO predictions for this cross 
section are calculated using jetrad. The jetrad predictions were smoothed by fitting 
them to an ansatz function of the form 



A-M~ a exp -[3M - ~fM 2 - 5M 3 V n (M) 



(4.16) 



where V n (M) is a polynomial of degree n < 6 and a, (3, 7, and 5 are fit parameters. The 
uncertainty due to the form of the ansatz function not being quite right is estimated to be 
< 2%. The uncertainties in the theoretical predictions are due to the choice of /i and PDF, 



and are approximately 40-50% with some dependence on M (see Fig. 12). 



B. Quark Compositeness 



The existence of three generations of quarks and leptons suggests that they may not 
be fundamental particles. For example, it has been proposed |L7j that they could be com- 
posed of "preons" which interact via a new strong interaction called metacolor. Below a 
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FIG. 12. The differences between the alternative predictions and the reference prediction 
(/j = 0.5Ef ax , CTEQ3M) of the inclusive dijet mass cross section (Eq. [Dj|) at ^/s = 1800 GeV 
for \rf ct \ < 1.0. Shown are the alternative predictions for the choices (a) [i = (0.25, 0.75, 1.0, 
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and (d) for a s = 0.110 — 0.122 using the CTEQ4A PDFs compared with the calculation using 
CTEQ4M, for which a s = 0.116. 
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characteristic energy scale A, the preons form metacolor singlets that are the quarks. The 
scale A characterizes both the strength of the preon coupling and the physical size of the 
composite state (A is defined so that g 2 /^ = 1). Limits are set assuming that all quarks are 
composite and A ^> a/I (where a/I is the center of mass energy of the colliding partons), so 
that quarks appear to be point-like. Hence, the substructure coupling can be approximated 
by a four-fermion contact interaction described by an effective Lagrangian |17| : 



C = &[VlL fe7^L)fe7^L) + 

Vlr {Ql^Ql) (qrI^r) + 
Vrl (QrI^Qr) (MlIpQl) + 

V°RR {QrYQr) {QrI^r) + 



Vll yiLT—qi I I QLl^—qL I + 
Vlr I <?L7 M y<?L j I ^^Y qR J + 
Vrl I <?r7 y<?/? I I QlI^Ql I + 

where rf^ H = 0, ±1, and H= L, R for left- or right-handed quarks, rj^j terms correspond 
to color-singlet (octet) contact interactions. These contact interactions modify the cross 
sections for quark-quark scattering. Limits are presented in Sections X.111 and XIV for the 

cases 

• a ll> wnere Vll = ±L 

• Ay, where rf LL = r)° RR = r]° RL = rf LR = ±1. 

• A % where Vll = Vrr = ~V°rl = ~Vlr = ±1- 

• Af y _ A) , where rj° LL = rj RR = 0; rj RL = rf LR = ±1. 

• A v 8 , wh ere rj\ L = r] RR = rj RL = r]\ R = ±1. 

• A^ g , where r] l LL = r] RR = Vrl = ~Vlr = ±1- 

• A fv-A) 8 , where Vll = Vrr = 0; Vrl = Vlr = ±1- 

Currently, there are no NLO compositeness calculations available; therefore LO calculations 
are used. The ratio of each LO prediction including compositeness to the LO prediction 
with no compositeness (A = oo) is used to scale the jetrad NLO prediction: 



a (A = X)t 

a (composite) = — t- \~ a ( = °°)nlo (4.18 

a (A = ooj LO 
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C. Coloron Limits 



A flavor-universal coloron model [Rq] inspired by technicolor has been proposed to explain 
the nominal excess in the inclusive jet cross section as measured by CDF [[F]. The model is 
minimal in its structure in that it involves the addition of one new interaction, one new scalar 
multiplet, and no new fermions. The QCD gauge group is extended to S77(3)i x 5 , C/(3) 2 . At 
low energies, due to symmetry breaking, this results in the existence of ordinary massless 
gluons and an octet of heavy coloron bosons. Below the mass of the colorons (M c ), coloron- 
exchange can be approximated by the effective four-fermion interaction: 

£ eff = —%mt [rfY q ) { q ^Y q ) (419) 

where cot 9 represents the mixing between colorons and gluons, and g\ = 4ira s - If 
Ay g y/c^ = M c j cot 9, this corresponds to Eq. |4.17| with r]\ L = r]\ R = T] RL = r)\ R = — 1 



and would represent new color-octet vector current-current interactions. Such interactions 
could arise from quark compositeness or from non-standard gluon interactions (e.g. gluon 
compositeness) [p9|. 



The phenomenology of the coloron has been studied [40] and limits have been placed on 



M r and cot 9. Constraints on the size of the radiative corrections of the weak-interaction 



p parameter require M c j cot 9 > 450 GeV [38], and a direct search for colorons in the dijet 



mass spectrum by the CDF Collaboration excludes colorons with mass below 1 TeV for 

cot9< i.5 grp. 
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V. TRIGGERING 

The D0 trigger was based on a multi-level system. The Level (L0) trigger consisted of 
two scintillating hodoscopes, one on each side of the interaction region. Coincident signals 
in the two hodoscopes indicate an inelastic collision and provide timing information for 
calculation of the position of the z- vertex of the interaction. The next level, Level 1 (LI), 
was a hardware trigger. The LI trigger required a specified number of calorimeter trigger 
tiles (At/ x A0 = 0.8 x 1.6) or towers (At/ x A0 = 0.2 x 0.2) above certain thresholds. 
Different trigger versions with slightly different LI requirements were instrumented during 
the run. If a LI rate was too large, a prescale was used to reduce the rate to an acceptable 
level. These prescale values were adjusted during the course of a beam store. A prescale 
of P allows only 1 out of every P events to be sent to the next level. Finally, Level 2 (L2) 
was a software trigger which selected the data to be written to tape. A fast jet algorithm 
used at L2 defines jet E^ as the sum of the transverse energy within a cone of opening angle 
1Z — 0.7 centered on the ^-weighted center of a LI trigger tile or tower. 

The L2 triggers used in the QCD analyses at yfs = 1800 GeV are called JET_30, JET_50, 
JET_85, and JET_115. The names follow the nomenclature that a JET_X trigger at L2 
requires at least one jet with E? greater than X GeV. During the running at y/s = 630 GeV 
the L2 triggers were JET_12, JET_2_12, and JET_30. A complete description of the LI and 
L2 trigger requirements is given in Table |. 



Trigger 


Level 1 (GeV) 


Level 2 (GeV) 


98% efficient 










= 1800 GeV 




JET.30 


1 tile 


> 


15 


1 jet with 


45 GeV 




& 1 tile 


> 


6 


E T > 30 




JET.50 


1 tile 


> 


15 


1 jet with 


75 GeV 




& 1 tile 


> 


6 


E T > 50 




JET.85 


1 tile 


> 


35 


1 jet with 


105 GeV 




& 2 tiles 


> 


6 


E T > 85 




JET.115 


1 tile 


> 


45 


1 jet with 


170 GeV 




& 1 tile 


> 


6 


E T > 115 










v^ 


= 630 GeV 




JET.12 


1 tower 


> 


2 


1 jet with 


20 GeV 










E T > 12 




JET_2_12 


2 towers 


> 


2 


1 jet with 


30 GeV 










E T > 12 




JET.30 


1 tile 


> 


15 


1 jet with 


45 GeV 










E T > 30 




TABLE I. Typical trigger configurations used in inclusive analyses. The LI and L2 require- 
ments are shown for each trigger. Also shown are the leading uncorrected jet Ei at which the 
average event trigger efficiency exceeds 98%. Redundant lower-Ii/p thresholds at LI were used to 
provide extended lists of seeds for jet clustering at L2. 



A study was performed to determine the trigger efficiency as a function of jet E^ for all 
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triggers used in D0 QCD analyses. There is an efficiency for an event to pass the LI trigger, 
and an efficiency for an event to pass the L2 trigger given that it passed LI. The combined 
efficiency to pass both LI and L2 is : 

total LI L2[L1 ^ 

c event c event c event V-"-- 1 -/ 

where e^in* is the efficiency for an event to pass L2 when it has passed LI. The LI and L2 
event efficiencies (e^ent an d Cevent) depend on the event topology (E^ and 77 of the jets in 
the event). The event trigger efficiency as a function of single jet efficiencies for an event 
with iVje ts is given by 

Wjets 

Cevent = 1" U [1 Tfc)] , (5-2) 

8=1 

where e» is the single jet efficiency for the ith jet. The product represents the probability 
that none of the jets in the event pass the trigger requirements. 

The efficiency of the L 1 trigger with the least restrictive requirements was measured using 
a data set that was required to pass only the L0 trigger. The single jet efficiency is given 
by the fraction of jets that satisfy the LI requirements at a given E^. The LI efficiencies for 
more restrictive LI triggers (MRT) were calculated using data samples that were required to 
pass a less restrictive LI trigger (LRT). This allows the LI efficiency of the more restrictive 
trigger to be calculated relative to the less restrictive trigger (given by 6mrt,lrt)- Hence the 
efficiency for a given LI trigger is given by the product of the efficiencies of all less restrictive 
triggers at a given E^. For example, the LI efficiency for the Jet_85 trigger is given by 

6 Jet_85 = 6 Jet_85,Jet_50 X 6 Jet_50,Jet_30 X 6 Jet_3O,L0- (5-3) 

The L2 trigger efficiencies for single jets are measured with respect to the LI trigger. 
The fraction of these events which have a L2 jet above threshold determines the L2 single 
jet efficiency. Figure [IB] shows the event efficiency for Jet_85 as a function of E^. 

Table | shows the typical trigger requirements and the E<? value for the leading jet 
at which each trigger averages an efficiency exceeding 98%. The leading jet's must be 
significantly higher than the L2 threshold in order for the trigger to be efficient. 
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1,1,1 






' 0.6 < |ri jet | < 1.6 

i,i,i 






1.6<|ri jet |<4.0 

i.i.i 



E T (GeV) 

FIG. 13. Average event efficiency for JET_85 as a function of leading jet Err and for three 
different pseudorapidity regions. 
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VI. LUMINOSITY 



The beam luminosity was calculated from the counting rate of the L0 counters and 
the cross section subtended by these counters. The cross section was determined using the 
geometric acceptance of the L0 hodoscopes, the L0 hardware efficiency, and the world 
average (WA) of the pp inelastic cross section measurements. The cross section of ob- 
served events in the L0 was found to be <7l0 = 43.1 ± 1.9 mb at \fs = 1800 GeV |42| and 
o"l0 = 32.9 ± 1.1 mb at y/s = 630 GeV (see j43| for a description of the method used). The 
effective luminosity was determined independently for each trigger on a run-by-run basis 
taking into account each trigger's prescale, the L0 inefficiency, and the detector deadtime. 
The WA pp cross section at \/s = 1800 GeV used in this paper is based on measurements by 
the E710 Collaboration JP|, the CDF Collaboration |E§, and the E811 Collaboration 
At y/s = 630 GeV there is no complete measurement of the pp cross section. Hence, the 
pp cross section was obtained by interpolating between the WA pp cross sections measured 
at y/s = 546 and 1800 GeV |43| . The WA pp cross section at yfs = 546 GeV is based on 



measurements by the UA4 [[47]] and CDF |45| Collaborations. The CDF collaboration only 
uses its measurement of the pp cross section to determine its luminosity. As a result the lu- 
minosities currently used by CDF are 6.1% lower than those used by D0. As a consequence 
of this all CDF cross sections are 6.1% lower than D0 cross sections |42]j . 

The integrated luminosities at y/s = 1800 GeV as measured using L0 for the Jet_30, 
Jet_50, Jet_85, and Jet_115 triggers are 0.368, 4.89, 56.7, and 95.7 pb -1 respectively, with 
an uncertainty of 5.1%. The luminosities at a/s = 630 GeV for the JET_12, JET_2_12, and, 
JET_30 were 5.12, 31.9 and 538 nb" 1 respectively with an uncertainty of 4.4%. 

The luminosity required corrections due to small discrepancies in the luminosity calcula- 
tion during different running periods at yfs = 1800 GeV. The initial luminosities for triggers 
Jet_85 and Jet_115 were taken from the luminosity calculation exclusively determined with 
the L0 counters. The inclusive jet cross sections calculated with the first 7.3 pb -1 of the 
data sample showed a 10% difference for Jet_115. The luminosity has been adjusted so that 
the dijet mass spectrum for the first 7.3 pb -1 matches that of the remaining data. This 
adjustment was also applied to Jet_85. Thus the luminosities for Jet_85 and Jet_115 are 56.5 
and 94.9 pb -1 respectively, a change of 0.7% from the value obtained using the L0 counters. 
This difference was added linearly to the 5.1% error on the initial luminosity value for a 
total error of 5.8%. 

In addition, for a part of the run the Jet_30 and Jet_50 triggers each required a sin- 
gle interaction at L0. The luminosities of Jet_30 and Jet_50 from the L0 calculation are 
estimated to be accurate only to 10% due to uncertainties in the efficiency of the single 
interaction requirement. The luminosity for the Jet_50 trigger was determined by matching 
the Jet_50 cross section to the Jet_85 cross section, and the Jet_30 cross section was matched 
to the Jet_50 cross section in regions of overlap. The trigger matching is analysis dependent; 
each analysis presented in this paper used the cross section of interest to match the triggers. 
The results obtained for the different measurements are consistent. The trigger matching 
errors are added in quadrature to the 5.8% error on Jet_85. The final luminosity and error 
for each trigger is shown in Table 



Since the analyses PJ25|| were first presented, the E811 Collaboration measurement of 
the total inelastic cross section was published |j46| . Including this measurement in the WA 
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Trigger 


Luminosity 


Error 


= 1800 GeV 


JET.30 


0.364 pb" 1 


7.8% 


JET.50 


4.84 pb- 1 


7.8% 


JET.85 


56.5 pb" 1 


5.8% 


JET.115 


94.9 pb" 1 


5.8% 


y/s = 630 GeV 


JET.12 


5.12 nb- 1 


4.4% 


JET_2_12 


31.9 nb" 1 


4.4% 


JET.30 


538. nb" 1 


4.4% 



TABLE II. Corrected luminosity and errors for the inclusive jet triggers. The trigger matching 
for Jet_30 and Jet_50 at \fs = 1800 GeV was carried out using the dijet mass cross section. 



changed the observed L0 cross section from o"l0 = 44.5 ± 2.4 mb to ctl0 — 43.1 ±1.9 mb at 
y/s = 1800 GeV. This changed the integrated luminosity of Jet_115 from 91.9 ± 5.6 pb -1 to 
94.9 ±4.7 pb -1 , an increase of 3.2%. Hence all cross sections at ^/s = 1800 GeV reported in 
this paper are reduced by 3.1% from the previously published results. It is worth noting that 
the inclusion of the E811 result had no perceptible impact on the cross section interpolation 
to 630 GeV. 

The luminosity calculation consists of three distinct ingredients: the geometric accep- 
tance of the L0 hodoscopes, the L0 hardware efficiency, and the pp inelastic cross section. 
The luminosity uncertainties are listed in Table |TJ. The largest contribution to the lumi- 
nosity uncertainty at y/s = 1800 GeV derives from the World Average (WA) pp total cross 
section. The pp cross section at \/s = 630 GeV was determined from a fit to the values at 
y/s = 546 and 1800 GeV p|. 



Source of 


Uncertainty (Percent) 




Uncertainty 


1800 GeV 


630 GeV 


World Average pp Cross Sections 


3.70 


2.75 


Hardware Efficiency 


2.32 


3.12 


Geometric Acceptance 


2.73 


1.51 


Time Dependencies 


0.70 


0.00 


All sources 


5.81 


4.43 



TABLE III. Uncertainties in the luminosity calculation excluding trigger matching. 
Two Monte Carlo minimum-bias event generators (mbr E3] and dtujet [JD]) were used 



to determine the geometric acceptance of the L0 hodoscopes. The difference in acceptance 
between the two MC results was taken as a source of systematic uncertainty for each y/s. 
The consistent behavior of each generator relative to the other between center-of-mass en- 
ergies indicates that the systematic uncertainty may be considered completely correlated. 
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(a) Total Cross Section 
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(b) Elastic Cross Section 



(c) Single Diffractive Cross Section 



500 1000 1500 2000 

Center of Mass Energy (GeV) 

FIG. 14. The three fits to the World Average pp cross sections. The stars depict the WA 
cross sections at y/s = 546 and 1800 GeV, and the closed square shows the interpolation to 
yj~s = 630 GeV. A fluctuation of the 1800 GeV point directly influences the interpolated value at 
630 GeV, particularly in the case of the total cross section (a). 
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Although the geometric acceptance of the L0 hodoscopes for diffract ive processes must 
be considered in luminosity calculations, the uncertainty in the non-diffractive acceptance 
dominates. 

A study of zero-bias events (a random sampling of the detector during a beam-beam 
crossing) determined the hardware efficiency of L0. Because the same estimation of the 
uncertainty appears in the calculation of the luminosities at both yfs values, the uncertain- 
ties are completely correlated. Table |TJ lists the systematic uncertainty in the hardware 
efficiency for both center-of-mass energies. 
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VII. THE EVENT VERTEX 



The location of the event vertex was determined using the central tracking system [28 



which provides charged particle tracking over the region | r\ \ < 3.2. It measures the trajectory 
of charged particles with a resolution of 2.5 mrad in and 28 mrad in 9. From these 
measurements the position of the interaction vertex along the beam direction (z) can be 
determined with a resolution of 8 mm. 

As the instantaneous luminosity increases, the average number of pp inelastic collisions 
per beam crossing increases. Hence there is the possibility of selecting the incorrect interac- 
tion vertex. If the incorrect vertex is chosen as the primary vertex, jet and event missing 
transverse energy ($t) wm be miscalculated. This may result in a significant contribution 
to the jet spectra at very high-Er since the high rate of jet production at lower Et can 
cause contamination in the lower rate regions. Visual scanning of the high-.&r jets shows 
that approximately 10% have misidentified interaction vertices. 



In order to study the effects of multiple interactions, a software tool called mitool foO 
was developed to provide information about the number of interactions. This tool uses the 
L0 hodoscopes, the calorimeter, and the central tracker in order to evaluate the number of 
interactions. A sum-of-times inconsistent with a single interaction from the L0 hodoscopes 
indicates the possibility of the presence of more than one interaction. The total energy in 
the calorimeter provides evidence of multiple interactions. If the total measured energy of an 
event is greater than 1.8 TeV, a multiple interaction is likely. Additional information from 
the number of vertices found with the central tracker is also used. Using this information 
the most probable number of pp interactions in the event is calculated. 

To a good approximation, the jet E T and r\ can be calculated for the second pp vertex 
using the measured vertex ^-position and a simple geometric conversion. Thus for all the jets 
in an event, the absolute magnitude of the vector sum of the jet Et, denoted St =\T,E 3 ^ | , can 
be calculated for each vertex. Except for soft radiation falling below the jet reconstruction 
threshold, St will be equal in magnitude to the $ T . Since QCD events should contain little 
$Ti the correct vertex was selected by choosing the vertex with the minimum St- 
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VIII. JET AND EVENT SELECTION 

The existence of random spurious energy deposits in the calorimeter may either fake 
or modify a real jet. Some sources of noise are electronic failures, cosmic ray showers, or 
accelerator losses due to Main Ring activity. A series of quality cuts was developed to remove 
this contamination. 

A. Removal of "Hot" Cells 

Before jet reconstruction, a cell suppression algorithm was implemented to suppress any 
cell with an unusually high deposition of energy relative to its longitudinal neighbors (a 
"hot" cell). Specifically, if a cell had more than 10 GeV of energy and more than 20 times 
the average energy of its immediate longitudinal neighbors, the cell energy was set to zero. 
This algorithm is successful in removing isolated high energy cells due to noise; however, 
the algorithm can also degrade the response to jets. 

Approximately 10% of the events have one or more suppressed cells. The rapidity distri- 
bution of the suppressed cells is very "jet-like" with a central plateau. A cell was restored 
to a jet if it was within /S.1Z = 0.7 of the original jet direction and if the cumulative total of 
hot cell Et was no more that 50% of the original jet E^. The jet rapidity and azimuth were 
then recalculated using the Snowmass definitions (Eq. |3~2j ). The event $ T was also adjusted 
if a cell was restored to a jet. 

The restoration algorithm has been shown to be 99% efficient by fitting the A1Z and 
restored cell fraction (the hot cell E^ divided by the jet's original E^) distributions and 
estimating the inefficiency in the cut regions. An event scan with restored jets (using relaxed 
restoration criteria) above 260 GeV showed no inefficiency Less than 5% of these new jets 
are contaminated. For those events with a single suppressed cell, the $ T is significantly 
reduced by the cell restoration. The kinematic variables (Ey, rj, and 0) of the high-Er jets 
which included a restored cell were compared to the kinematic variables calculated with a 
full reconstruction in which the suppression algorithm was disabled. The differences were 
small and well within the characteristic resolutions of the variables. 

B. Quality Cuts 

Even after the removal of isolated anomalously large cell energies, there still remain 
spurious jets. Quality cuts were developed to remove these fake jets. The quality cuts were 
applied on either the jet or to the event. 

The jet quality cuts are based on the distribution of energy within the jet. Three standard 
variables are used: 

1) Electromagnetic Fraction (EMF) — the fraction of the jet energy contained in the 
electromagnetic section of the calorimeter. Jets are retained if: 

EMF<0.95 (1.2 <|%et|< 1-6), 
0.05 < EMF< 0.95 (otherwise), (8.1) 
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where r] det is the pseudorapidity of the jet calculated using a vertex position of z = 0. 
The cut EMF > 0.05 is not applied for 1.2 <| r] dct |< 1.6 because of the gap between 
the CC and EC calorimeters (Section ||). 

2) Coarse Hadronic Fraction (CHF) — the fraction of the jet energy contained in the 
coarse hadronic section of the calorimeter. This cut is designed to remove fake jets 
introduced by Main Ring particles depositing energy in the calorimeter. Jets are 
retained if: 

CHF < 0.4. (8.2) 



3) Hot Cell Fraction (HCF) — the ratio of the most energetic cell of a jet to the second 
most energetic cell. Jets are retained if: 

HCF < 10.0. (8.3) 



A cut on $ T is also used to remove bad events. Since QCD events are expected to have 
no intrinsic $ T , a cut on events with large $ T typically used: 



T 



E 



jetl 



< 0.7, (8.4) 



Jrp 

where E^ 1 is the transverse energy of the highest Et jet in the event. In the case of the 
inclusive jet analysis, the measurement is more susceptible to contamination from events in 
which the primary vertex is located outside of the tracking detector and a vertex due to an 
additional minimum bias event is identified as the primary vertex, leading to an overestimate 
of the jet .Ex- In this case it is found that a $ T cut of 

E 

^Zr < 0.3 if E£ tl > 100 GeV, 



E T 



jetl ^ u u 11 
$ T < 30 GeV if E^ 1 < 100 GeV (8.5) 



removes the contamination. 

Since the data collected at ^Js = 630 GeV were taken at low instantaneous luminosity, 
there were fewer events with multiple interactions and incorrectly identified vertices. As a 
result, the quality cuts on the data were adjusted to maximize efficiency without increasing 
contamination. The resulting cuts are: 

EMF < 0.90, (8.6) 
CHF < 0.4, (8.7) 
HCF < 20.0, (8.8) 
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C. Efficiency 



The efficiencies of the quality cuts were measured. The data sample used to calculate 
the efficiencies was selected by making cuts in 77 and <fi. We verified that the changes in 
shape of the EMF, HCF, and CHF distributions due to the $ T were negligible. 

To calculate the total efficiency, each individual efficiency is measured. First the $ T cut 



is applied to the data and the efficiency of the EMF cut is calculated. Figure [15| illustrates 
the EMF distribution after the $ T cut is applied. Contamination is visible as small peaks 
near EMF m or 1. A Gaussian-like curve is projected under the noise signal and used to 
estimate the data signal lost due to the cut. After the $ T and EMF cuts are applied, the 
HCF efficiency is measured (Fig. 0(a) shows the HCF distribution). Then after both the 
EMF and HCF cuts were made, the efficiency of the CHF cut is measured. The total jet 
efficiency is calculated by multiplying the individual cut efficiencies together. 

The standard jet cuts remove most of the noise from the sample; however, there is still 
some contamination at high due to cosmic rays and "Main Ring Events". The $ T cut 
removes this remaining noise. By fitting this distribution in regions of the calorimeter where 
the noise effects are negligible, an extrapolation can be used to determine the efficiency. 
Figure |l6|(b) shows the $ T distribution used to calculate the efficiency. The inclusive jet 
efficiencies as a function of jet in the central region at ^/s = 1800 GeV are shown in 
Fig. |17[ Figure shows the efficiency of the $ T cut for dijet events. The efficiencies of the 
quality cuts used at yfs = 630 GeV are given in Fig. [19]. The efficiency of the $ T cut at 
fs = 630 GeV is > 99%. 



a 

> 




EMF Distribution 

FIG. 15. The measured EMF distributions for different ranges. The lower plots show the 
cut values and the fit used to calculate the efficiency of the cut. The dashed histogram shows the 
full data sample and the solid histogram shows a data sample with minimal noise contamination. 
The arrows indicate the cut values. The peaks at EMF ~ or 1 are due to contamination. 
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FIG. 16. a) The 1/HCF distribution. The arrow shows the location of the cut for 
y/a = 1800 GeV. b) The distribution of £^ etl /#r- The arrow at E^ tl = 1.43 corresponds 
to the ^j'/E^ tl cut of 0.7. The peak at 1.0 is due to contamination from cosmic rays and the Main 
Ring. The dashed histograms show the distributions for the full data samples. 
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FIG. 17. Top: The efficiency of the standard jet quality cuts for 1 77 1 < 0.5 (Eqs. |8.1| , |3^ , and 
at yfs = 1800 GeV. Bottom: The efficiency of the $t cu t used in the inclusive jet analysis 
(Eq. |8.5p at y/s = 1800 GeV. The dotted curves show fits to the measured efficiencies. 
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FIG. 18. The efficiency of the $ T cut used in the dijet analyses at y/s = 1800 GeV (Eq. |0) 




20 30 40 50 

E T (GeV) 

FIG. 19. The efficiency of the standard jet quality cuts for | r/ 1 < 0.8 at yfs = 630 GeV (Eqs. 
8.7, and 3.8). The three curves show the fit to the efficiencies and the uncertainty in the fit. 
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D. Contamination 

In order to measure the remaining contamination after all quality cuts have been im- 
plemented, two separate studies were performed. Residual contamination was estimated 
by overlapping the observed hot cell distribution on a simulated inclusive jet sample. The 
simulated cross section changed by less that 1% after imposition of the jet quality cuts. 
The simulation also indicated that the jet quality cuts reject > 99% of the "fake" jets with 
E T = 500 GeV. 

To measure the contamination due to misvertexing, events at high-i?T were visually 
inspected. Misvertexing tends to cause lower E T jets to migrate to higher E T . Since the 
cross section is steeply falling, this can corrupt the high-.&r cross section. This study shows 
that after the vertex selection procedure has been applied, less than 1% of the events are 
contaminated at high E T . 
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IX. JET ENERGY SCALE 



The in-situ jet energy calibration uses reconstructed collider data, and is described in 
more detail in |H| . The measured energy of a jet -Ej™ eas depends strongly on the jet definition. 



The particle-level (true) jet energy Ef^ is defined as the energy of a jet consisting of 
final-state particles produced by the high-p T parton-parton scattering, and found using the 
Snowmass algorithm. The jet should not include the particles produced by the underlying 
event (Section |III C|) . The jet energy scale corrects the measured jet energy, on average, 
back to the energy of the final-state particle-level jet. E^ 1 is determined as: 



pmeas zp 

< = % q ° (9.1) 



where: 



• Eq is an offset, which includes the underlying event, noise from radioactive decays in 
the uranium absorber, the effects of previous interactions (pile-up), and the contribu- 
tion from additional pp interactions in the event. 

• -Rj et is the calorimeter energy response to jets. i?j et is typically less than unity due 
to energy deposited in uninstrumented regions of the detector, and differences in the 
response to electromagnetic and hadronic particles (e/h > 1). 

• Sh is the fraction of the jet energy that showered inside the algorithm cone at the 
calorimeter level. 

The calibration is performed using data taken in pp collisions at yfs =1800 GeV and 630 GeV. 

A. Offset Correction, Eq 

The total offset correction is measured as a transverse energy density in rj-cf) space and 
subdivided as Do = D ue + D®. D ue represents the contribution due to the underlying event, 
i.e. energy associated with the spectator partons in a pp event. D® accounts for uranium 
noise, pile-up, and energy from additional pp interactions. The offset correction Eo is given 
by Do multiplied by the r]-<f> area of the jet. 

D ue is measured using the average transverse energy density in minimum-bias events 
(where a pp interaction has occurred, usually inelastic scattering). D® is determined from a 
zero-bias sample (a random sampling of the detector during a beam-beam crossing). The rj 



dependencies of both quantities and the luminosity dependence of D® are shown in Figs. £0 
and The statistical and systematic errors of the offset correction are 8% and 0.25 GeV 
respectively. 



B. Response Correction, i?j et 



D0 makes a direct measurement of the jet energy response using conservation of Pt in 
photon-jet (7-jet) events [ 5T| . The electromagnetic energy scale is determined from the D0 
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FIG. 20. Physics underlying event density D ue versus r\ for events with y/s = 1.8 TeV and 
fs = 630 GeV. 
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Z{—* e + e ), J/if>, and 7r° data samples, using the masses of these known resonances. In the 
case of a 7-jet two body process, the jet response can be measured through: 

= 1 + %^ , (9-2) 

where and n^ 7 are the transverse energy and direction of the photon. To avoid response 
and trigger biases, i?j et is binned in terms of E' = E Tl x cosh (77^) and then mapped onto 
i?™ t eas . E' depends only on photon variables and jet pseudorapidity, which are both measured 
with very good resolution. 



1. rj-Dependent Corrections 



Most measurements need a high degree of accuracy in the jet energy scale at all rapidities. 
An ^-dependent correction becomes necessary. The cryostat factor F CTy is defined as the ratio 
Rf^/R^. The measured factor 0.977 ± 0.005 is constant as a function of E' . This was 
expected because the CC and the EC calorimeters use the same technology. 

The intercryostat region (IC), which covers the pseudorapidity range 0.8 < \ r} \ < 1.6, is 
the least well-instrumented region of the calorimeter system. A substantial amount of energy 
is lost in the cryostat walls, module end plates, and support structures. An IC correction 
is performed after the F cry correction and before the energy-dependent response correction. 
Because the energy dependence of R ]et is included in R )et as a function of rj, this function 
is not a constant, but should be smooth. The IC correction is set so that the response as a 
function of rj agrees with the fit to the functional form, Rj Ct = a + b ■ In [cosh (77)], of the CC 



and EC response, as shown in Fig. 22 
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FIG. 22. Response versus rj for 7-jet data before the ry-dependent correction. The dashed line 
is the fit to the expected IC response. 
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2. Energy- Dependent Correction 



Following the above procedure, the energy dependence of _Rj e t is then determined as a 
function of E' as illustrated in Fig. Uniformity of the calorimeters allows the use of data 
from both the CC and the EC to measure Rj et . The rapidly falling photon cross sections 
limit the use of CC data to energies <120 GeV. EC data are used to extend the energy 
reach to ~300 GeV. We exploit the fact that jet energy in the EC is larger than in the CC 
for the same Et- Monte Carlo data are also included at the highest energy to constrain 
the extrapolation. A set of 7-jet events is generated using herwig |33|], processed through 
the D0GEANT [|34| detector simulation, and reconstructed with the standard photon and jet 
algorithms. The Monte Carlo simulation is improved by incorporating the single particle 
response of the calorimeter as measured in test beam. 

The response versus energy for the 1Z = 0.7 cone algorithm is shown in Fig. |2^. The CC 
and EC data, and the expected response from MC at E^ = 500 GeV are fit with the func- 
tional form Rj et (E) = a + b ■ \n(E) + c • ln(E) 2 (see Fig. |24|). This function is motivated by 
the hadronic shower becoming gradually more "electromagnetic" (EM) with increasing en- 
ergy [52]. If e and h are the responses of the calorimeter to the EM and non-EM components 
of a hadronic shower, and tt is the response to charged pions, then e/ir = l/[~ — (/em)(- — 1)]- 
The functional form for the mean electromagnetic fraction of the jet (Jem) is ~ a • ln(E), 
giving the expected logarithmic dependence for energy carried by the charged pions and, 
therefore, jets. 
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FIG. 23. Rj c t versus E' measured in the CC, IC and EC calorimeter regions after r\ dependent 
corrections. 

In addition to the uncertainty from the fit (1.5%, 0.5%, 1.6% for 20, 100, 450 GeV jets 
respectively), there is also a ~0.5% uncertainty from the W boson background in the photon 
sample. Some of the events in the 7-jet sample are not two-body processes. In the IC region, 
the ^-dependent corrections contribute an additional ~1% uncertainty. 
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C. Showering Correction, Sh 

As a jet of particles strikes the detector, it interacts with the calorimeter material pro- 
ducing a wide shower of particles. Some particles directed inside the cone deposit a fraction 
of their energy outside the cone (and vice versa) as the shower develops inside the calorime- 
ter. We do not correct for any QCD radiation or particles that are radiated from the cone; 
we only correct for the effects of the detector. 

The correction for this showering is determined using jet energy density profiles from 
data and particle-level herwig |J3| Monte Carlo. The data contains the contributions of 
both gluon radiation and showering effects outside the cone. The former contribution is 
subtracted using the particle-level Monte Carlo profiles. Sh is defined as the inverse of the 
measured correction factor; that means Sh is the fraction of the jet energy showered inside 
the algorithm cone in the calorimeter (Eq. |9.1|) . The showering correction is negligible for 
7Z = 0.7 cone jets above ~100 GeV in the central region (\r)\ < 1.0) with an uncertainty 
of ~ 1%. Both the correction and uncertainty are larger for lower energies, higher t], and 
smaller cone sizes. 



D. Correlations of the Uncertainties 

The uncertainties in the jet energy scale can be separated into five sources: offset, in- 
dependent corrections, response corrections, method, and showering corrections. The corre- 
lations of these uncertainties as a function of Et and r\ have been studied: 

1. Offset: This is the dominant uncertainty at low Et but is unimportant at high Et- 
The uncertainty due to the offset correction is divided into two parts: a systematic 
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error related to uncertainties in the method which is correlated as a function of E^, 
and a statistical error that is uncorrelated as a function of pseudorapidity due to the 
finite size of the data sample used to determine the offset. 

2. 77-Dependent Correction: The uncertainty due to this correction was separated 
into two parts. The first is due to the cryostat factor and is correlated as a function 
of Et and rj. The second is the IC correction, which is uncorrelated as a function of 
Et and rj. 

3. Response Correction: The uncertainty associated with the hadronic response is 
unimportant at low Ey but dominant at high E^. As a result of using a fit, the 
uncertainty is partially correlated as a function of E^. The correlation matrices for 
various jet cone sizes can be found in Ref. [[Bill . 

4. Method: The uncertainty in the method used to determine the energy scale correction 
arises from the data selection requirements, and punch-through at very high energies. 
The method uncertainty is correlated as a function of E^. 

5. Showering Correction: The uncertainty due to this correction is small except at 
very low E T and is considered to be fully correlated as a function of E T . 



E. Summary and Verification Studies 



Figure ^ shows the magnitude of the correction and uncertainty for 7Z — 0.7 cone jets 
with rj =0. The overall correction factor to jet energy in the central calorimeter is 1.160 ± 
0.018 and 1.120 ± 0.025 at 70 GeV and 400 GeV, respectively. Point-to-point correlations 
in the energy uncertainty are very high for jets with 200 < E^ < 450 GeV. 

The accuracy of the jet energy scale correction is verified using a herwig 7-jet sample 
and the D0GEANT detector simulation. A Monte Carlo jet energy scale is derived and 
the corrected jet energy is compared directly to the energy of the associated particle jet. 
Figure |26| shows the ratio of calorimeter to particle jet energy before (open circles) and after 
(full circles) the jet scale correction in the CC. The ratio is consistent with unity to within 
~ 0.5%. 
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sample, and the outer error bars represent the systematic uncertainty on the energy scale. 
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X. JET RESOLUTIONS 



The observed energy distributions are smeared due to resolution effects. The fractional 
energy resolution ue/E may be parameterized as: 



The nature of the incident particles, sampling fluctuations, and showering fluctuations, 
contribute mostly to the sampling term, S. Detector imperfections and deviations from an 
electron/hadron single particle response of unity, limit the resolution at high energies and 
are described by the constant term, C. Noise fluctuations (including the effects of multiple 
interactions) affect the low energy range and are given by the noise term, N. 

In the analyses reported here, we measure the E T of the jets; hence we need to measure 
the resolution of Et, which will have the same form: 



The relationship between oeIE and oe t IE^ depends on the 77 resolution, a v . Using E^ = 
E I cosh 77 and assuming that <je t and <r v are uncorrelated then: 



In addition to the detector resolution, other contributions must be folded into the res- 
olutions used for physics analyses. These are, for example, fluctuations of the out-of-cone 
losses, and the fluctuations of the vertex z-position about its measured value. 

Using D0 dijet data we derive the energy resolutions using energy conservation in the 
transverse plane. The following criteria are applied to dijet events in order to eliminate 
sources of contamination due to additional low-£r jets: 

• The z-coordinate of the interaction vertex must be within 50 cm of the center of the 
detector. 

• The two leading- E T jets must be back-to-back (A0 > 175°). 

• If there is a third jet in the event, it must have E^ t3 less than a specified value. 

• All jets in the event must satisfy the jet quality cuts. 

• Both leading jets are required to be in the same 77 region so that their resolutions are 
similar, i.e. |77j etl | |r]j et2 |. 

The dijet balance method is based on the asymmetry variable A, which is defined as: 




(10.1) 




(10.2) 
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where E^ tl and E^ t2 are the randomly ordered transverse energies of the two leading-.&r 
jets in an event. The variance of the asymmetry distribution can be written as: 



dA 
dE$ tl 



U jetl + 



dA 
dE^ 2 



2 

U p jet2 . 



(10.5) 



Assuming Et = E^ tl = E^ z and <je t = <? E icti = a E iet2, the fractional transverse energy can 
be expressed as a function of a a in the following way: 



njct2 



<?E 7 
Enp 



y/2, 



OA- 



;io.6) 



Figure [27] shows the asymmetry distributions A for different Et bins. The asymmetry 
distributions show minimal tails (-C 1%) and are well-described by a Gaussian distribution. 
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FIG. 27. Asymmetry distribution in several Et bins for jets with \ rj \ < 0.5 and E^ t3 < 8 GeV. 



A. Soft Radiation Correction 



Although the A(f> and third-jet E^ cuts (A0 > 175° and E^ t3 < 8 GeV) are designed to 
remove events with more than two reconstructed jets, events may still contain soft radiation 
that prevents the two leading-E^ jets from balancing in the transverse plane; therefore the 
measured resolutions are overestimates of the hypothetical "true resolutions." To evaluate 
this bias, the resolutions were determined from samples with different E^ t3 cuts: 8, 10, 12, 
15, and 20 GeV. The resolutions are then extrapolated to a "true" dijet system with E^ t3 



0. Figure shows the fractional jet resolutions as a function of E^ t3 cut for several E^ 



bins. 
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FIG. 28. Resolutions as a function of the cut on Ejf° for different E^ bins 



The 

= o, 



This procedure is repeated for every bin. We expect the correction for additional 
radiation in the event to be continuous as a function of jet E^ and to be given by a function 
K(Et). Because the soft radiation bias should primarily affect small values of E^ but be 
negligible at high Et, we parameterize the soft radiation correction K(E^) with the function: 



K(E 1 



1 — exp (— a — cliEt) 



;io.7) 



For the pseudorapidity bin \ rj \ < 0.5, a = 2.20 and a\ = 0.0055 (Fig. |29|). This parameter- 
ization corrects the resolutions of each E^ bin for the effects of soft radiation. 



Note that the point-to-point correlations in Fig. are very large because each data 
point represents a subsample of the data point to its immediate right. In addition, it is 
not clear that the linear trend continues down to E^ t3 = 0; hence we do not use the errors 
obtained from the fits to calculate the error on the corrected resolutions. The uncertainty in 
the extrapolation is the sum in quadrature of the following: the uncertainty in the resolution 
at E^ t3 > 8 GeV, the difference in the extrapolation to E^ t3 = including and excluding 
the sample with the E^ t3 cut of 8 GeV, and the uncertainty in the fit to the point-to-point 
correlations. 



B. Particle Jet Imbalance 

Since we are correcting our measurements to the particle level, we must not include the 
effects of hadronization of the quarks and gluons in the resolutions. The energy carried 
by particles emitted outside the particle-level cone does not belong to the particle jet. In 
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FIG. 29. The soft radiation correction, K(E^), as a function of (\rj\ < 0.5). The error 
bars show the total uncertainty in the point-to-point correlations. The inner error bars show the 
uncertainty in the resolutions measured with E^ t3 > 8GeV. 



other words, at LO the total pt of a dijet event at the particle level is zero, but the two 
reconstructed particle jets do not necessarily balance, since there could be particles emitted 
outside the cones. The asymmetry distribution measures the detector resolution convoluted 
with the contribution of the dijet imbalance at the particle level. The latter must be removed. 

The particle-level resolution is obtained by applying the same techniques as used on the 
data to a herwig |33| Monte Carlo sample, e.g. no energy fluctuations. The calorimeter 
resolution is obtained by removing the particle-level resolution using: 
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The fractional E^ resolutions before the particle jet imbalance correction are shown in 
Fig. along with the MC data used to calculate the particle jet imbalance correction. The 
fully-corrected resolutions are given in Table [TV]. 



C. Studies of Systematic Uncertainties 



In principle, the soft radiation correction should remove the effects of additional gluon 
radiation in the data sample; however, this may not be the case because not all the particles 
present in the detector appear in reconstructed jets. It is also possible that the requirement 
that jets be back-to-back (A</> > 175°) preferentially selects events with better-than-average 
resolution. The possible size of these effects is studied by changing the back-to-back require- 
ment to A0 > 165° and repeating the determination of the resolutions. The result of this 
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FIG. 30. 

a E^/E T as a function of average E? for 1 77 1 < 0.5. The data points (squares) indicate 
the resolutions after the soft radiation correction and the solid curve shows the fit to the resolutions. 
The dash-dot lines show the systematic uncertainty due to the method. The dashed line is a fit to 
the particle-level resolutions obtained from MC points (circles). 



study is shown in Fig. The resulting resolutions are slightly higher than the resolutions 
calculated with a cut of A<ft > 175° and this difference is included in the overall systematic 
error. 

Some analyses require a tighter cut on the ]p T than the standard cut. In particular, the 
measurement of the inclusive jet cross section requires a $ T cut of $ T /E^ t:L < 0.3 when 
E^ n > 100 GeV, or $ T < 30 GeV when E^ 1 < 100 GeV. Any strengthening of the $ T 
requirement will implicitly reduce the difference between the E^'s of the two jets selected 
and also reduce the amount of soft gluon radiation; hence the resolutions should improve. 
The resolution parameterization using this $ T cut is depicted in Fig. 

The fractional Et resolutions are parameterized using Eq. 10.2 for all rapidities (| rj | < 1) 
and are given in Table |V| and are plotted in Figs. |3T] and [32|. 



D. Jet Resolutions at yfs = 630 GeV 

The jet resolutions at <Js = 630 GeV are measured using the same techniques as the 
resolutions at ^Js = 1800 GeV. These resolutions are supplemented at low values of jet E T 
by resolutions measured using photon-jet events. 

The energy resolution for photons is approximately 10 times better than that for a jet, 
allowing a convenient redefinition of Eq. 10. 4| . The photon-jet asymmetry is defined as 



A 7jet = T ^ , (10.9) 
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FIG. 31. Fully corrected oe t /Et as a function of average Et for 1 77 1 < 0.5 (i.e. the soft 
radiation correction and the particle- level dijet imbalance corrections have been applied). The data 
points (solid curve) show the resolution as calculated with cuts I^t/E^ 1 < 0.7 and A</> > 175°. 
The dashed line shows the effect of using a cut of Acft > 165°. In addition, the effects of using a $ T 
cut of $ T /E^ tl < 0.3 when E^ tl > 100 GeV, or $ T < 30 GeV when E^ 1 < 100 GeV are shown 
(dash-dot and solid-dots lines). 




57 



TABLE IV. The measured jet resolutions at y/s = 1800 GeV and their uncertainties. 
(£t> (GeV) a(E T )/E T A(a(E T )/E T ) 
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n 1 9n 

u. 1ZU 


u.uu^t 




n mfi 

U. 1UD 


U.UUo 


fi7 7n 

D 1 . 1 U 


u.uyo 


U.UUo 


86.43 


0.088 


0.001 


105.08 


0.078 


0.002 


130.42 


0.070 


0.001 


155.54 


0.068 


0.001 


182.40 


0.062 


0.001 


213.44 


0.056 


0.002 


241.69 


0.059 


0.003 


295.10 


0.050 


0.003 



TABLE V. The resolution fit parameters at y/s = 1800 GeV. 





n 


Fit Variables for a $ T cut of $ T / E^ 1 


< 0.7. 


C 




N 


\v\ 

0.5 < 
0.1 < 


< 0.5 

7/| < 1.0 
77 | < 0.7 


0.033 ± 0.006 
0.047 ± 0.008 
0.040 ± 0.013 


0.686 ± 0.065 
0.783 ± 0.137 
0.641 ± 0.160 


2.621 ± 0.810 
0.590 ± 9.334 
2.891 ± 1.413 




V 


Fit Variables for a $t cu * °f ^t/^t" 


< 0.3. 


C 


S 


N 


\v\ 


< 0.5 


0.037 ± 0.002 


0.514 ± 0.027 


4.009 ± 0.202 


0.5 < 


7/| < 1.0 


0.036 ± 0.006 


0.736 ± 0.059 


1.972 ± 0.904 


0.1 < 


7/| < 0.7 


0.038 ± 0.005 


0.550 ± 0.074 


3.654 ± 0.487 



where and are the fully corrected photon and jet transverse energies, respectively. 
If one approximates E^ ~ Ej? = E^ as before, and takes 8Ej~ ~ 0, the standard deviation 
of the photon-jet asymmetry identically becomes the fractional jet resolution: 



&E 7 



E'Y 

Figure I33J displays a typical distribution of photon-jet asymmetry. 



^ef (10-10) 



As described in previous sections, the measured resolution is adjusted to reflect third-jet 
biases and the particle-jet asymmetry. The results bolster the low-statistics dijet results at 
y/s = 630 GeV. The resulting resolutions are given in Table |VT| and are compared with the 
resolutions at y/s = 1800 GeV in Fig. |34[ It is clear that the resolutions at the different 
center-of-mass energies are significantly different (a probability of agreement of 0.0007). 
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FIG. 33. Distribution of photon-jet asymmetry for jet E^ between 15 and 20 GeV in the central 
region. 



TABLE VI. The measured jet resolutions (and uncertainties) at i/s = 630 GeV. 



Data Set 


(Et) (GeV) 


a(E T )/E T 


A(a(E T )/E T ) 


7- jet 


13.51 


0.205 


0.023 


7-jet 


17.81 


0.217 


0.048 


7-jet 


21.52 


0.175 


0.016 


7-jet 


24.27 


0.169 


0.019 


jet-jet 


26.28 


0.148 


0.012 


jet-jet 


34.35 


0.117 


0.015 


jet-jet 


40.87 


0.114 


0.010 


jet-jet 


52.27 


0.097 


0.009 


jet-jet 


59.12 


0.079 


0.007 


jet-jet 


70.53 


0.075 


0.006 
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FIG. 34. The single jet resolutions at yfs = 630 GeV (triangles) and 1800 GeV (circles). The 
resolutions at the two center-of-mass energies have been fitted separately to Eq. |10.2| . The fit to 
the yfs = 1800 GeV data is the solid line, and the fit to the = 630 GeV data is the dashed line. 

1. Parameterization of the Jet Resolutions. 

There are several parameterization choices that can be used to fit the data at both 
center-of-mass energies. We considered five alternative parameterizations of the resolutions: 



1. Fit the data simultaneously with Eq. |10.2[ the CSN model. 



2. Fit the data with common C and S terms and different noise terms (Ni 800 , -^630) a ^ 
the two CM energies: the CSNN model. 

3. Fit the data with common C and N terms and different sampling terms (Sisoo, 563o) 
at the two CM energies: the CSSN model (Fig. |35|) . 

4. Fit the data with a common C term and different sampling and noise terms at the 
two CM energies: the CSSNN model. 

5. Fit the data with no common terms: the CCSSNN model. 

A model where only the C term was allowed to vary between the two CM energies was not 
considered because C depends on the physical structure of the calorimeter, and hence should 
not change. The \ 2 and numbers of degrees of freedom for these five models are calculated 



and compared in Table [VII| . The fit parameters are given in Table [VIII . 

It is clear from the x 2 °f t ne parameterizations that the data cannot be represented by a 
single fit with common C, S, and iV (CSN model). Of the other models, the CSSN model 
gives the best fit to the data. If we allow additional parameters to be included in the fit, 
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FIG. 35. The single jet resolutions at y/s = 630 GeV (triangles) and 1800 GeV (circles). The 
resolutions at the two center of mass energies have been fitted using the CSSN model (solid lines). 



TABLE VII. 
lutions. 



X 2 for the different models that can be used to parameterize the single jet reso- 



Model 
CSN 
CSNN 
CSSN 
CSSNN 
CCSSNN 



X 
44.9 
25.5 
18.7 
17.9 
17.9 



Degrees of freedom 
19 
18 
18 
17 
16 



Probability 
0.0007 
0.11 
0.41 
0.35 
0.33 
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TABLE VIII. The fit parameters for all models used to fit the resolution data. The correlation 
matrix for the CSSN model is also given. 



Model 


Parameter 


Value 


Statistical Error 




N 


1.098 


1.128 


CSN 


S 


0.745 


0.038 




C 


0.028 


0.004 




N 


2.571 


0.309 


CSSN 


"51800 


0.691 


0.027 




S630 


0.510 


0.057 




c 


0.032 


0.003 


Correlation Matrix 


1.000 


-0.812 


-0.838 


0.575 


-0.812 


1.000 


0.751 


-0.902 


-0.838 


0.751 


1.000 


-0.589 


0.575 


-0.902 


-0.589 


1.000 




^1800 


3.543 


0.399 


CSNN 


^630 


1.907 


0.437 




S 


0.590 


0.049 




c 


0.040 


0.003 




^1800 


2.510 


0.893 




^630 


2.587 


0.374 


CSSNN 


Sl800 


0.696 


0.068 




£630 


0.509 


0.063 




c 


0.031 


0.007 
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the x 2 does not improve significantly. The CSNN model does not fit the data as well. The 
noise distribution in the calorimeter is similar at the two different CM energies (Fig. |36|) ; 
hence the CSSN parameterization model was chosen to fit the resolutions. The cause of the 
change in sampling term as a function of CM energy is not known. One possible explanation 
is that a jet sample at y/s = 630 GeV will have relatively more "quark" jets than a sample 
at yfs = 1800 GeV with the same mean Ei [f22],[53|1 . 
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FIG. 36. a) The average value and RMS width of the calorimeter noise distributions is given 
by the two lower points. The two upper points are the values of the N parameter obtained in the 
fit to the resolution data using the CSNN model, b) and c): The noise distribution found within 
a standard jet cone at each center-of-mass energy (the cone used to measure the noise is required 
to be at least 90° in <ft from any other jet in the event). 



E. Monte Carlo Consistency Tests 

To verify the resolution extraction methods, a Monte Carlo study compared events with 
and without the detector simulation. The jet resolutions of the MC sample are measured in 
two ways; the first is the asymmetry method, and the second is a direct measurement of the 
resolutions. If the Et of a jet as measured by the calorimeter is simply denoted by Et, and 
the Et as measured at the particle-level is denoted by E T tcl , then the jet resolution can be 
derived from the ratio: 



E* td - Eh 



Erp 



ptcl 



(10.11) 



Figure 37 shows the differences in the resolutions as measured by the two methods. The 



differences between the two are scattered about zero, indicating lack of bias in the method. 
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The differences between the two methods, less than 1%, indicates the magnitude of the 
systematic uncertainty, which can be parameterized as 

A (^) = i? +a0L (1012) 



C 05 

<u n 

a 

Q -0.05 



Vs = 1800 GeV 



100 



400 



500 



200 300 

Average E T (GeV) 

FIG. 37. Resolution closure from herwig Monte Carlo simulation; the difference in resolution 
obtained using the two techniques. The degree of closure is within 1% for all data points above 25 
GeV. 

Subsequent to the publication by D0 of the inclusive jet cross section || (Section [XI]) 
and the dijet mass spectrum [^] (Section |XIV|) at \fs = 1800 GeV the MC closure data for 
the resolutions were reexamined (see Fig. |38|). As a result of this, the MC closure error at 
y/s = 1800 GeV was reduced for E T > 40 GeV: 



A 
A 



E t / ^=630 GeV 
V Et J v^=i800 GeV 



2.23 

—2 + 0.0021. 

E'y 

14.1 

+ 0.0024. 



Et^ 



10.13) 



The effect of reducing the error on the inclusive jet cross section and dijet mass spectrum 
was negligible, and hence the results were not updated. The reduced errors are important 
for the analysis of the ratio of inclusive jet cross sections at a/s = 630 and 1800 GeV |2~I 
(Section pOl. 



In Fig. the measured resolutions are compared with the CSSN fit. The shaded region 



shows the size of the fit uncertainty, and the hatched region shows the size of the fit and 
MC closure uncertainties added in quadrature. Also shown are the other models. It is clear 
that the combined fit and MC closure uncertainties are of reasonable size and that the total 
uncertainties are not underestimated. 



F. ij and Resolutions 



After the 77-bias correction is applied, the average rj of the reconstructed jet is equal to 
the rj of a particle-level jet, but due to calorimeter showering effects, both i] and resolutions 
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FIG. 38. The improved resolution closure obtained using the herwig Monte Carlo simulation, 
for both center-of-mass energies. For most of the kinematic range, the degree of closure lies within 
a fraction of a percent. 



remain non-zero. The 77 resolution is obtained by using herwig Monte Carlo and studying 
^ptci — V as a function of jet energy and 77. Figure [H] shows the ^-resolution as a function of 
jet energy for different energy regions. The distributions show no tails and are well-described 
by a Gaussians. The resolution is determined by measuring ptcl — as a function of jet 
energy and 77. Figure f|l] shows the resolutions which are similar in magnitude to the 77 
resolutions. 
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FIG. 39. A comparison of the measured jet resolutions and the fit using the CSSN model. 
Also shown are curves representing comparisons between the different models and CSSN. The 
shaded regions show the uncertainty in the fit. The hatched region shows the magnitude of the fit 
and MC closure uncertainties added in quadrature. 
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FIG. 41. (p resolution as a function of the particle-level jet energy using a herwig simulation. 
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XI. INCLUSIVE JET CROSS SECTION AT y/s = 1800 GeV 

In this section we describe the measurement of the inclusive jet cross section in the 
pseudorapidity ranges | r\ | < 0.5 and 0.1 < | rj \ < 0.7. The inclusive jet cross section is given 
by: 

d 2 a Nid 

= — (11.1) 

dE T dr] C^AEtAt] v ; 

where Ni is the number of accepted jets in Et bin i of width AE^; Ci is the resolution 
unsmearing correction; Li is the integrated luminosity; e, is the efficiency of the trigger, 
vertex selection, and the jet quality cuts; and Ar\ is the width of the pseudorapidity bin. 



A. Data Selection 



The selected data are events with one or more jets which satisfy the requirements of the 
inclusive jet triggers. Jets are required to pass the standard jet quality criteria to be included 
in the cross section sample (Section |VH1| ). The $ T of the event is required to satisfy Eq. |875| . 
The vertex of the event must be within 50 cm of z = 0. The efficiency for each jet is then 
given by the product of the efficiencies of the jet quality cuts (ej et ), the efficiency of the cut 
on $ T (e me t), the efficiency for an event to pass the trigger (et r i g ger), and the efficiency for 
passing the vertex cut (e ve rtex): 

Q ^jet^met^trigger^vertex- (H""^) 



The values of ej et and e me t are plotted in Fig. [L7]. The efficiency of the vertex requirement is 
90 ± 1%. 



B. Filter Efficiency and Luminosity Matching 

Figure (42] shows the cross section ratios for Jet_50/Jet_30, Jet_85/Jet_50 and 
Jet_115/Jet_85. Since the denominator in each ratio represents a less restrictive trigger 
than the numerator, the numerator trigger is efficient where the ratio stabilizes at a con- 
stant value. Thus Jet_50, Jet_85, and Jet_115 are efficient above 90, 130, and 170 GeV, 
respectively. The efficiency for Jet_30 was determined to be 100% at 50 GeV (Section [V]). 

The determination of the integrated luminosity for each of the jet triggers is described 
in detail in Section [VT|. The luminosity used for Jet_50 is determined by matching the 
Jet_50 inclusive jet cross section to the Jet_85 cross section above 130 GeV, introducing a 
1.1% statistical error. The Jet_30 luminosity is determined by matching to the Jet_50 cross 
section above 90 GeV, which results in a 1.4% statistical error. Hence the matching error 
for Jet_30 is given by 1.1% and 1.4% added in quadrature, or 1.7%. These errors are added 
to the 5.8% error on Jet_85. The final Jet_30 and Jet_50 luminosities are then 0.350 pb _1 
and 4.76 pb -1 with errors of 6.1% and 5.9%, respectively. 

Figure ^ shows E^ spectra for the four jet triggers, without luminosity normalization, 
in the central rapidity region < 0.5) after efficiency and energy corrections. 
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FIG. 42. Inclusive cross section ratios. The arrows signify the Et above which the higher 
threshold trigger is used. 
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FIG. 43. Energy-corrected Et spectra for Jet_115 (solid line), Jet_85 (dashed), Jet_50 (dotted), 
and Jet_30 (dot-dashed). The arrows signify the Et range in which each trigger's spectrum is used. 
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C. Observed Cross Section 



Figure |44| shows the central cross section compiled from the four triggers. As suggested by 
the cross section ratios, and in order to maximize statistics, the spectrum from 60 < Ey < 90 
GeV is taken from the Jet_30 data, 90-130 GeV from Jet_50, 130-170 GeV from Jet_85, and 
above 170 GeV from Jet_115. The three data sets in Fig. [H] correspond to the low, nominal, 
and high energy scale corrections. The differences can be considered to be an error estimate 
on the cross section which dominates all other sources of error (luminosity, jet, missing E^, 
and vertex cuts). 
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FIG. 44. Energy-corrected and luminosity-normalized E^ spectra. The points with error bars 
correspond to the nominal energy scale correction. The dashed (dotted) histogram corresponds to 
the high (low) energy scales corrections. 



D. Highest Et Event Scanning 



Since the cross section decreases rapidly as the E^ increases, a small amount of contam- 
ination can have a significant effect on the measured cross sections at large E^. The data 
set included 46 events that passed selection cuts and contained a central jet (| r] \ < 0.7) with 
transverse energy greater than 375 GeV. These events were visually scanned for contamina- 
tion. We defined an event to be "good" if it had at least two jets with well-contained energy, 
if there were no isolated cells forming jets, and if there was no activity in the muon chambers 
consistent with cosmic ray interactions associated with the event. These conditions were in- 
tended to reject high-i?T jets arising from noisy calorimeter cells, cosmic rays, or beam halo 
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from the main ring, which passes through the D0 detector [pqj . The 46 events contained 62 
jets with £"x greater than 375 GeV. Seven of these jets included restored cells and seven of 
the events preferred the second vertex. All of the events passed visual inspection. 

E. Resolution Unfolding 

The steep E^ spectrum is distorted by jet energy resolution. The distortion was corrected 
by using an ansatz function for the cross section, 

exp(A)E^ (l-^j^ , (11.3) 

smearing it with the measured resolution (Table |V|), and comparing the smeared result with 
the measured cross section. The parameters A, a, and j3 were varied until the best fit was 
found between the observed cross section and the smeared trial spectrum. The x 2 for the fit 
is 21.2 for 24 bins and three parameters, corresponding to 21 degrees of freedom (Table |TX|) . 



Figure 45 shows an example of the energy scale corrected data with the best-fit smeared 



and unsmeared ansatz functions. Simulations have shown that 77-smearing causes negligible 



changes in the inclusive cross sections [2G 



TABLE IX. Unsmearing ansatz function parameters for the inclusive jet cross section (in fb) 
at ,r S = 1800 GeV. 



Rapidity Range 


Parameter 


Value 




A 


37.28 


|t7 < 0.5 


a 


-5.04 






8.23 




A 


37.30 


0.1 < |t/| < 0.7 


a 


-5.05 




P 


8.37 



Figure |46| shows the unsmearing correction as a function of transverse energy. The 
observed cross section is multiplied by this correction. The central curve shows the nominal 
correction. The change in cross section is greatest at low E T due to the steepness of the 
inclusive spectra and the relatively poor, rapidly changing jet resolution. The magnitude of 
the correction is —13% at 64.6 GeV, drops to ~ —6% at 205 GeV, and then rises to —12% 
at 461 GeV. 

The two outer curves of Fig. [P] show the extent of the uncertainties in the nominal 
correction due to the resolution uncertainties. This error was estimated directly with the 
data by unfolding with the upper and lower estimates of the resolution curves. For \ < 0.5 
the maximum error is 3%. Varying the fit parameters by up to 3 standard deviations results 
in negligible changes in the resolution correction. 

The resolution correction errors due to the fit procedure and statistical fluctuations of 
the data were estimated by performing the unfolding procedure on distributions simulated 
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FIG. 45. Data with smeared and unsmeared fit hypotheses (Eq. |11,3 ). The lower pane shows 
the smeared fit residuals, (data — smeared fit)/smeared fit. 
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with jetrad. A generated theoretical distribution was smeared with a resolution function. 
The ratio of the generated theoretical distribution to the smeared theoretical distribution 
was taken as the "true" unsmearing correction. Next, the previously described unfolding 
procedure was applied to the "smeared theory" and the resulting unsmearing correction was 
compared with the "true" unsmearing correction. The difference between the two corrections 
provided a measurement of the unfolding error. Above Et = 50 GeV, the differences were 
less than 1%. The error due to statistical fluctuation was estimated by simulating many 
jet samples containing the same total number of jets as the data sample. The statistical 
fluctuations between the different simulated samples lead to an error below 0.25% in any Et 
bin. A detailed description of this unfolding, and the unfolding error estimation procedures 
can be found in Ref. |H| . 



F. Unfolded Cross Section 

The central inclusive jet cross section is shown in Fig. H7I. The cross section values are 



plotted in each bin at the Et value for which the average integrated cross section is equal 
to the value of the analytical function (Eq. p.l.3| ) fitted to the cross-section The error 



bars are purely statistical and are visible only for the highest Et value. The error band 
indicates a one standard deviation variation of all systematic uncertainties, except the 5.8% 
uncertainty on the absolute normalization. The measured cross-section is compared to the 
inclusive cross section for the same Et values calculated with the jetrad program using 
the CTEQ3M PDF and the scale p = 0.5i?™ ax . This prediction lies within the error band 
for all Et bins. Table [X] lists the cross sections for 1 77 1 < 0.5 and 0.1 < 1 77 1 < 0.7. 

G. Cross Section Uncertainties 

The cross section uncertainties are dominated by the uncertainties in the energy scale 



correction. Table XI summarizes the uncertainties in the unfolded cross section. A detailed 



list of the uncertainties and their magnitudes is given in Tables [XI1| and |XIII| . Figure f|8 



shows the various uncertainties for the 1 77 1 < 0.5 cross section. The second uppermost curve 
shows the uncertainty in the energy scale, which varies from 8% at low Et to 30% at 450 
GeV. Clearly, this contribution dominates all other sources of error except at low Et where 
the 5.8% luminosity error is of comparable magnitude. The other sources of error (jet and 
event selection, trigger matching, and jet resolution) are relatively small. 

Most of the systematic uncertainties in the inclusive jet cross section are highly correlated 
as a function of Et- The uncertainties are separated into three "types," depending on the 
correlation (p) between two bins: 

• p = 1 : "Completely correlated," indicating that a la fluctuation in an error at a 
particular Et bin is accompanied by a la fluctuation at all other Et bins (Fig. |49|). 

• p = p(Eti,Et2) = [— 1,1} '■ "Partially correlated," possessing a varying degree of 
correlation in Et- A la fluctuation thus implies a less than la fluctuation elsewhere 
(Fig. |50D; negative p indicates the shifts will have opposite directions at the two points. 
This type of error is the most complicated to calculate and propagate. 
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TABLE X. The \ rj\ < 0.5 and 0.1 < | rj \ < 0.7 cross sections (Eq. 11.1). Also given is the value 
of the fit to the cross section using Eq. 11.3| . 



\y\ < 0.5 



0.1 < |r?| < 0.7 



Bin Range 


Plotted Et 


Cross Sec. 




Svs 


Fitted 


Plotted _Et 


Cross Sec. 




Svs 


Fitted 


(GeV) 


(GeV) 


± Stat 






Uncer.(%) 


Cross Sec. 


(GeV) 


± Stat 






Uncer.(%) 


Cross Sec. 








(fb/GeV/Ar?) 




(fb/GeV/A?j) 




(fb/GeV/Ar?) 




(fb/GeV/Ar;) 


60 - 


70 


64.6 


(6.39 ± 


04) 




10° 


±10 


6.27 


y 


10 B 


64.6 


( 6 26 + OA) 




10° 


±10 


6.13 


x 10 tt 


70 - 


80 


74.6 


(2.80 ± 


03) 




10 6 


±10 


2.74 




10° 


74.6 






10 6 


±10 


2.67 


x 10 6 


80 - 


90 


84.7 


( 1.36 ± 


02) 




10° 


±10 


1.31 




10 6 


84.7 






10 6 


±10 


1 28 


x 10 6 


90 - 


100 


94.7 


(6.84 ± 


04) 




10 5 


±10 


6.74 




10 5 


94.7 


iu.uu _i_ \j.\j'-t) 




10 5 


±10 


6.53 


x 10 5 


100 - 


110 


104.7 


(3.76 ± 


03) 


X 


10 5 


±10 


3.67 


X 


10 5 


104.7 


fa K3 + n n3~i 




10 s 


±10 


3.54 


x 10 5 


110 - 


120 


114.8 


(2.14 ± 


02) 


X 


10 5 


±10 


2.08 


X 


10 5 


114.8 


(2.07 ± 0.02) 


X 


10 s 


±10 


2.01 


x 10 5 


120 - 


130 


124.8 


(1.23 ±0 


02) 


X 


10"' 


±10 


1.23 


X 


10 5 


124.8 


(1.19 ± 0.01) 


X 


10 s 


±10 


1.18 


x 10 s 


130 - 


110 


134.8 


(7.46 ±0 


04) 


X 


10 4 


±10 


7.49 


X 


10 4 


134.8 


(7.16 ± 0.03) 


X 


10 4 


±10 


7.18 


x 10 4 


140 - 


150 


144.8 


(4.71 ± 


03) 


X 


10 4 


±10 


4.69 


X 


10 4 


144.8 


(4.51 ± 0.03) 


X 


If) 4 


±10 


4.48 


x 10 4 


150 - 


160 


154.8 


(2.97 ± 


02) 


X 


10 4 


±10 


3.00 


X 


10 4 


154.8 


(2.83 ± 0.02) 


X 


If) 4 


+11,-10 


2.86 


x 10 4 


160 - 


170 


164.8 


(1.94 ±0 


02) 


X 


10 4 


+11,-10 


1.96 


X 


10 4 


164.8 


(1.83 ±0.02) 


X 


10 4 


±11 


1.86 


x 10 4 


170 - 


180 


174.8 


(1.30 ±0 


01) 


X 


10 4 


±11 


1.30 


X 


10 4 


174.8 


(1.23 ± 0.01) 


X 


10 4 


±11 


1.23 


x 10 4 


180 - 


190 


184.8 


(8.83 ±0 


10) 


X 


10 3 


±11 


8.75 


X 


10 3 


184.8 


(8.38 ± 0.09) 


X 


10 3 


±11 


8.28 


x 10 3 


190 - 


200 


194.8 


(5.95 ±0 


08) 


X 


10 3 


±11 


5.98 


X 


10 3 


194.8 


(5.64 ± 0.07) 


X 


10 3 


±12,-11 


5.64 


x 10 3 


200 - 


210 


204.8 


(4.15 ±0 


07) 


X 


10 3 


±12,-11 


4.13 


X 


10 3 


204.8 


(3.93 ± 0.06) 


X 


10 3 


±12,-11 


3.88 


x 10 3 


210 - 


220 


214.8 


(2.84 ± 


06) 


X 


10 3 


+12,-11 


2.88 


X 


10 :i 


214.8 


(2.67 ±0.05) 


X 


10 3 


±12 


2.70 


x 10 3 


220 - 


230 


224.8 


(2.08 ±0 


05) 


X 


10 3 


±12 


2.03 


X 


10 3 


224.8 


(1.95 ±0.04) 


X 


10 3 


±13,-12 


1.90 


x 10 3 


230 - 


250 


239.4 


(1.26 ±0 


03) 


X 


10 3 


+13,-12 


1.24 


X 


10 3 


239.4 


(1.17 ±0.02) 


X 


10 3 


±13,-12 


1.15 


x 10 3 


250 - 


270 


259.4 


(6.34 ± 


19) 


X 


10 2 


+14,-13 


6.40 


X 


10 2 


259.4 


(5.84 ±0.17) 


X 


10 2 


±14,-13 


5.94 


x 10 2 


270 - 


290 


279.5 


(3.65 ±0 


15) 


X 


10 2 


+14,-13 


3.39 


X 


10 2 


279.5 


(3.21 ± 0.12) 


X 


10 2 


±15,-14 


3.13 


x 10 2 


290 - 


320 


303.9 


(1.73 ±0 


08) 


X 


10 2 


+16,-14 


1.60 


X 


10 2 


303.9 


(1.56 ±0.07) 


X 


10-' 


±16,-15 


1.47 


x 10 2 


320 - 


350 


333.9 


(6.60 ±0 


50) 


X 


10 1 


+17,-16 


6.50 


X 


10 1 


333.9 


(6.06 ± 0.44) 


X 


10 1 


±18,-16 


5.91 


x 10 1 


350 - 


410 


375.7 


(1.83 ±0 


19) 


X 


10 1 


+21,-18 


1.91 


X 


10 1 


375.7 


(1.48 ±0.15) 


X 


10 1 


±22,-19 


1.72 


x 10 1 


410 - 


560 


461.1 


(1.20 ±0 


30) 


X 


10° 


+30,-25 


1.57 


X 


10° 


460.9 


(1.05 ±0.25) 


X 


10° 


±31,-26 


1.39 


x 10° 



TABLE XI. Unfolded cross section errors 



Source 


Percentage 


Comment 


Jet and Event Selection 


< 2 


Correlated 


Luminosity 


5.8 


Correlated 


Luminosity Match 






60- 90 GeV 


1.7 


Statistical, Correlated 


90-130 GeV 


1.1 


Trigger-to- Trigger 


Energy Scale 


15-30 


Mostly Correlated 


Unfolding 






Resolution Function 


1-3 


Correlated 


Closure 


1-2 


Correlated 
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FIG. 47. The \rj\ < 0.5 inclusive jet cross section. Statistical uncertainties are invisible on 
this scale except for the highest bin. The histogram represents the jetrad prediction and the 
shaded band represents the ±lcx systematic uncertainty band about the prediction excluding the 
5.8% luminosity uncertainty. 
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• p = : "Uncorrelated," statistical in nature or otherwise independent of one another. 
Some small errors with unknown (but probably positive) Et correlation are treated as 
uncorrelated for simplicity. Such treatment is conservative. 

The uncertainties due to jet selection are correlated as a function of E^. The uncertainties 
due to unsmearing are also correlated. The luminosity uncertainty is correlated as a function 
of i?T- The trigger matching uncertainties are correlated as a function of E^ for bins that 
are derived from the same trigger sample and uncorrelated for all other bins. The energy 
scale errors are partially correlated as a function of Ey and are discussed below. 

The energy scale calibration is implemented as a series of corrections, each with its own 
uncertainty (Section |TX|). The uncertainty due to the energy scale is separated into several 
components so that the correlations as a function of E T can be studied (Fig. ^T|). The 
energy scale uncertainties were calculated with a Monte Carlo simulation of the inclusive 
jet cross section. At each uncorrected Ei the simulation generated an ensemble of jets 
with rapidity, vertex position, luminosity, and variable correlations derived from the data. 
Figure [51] shows the components of energy scale uncertainty as a function of E^. The E^ of 
each of the simulated jets was then corrected and the resulting uncertainty in the jet cross 
section calculated. These uncertainties are in good agreement with the uncertainties derived 
from the data. 

The uncertainties due to the offset correction, the 77-dependent correction, the showering 
correction, and the method are all correlated as a function of E^. The hadronic response 
uncertainty is partially correlated as a function of E T (Section |TX| ). The hadronic response 
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FIG. 49. Example of an error band relative to some nominal distribution (illustrated here with 
a flat line). If the errors at points 1 and 2 are completely correlated, then a one standard deviation 
(la ) Ai at the first position necessarily results in a la A2 at the second position. 



A 


A 
A 2 , 

( Nominal ^ 




2 



FIG. 50. If the errors at points 1 and 2 are partially correlated, then a full la Ai at the first 
position results in a smaller than la A 2 ' at the second position. The correlation factor illustrated 
here is 0.55. 
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FIG. 51. Percentage cross section errors for 
scale correction. 
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1 77 1 < 0.5 associated with the components of energy 



correlations are illustrated in Fig. |52| and are given in Table |XIV| . In addition, the response 
uncertainties are only approximated by a Gaussian uncertainty distribution. Tables |XV| and 
XVI| give the actual uncertainties for a given percentage confidence level (C.L.). i.e. if one 



has a +20% error in the cross section at a given corresponding to 95% C.L., then with 
95% probability the response errors will cause a deviation in the cross section of < 20%. 
The correlations for the total systematic uncertainties are given in Table |XV1J . 



H. Comparison of the Data to Theory 

Figures |53] and |54| show the fractional difference between the data, D, and a jetrad 
theoretical prediction, T, normalized by the prediction, ((D — T)/T), for \rj\ < 0.5 and 
0.1 < 1 77 1 < 0.7 respectively. The jetrad prediction was generated with /1 = 0.5i?™ ax , 
TZ sep = 1.3, and several different choices of PDF. The error bars represent statistical errors 
only. The outer bands represent the total cross section error excluding the 5.8% luminosity 
uncertainty. Given the experimental and theoretical uncertainties, the predictions are in 
agreement with the data; in particular, the data above = 350 GeV show no indication 
of an excess relative to QCD. 

The data and theory can be compared quantitatively with a x 2 test incorporating the 
uncertainty covariance matrix (Table |XVII| , ||55| ). The \ 2 is given by: 



x 2 = E^7%- 



(11.4) 
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FIG. 52. The correlations of the uncertainty due to the hadronic response correction as a 
function of Et. The solid curve shows the correlations relative to the 461 GeV bin, the dashed 
curve with respect to the 205 GeV bin, and the dotted curve with respect to the 105 GeV bin. 



TABLE XII. Percentage \ rj\ < 0.5 cross section uncertainties. The last row gives the nature of 
the Et bin-to-bin correlations: signifies uncorrelated uncertainties, 1 correlated, and p partially- 



correlated. 



£7t 


Stat 


Jet 


Lumin 


Lumin Unsmcaring 










Energy Scale 








GeV 


Error 


Sel 




Match 


High 


Low 


Underlying 


V 




Method 


Shower 


Response 
















High 


Low 


High 


Low 


High 


Low 


High 


Low 


High 


Low 


64.6 


0.7 


0.5 


5.8 


1.8 


2.5 


-2.6 


5.3 


-5.0 


0.2 


-0.2 


3.4 


-3.2 


4.7 


-4.6 


1.1 


-1.0 


74.6 


1.0 


0.5 


5.8 


1.8 


2.3 


-2.4 


4.7 


-4.5 


0.2 


-0.2 


3.5 


-3.3 


4.8 


-4.7 


1.2 


-1.1 


84.7 


1.5 


0.5 


5.8 


1.8 


2.1 


-2.2 


4.2 


-4.1 


0.2 


-0.2 


3.5 


-3.5 


4.9 


-4.8 


1.5 


-1.3 


94.7 


0.6 


0.5 


5.8 


1.1 


1.9 


-2.1 


3.8 


-3.7 


0.2 


-0.2 


3.6 


-3.5 


5.0 


-4.9 


1.9 


-1.9 


104.7 


0.8 


0.5 


5.8 


1.1 


1.9 


-2.1 


3.5 


-3.4 


0.2 


-0.2 


3.6 


-3.6 


5.1 


-5.1 


2.2 


-2.2 


114.8 


1.0 


0.5 


5.8 


1.1 


1.8 


-2.0 


3.2 


-3.1 


0.2 


-0.2 


3.7 


-3.6 


5.3 


-5.1 


2.6 


-2.5 


124.8 


1.4 


0.5 


5.8 


1.1 


1.8 


-2.0 


3.0 


-2.9 


0.2 


-0.2 


3.9 


-3.7 


5.4 


-5.2 


2.8 


-2.8 


134.8 


0.5 


0.5 


5.8 


0.0 


1.7 


-2.0 


2.9 


-2.8 


0.2 


-0.2 


4.1 


-3.9 


5.6 


-5.4 


3.1 


-3.0 


144.8 


0.6 


0.5 


5.8 


0.0 


1.7 


-2.0 


2.7 


-2.6 


0.2 


-0.2 


4.1 


-3.9 


5.6 


-5.5 


3.3 


-3.2 


154.8 


0.8 


0.5 


5.8 


0.0 


1.7 


-2.0 


2.5 


-2.5 


0.3 


-0.3 


4.3 


-3.9 


5.9 


-5.6 


3.7 


-3.4 


164.8 


1.0 


0.5 


5.8 


0.0 


1.7 


-2.1 


2.4 


-2.4 


0.3 


-0.3 


4.5 


-4.1 


5.9 


-5.7 


3.9 


-3.6 


174.8 


0.9 


0.5 


5.8 


0.0 


1.7 


-2.1 


2.3 


-2.3 


0.3 


-0.3 


4.7 


-4.2 


6.1 


-5.9 


4.1 


-3.9 


184.8 


1.1 


0.5 


5.8 


0.0 


1.8 


-2.1 


2.2 


-2.2 


0.3 


-0.3 


4.9 


-4.3 


6.2 


-6.0 


4.3 


-4.2 


194.8 


1.4 


0.5 


5.8 


0.0 


1.8 


-2.2 


2.1 


-2.1 


0.3 


-0.3 


5.1 


-4.4 


6.4 


-6.2 


4.5 


-4.4 


204.8 


1.7 


0.5 


5.8 


0.0 


1.8 


-2.2 


2.1 


-2.0 


0.3 


-0.3 


5.3 


-4.5 


6.5 


-6.3 


4.8 


-4.6 


214.8 


2.0 


0.5 


5.8 


0.0 


1.9 


-2.3 


2.0 


-2.0 


0.3 


-0.3 


5.6 


-4.5 


6.7 


-6.4 


5.2 


-5.0 


224.8 


2.4 


0.5 


5.8 


0.0 


1.9 


-2.4 


2.0 


-1.9 


0.3 


-0.3 


5.8 


-4.7 


6.9 


-6.6 


5.4 


-5.2 


239.4 


2.1 


0.5 


5.8 


0.0 


2.0 


-2.5 


1.9 


-1.9 


0.3 


-0.3 


6.1 


-4.9 


7.1 


-6.8 


5.8 


-5.7 


259.4 


3.0 


0.5 


5.8 


0.0 


2.1 


-2.6 


1.8 


-1.8 


0.3 


-0.3 


6.6 


-5.1 


7.5 


-7.2 


6.4 


-6.3 


279.5 


4.0 


0.6 


5.8 


0.0 


2.2 


-2.8 


1.8 


-1.8 


0.4 


-0.3 


7.1 


-5.4 


7.9 


-7.5 


7.1 


-6.9 


303.9 


4.7 


0.6 


5.8 


0.0 


2.4 


-3.1 


1.8 


-1.8 


0.4 


-0.4 


7.9 


-5.7 


8.4 


-7.9 


8.3 


-7.8 


333.9 


7.6 


0.7 


5.8 


0.0 


2.7 


-3.5 


1.8 


-1.7 


0.4 


-0.4 


9.0 


-6.1 


9.1 


-8.6 


9.8 


-9.3 


375.7 


10.2 


1.0 


5.8 


0.0 


3.2 


-4.2 


1.7 


-1.7 


0.5 


-0.5 


10.8 


-6.8 


10.2 


-9.6 


12.4 


-11.7 


461.1 


25.0 


2.1 


5.8 


0.0 


4.6 


-5.9 


1.7 


-1.7 


0.6 


-0.6 


15.0 


-8.6 


13.2 


-12.0 


20.3 


-18.2 


Corrcl. 





1 


1 


p 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


P 


P 



78 



0.25 


-0.25 

0.25 


-0.25 

5 

M 



o 
o 

I 

03 

Q 



0.25 


-0.25 

0.25 


-0.25 

0.25 


-0.25 

0.25 


-0.25 



CTEQ3M 



MjJ<0-5 
I 



CTEQ4M 



CTEQ4HJ 



MRS(A') 



MRST 








MRST(gi) 
I 



100 



200 



300 



400 



E T (GeV) 



FIG. 53. The difference between data and jetrad QCD predictions normalized to predictions 
for 1 77 1 < 0.5. The shaded region represents the ±la systematic uncertainties about the prediction. 
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FIG. 54. The difference between data and jetrad QCD predictions normalized to predictions 
for 0.1 < 1 17 1 < 0.7. The error bars show the statistical uncertainties. The shaded region represents 
the ±1(7 systematic uncertainties about the prediction. 
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TABLE XIII. Percentage 0.1 < \rj\ < 0.7 cross section uncertainties. The last row gives the 
nature of the bin-to-bin correlations: signifies uncorrelated uncertainties, 1 correlated, and p 



partially correlated. 





Stat 


Jot 


Lumin 


Lumin Unsmearing 










Enerj 


$y Scale 








GeV 


Error 


Sel 




Match 


High 


Low 


Underlying 


V 




Method 


Shower 


Response 
















High 


Low 


High 


Low 


High 


Low 


High 


Low 


High 


Low 


64.6 


0.6 


0.5 


5.8 


1 s 
l .0 


9 °, 


9 fi 
-z.u 


o.o 


-o.u 


0.6 


-0.6 


A 

0.1 


9 

-o.z 


A 8 
-l.o 


A 8 


i i 
i . i 


i 1 

-i . i 


74.6 


0.9 


0.5 


5.8 


1 S 


9 1 
Z . 1 


9 A 


A 7 
'1. t 


A ^ 


0.6 


-0.6 


. > . ■ ) 


^ A 


1 
i. y 


A 7 


l.O 


i 1 


84.7 


1.4 


0.5 


5.8 


1 Si 
1 .0 


9 n 

z . u 


9 'X 
-z.o 


1 9 
1. Z 


A 1 


0.6 


-0.6 


O . 


-O.O 


r l ft 
o.u 


A Q 


"1 fi 


-l.O 


94.7 


0.5 


0.5 


5.8 


2 i 




-2 2 


' \ si 
o.o 


-o. / 


0.6 


-0.6 


'3 fi 


-o.o 


0. 1 


-o.u 


9 "1 
Z. 1 


9 1 
-Z. 1 


104.7 


0.7 


0.5 


5.8 


i i 


1 .0 


9 9 
-Z. Z 


O.O 




0.7 


-0.7 


O. ( 


^ fi 
-o.o 


o.z 


9 

-o.z 


9 '3. 
Z.O 


9 1 
-Z.O 


114.7 


1.0 


0.5 


5.8 


j 2 


1 .0 


-2 1 


,; t 9 

o. z 


-o. 1 


0.7 


-0.7 


O . o 


^ fi 

-O.D 


- r. 
0.0 


^ 9 

-o.z 


9 7 
Z. / 


9 fi 
-Z.U 


124.8 


1.3 


0.5 


5.8 


i i 
i . i 


1 .0 


9 1 
-Z. 1 


o. U 


9 O 

-z.y 


0.7 


-0.7 


o.y 


Q T 
-O. f 


0.0 


-0.4 


9 O 

Z.y 


9 O 

-Z.y 


134.8 


0.5 


0.5 


5.8 


u.u 


1 .0 


9 1 
-Z. L 


9 ^ 
Z.o 


9 £ 
-Z.o 


0.7 


-0.7 


-i 1 


-o.y 


^ 7 
0. i 


-0.0 


9 

o. z 


-0. ± 


144.8 


0.6 


0.5 


5.8 


o n 

u.u 


1 .0 


-2 2 


9 7 
z. / 


9 fi 
-z.u 


0.7 


-0.7 


4 2 


-■41. U 


O.O 


^ fi 

-o.u 


fi 
o.o 


-o.o 


154.8 


0.8 


0.5 


5.8 


o n 

u.u 


1 .0 


-2 2 


2.5 


-2.5 


0.8 


-0.8 


4.4 


-4.0 


5.9 


-5.8 


3.7 


-3.5 


164.8 


0.9 


0.5 


5.8 


0.0 


1.8 


-2.2 


2.4 


-2.4 


0.8 


-0.8 


4.6 


-4.2 


6.0 


-5.9 


4.0 


-3.8 


174.8 


0.9 


0.5 


5.8 


0.0 


1.8 


-2.3 


2.3 


-2.3 


0.8 


-0.8 


4.8 


-4.2 


6.2 


-6.0 


4.2 


-4.0 


184.8 


1.1 


0.5 


5.8 


0.0 


1.9 


-2.4 


2.2 


-2.2 


0.8 


-0.8 


5.0 


-4.3 


6.4 


-6.2 


4.4 


-4.4 


194.8 


1.3 


0.5 


5.8 


0.0 


1.9 


-2.4 


2.1 


-2.1 


0.8 


-0.8 


5.2 


-4.4 


6.5 


-6.3 


4.7 


-4.4 


204.8 


1.6 


0.5 


5.8 


0.0 


1.9 


-2.5 


2.1 


-2.0 


0.9 


-0.9 


5.4 


-4.5 


6.7 


-6.5 


5.0 


-4.8 


214.8 


1.9 


0.5 


5.8 


0.0 


2.0 


-2.6 


2.0 


-2.0 


0.9 


-0.9 


5.7 


-4.6 


6.9 


-6.6 


5.3 


-5.0 


224.8 


2.2 


0.5 


5.8 


0.0 


2.0 


-2.7 


2.0 


-2.0 


0.9 


-0.9 


6.0 


-4.8 


7.0 


-6.7 


5.6 


-5.5 


239.4 


2.0 


0.5 


5.8 


0.0 


2.1 


-2.8 


1.9 


-1.9 


1.0 


-0.9 


6.2 


-4.9 


7.3 


-6.9 


6.0 


-5.8 


259.4 


2.9 


0.5 


5.8 


0.0 


2.3 


-3.0 


1.9 


-1.9 


1.0 


-1.0 


6.8 


-5.1 


7.7 


-7.2 


6.7 


-6.4 


279.5 


3.9 


0.6 


5.8 


0.0 


2.4 


-3.3 


1.8 


-1.8 


1.1 


-1.0 


7.4 


-5.4 


8.0 


-7.6 


7.5 


-7.2 


303.9 


1.5 


0.6 


5.8 


0.0 


2.7 


-3.6 


1.8 


-1.8 


1.1 


-1.1 


8.1 


-5.7 


8.5 


-8.0 


8.7 


-8.3 


333.9 


7.2 


0.7 


5.8 


0.0 


3.0 


-4.1 


1.8 


-1.8 


1.2 


-1.2 


9.3 


-6.2 


9.2 


-8.6 


10.2 


-9.8 


375.7 


10.3 


1.0 


5.8 


0.0 


3.6 


-4.9 


1.8 


-1.7 


1.3 


-1.3 


11.1 


-6.9 


10.3 


-9.7 


13.5 


-12.5 


460.9 


24.3 


2.1 


5.8 


0.0 


5.2 


-7.0 


1.7 


-1.7 


1.6 


-1.6 


15.1 


-8.7 


13.3 


-12.2 


21.8 


-19.4 


Correl. 





1 


1 


p 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


P 


P 



where <5j is the difference between the data and theory for a given bin, and is element 
i, j of the covariance matrix: 

Vij = Pij ■ AcTj • Act,-. (If. 5) 

where Act is the sum of the systematic error and the statistical error added in quadrature 
if % — j and the systematic error if i ^ j, and is the correlation between the systematic 
uncertainties of bins as given in Table [XVI]] . The systematic uncertainty is given by 



the percentage uncertainty times the theoretical prediction (see Appendix for a discussion 
of the x 2 )- The resulting x 2 values are given in Table |XVHI| for all of the theoretical choices 
described above. The choice of PDF and renormalization scale is varied. Each comparison 
has 24 degrees of freedom. 

All but one of the jetrad predictions adequately describe the \r]\ < 0.5 and 0.1 < 
< 0.7 cross sections. For these, the probabilities for x 2 t° exceed the listed values are 
between 10% and 86%. The prediction using CTEQ4HJ and /i = 0.5-E™ ax produces the 
highest probability for both measurements. The prediction with the MRST(g j) PDF has a 
probability of agreement with the data of 0.3%, and thus is incompatible with our data. 

Comparisons between the data and eks calculations using various PDFs, 7Z scp = 1.3, 
and with renormalization scales fi = (0.25, 0.50, 0.75, 1.00, 1.25, 1.50, 1.75, 2.00)i?T (where 
E T = £™ ax and E^ 1 ) are also made (Table P0X| ). The eks predictions give a reasonable 
description of the \ \ < 0.5 cross section. However, unlike the jetrad predictions, the eks 
calculation using CTEQ4M and fi = 0.5E^ ax has the highest probability of agreement. The 
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TABLE XIV. The correlations for the uncertainty due to the energy scale response for 
{rf^l < 0.5, and 0.1 < \rf ,ct \ < 0.7. The correlation values above the diagonal are the corre- 
lations for | if ct | < 0.5 and the correlations below the diagonal correspond to 0.1 < | rf ct \ < 0.7. In 



both cases the correlation matrices are symmetric. 
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TABLE XV. The percentage cross section uncertainties due to the energy scale response cor- 
rection that correspond to a given percentage confidence level for 1 77 1 < 0.5. 



Upper Lower 



Bin 


40% 


68.3% 


86% 


90% 


95% 


99% 


40% 


68.3% 


86% 


90% 


95% 


99% 


1 


0.7 


1.1 


1.3 


1.3 


1.6 


1.9 


-0.7 


-1.0 


-1.3 


-1.3 


-1.6 


-1.8 


2 


0.7 


1.1 


1.5 


1.4 


1.8 


2.2 


-0.7 


-1.1 


-1.5 


-1.4 


-1.7 


-2.1 


3 


1.0 


1.5 


1.9 


1.9 


2.2 


2.6 


-0.9 


-1.3 


-1.8 


-1.8 


-2.2 


-2.6 


4 


1.3 


1.9 


2.5 


2.6 


3.0 


3.6 


-1.3 


-1.8 


-2.4 


-2.6 


-3.0 


-3.5 


5 


1.4 


2.2 


3.0 


3.1 


3.4 


4.0 


-1.4 


-2.2 


-2.9 


-3.0 


-3.3 


-4.0 


6 


1.6 


2.6 


3.3 


3.4 


3.8 


4.6 


-1.5 


-2.5 


-3.1 


-3.3 


-3.6 


-4.4 


7 


1.6 


2.9 


3.6 


3.8 


4.1 


5.0 


-1.6 


-2.8 


-3.5 


-3.7 


-4.0 


-4.8 


8 


1.6 


3.0 


3.9 


4.0 


4.2 


5.3 


-1.7 


-3.0 


-3.8 


-3.0 


-4.2 


-5.3 


9 


1.7 


3.3 


4.1 


4.2 


4.6 


5.7 


-1.7 


-3.2 


-4.1 


-4.2 


-4.5 


-5.6 


10 


1.9 


3.6 


4.5 


4.6 


5.0 


6.2 


-1.8 


-3.4 


-4.3 


-4.4 


-4.7 


-5.9 


11 


2.0 


3.9 


4.8 


4.9 


5.1 


6.7 


-1.8 


-3.6 


-4.6 


-4.6 


-4.9 


-6.3 


12 


2.1 


4.1 


5.1 


5.0 


5.4 


7.0 


-2.0 


-3.9 


-5.0 


-4.9 


-5.2 


-6.8 


13 


2.1 


4.3 


5.3 


5.2 


5.7 


7.4 


-2.1 


-4.1 


-5.2 


-5.2 


-5.6 


-7.2 


14 


2.4 


4.5 


5.6 


5.6 


6.1 


7.8 


-2.3 


-4.5 


-5.5 


-5.5 


-6.0 


-7.5 


15 


2.6 


4.8 


5.9 


5.9 


6.5 


8.3 


-2.5 


-4.6 


-5.7 


-5.8 


-6.4 


-7.9 


16 


3.1 


5.2 


6.4 


6.5 


7.3 


9.0 


-2.9 


-5.0 


-6.2 


-6.3 


-7.0 


-8.6 


17 


3.3 


5.4 


6.7 


6.9 


7.7 


9.4 


-3.1 


-5.2 


-6.4 


-6.6 


-7.5 


-8.9 


18 


3.8 


5.8 


7.1 


7.5 


8.8 


10.2 


-3.7 


-5.6 


-7.0 


-7.4 


-8.4 


-9.8 


19 


4.4 


6.4 


8.2 


8.5 


10.0 


11.7 


-4.3 


-6.3 


-7.8 


-8.1 


-9.4 


-11.0 


20 


5.0 


7.1 


9.0 


9.4 


11.1 


13.0 


-4.9 


-6.9 


-8.8 


-9.0 


-10.6 


-12.3 


21 


6.0 


8.3 


10.8 


11.1 


13.3 


15.6 


-5.7 


-7.8 


-10.0 


-10.2 


-12.0 


-14.0 


22 


7.4 


10.0 


13.3 


13.3 


15.8 


18.9 


-7.0 


-9.4 


-12.3 


-12.2 


-14.5 


-17.0 


23 


9.4 


12.6 


17.1 


16.8 


20.0 


24.4 


-9.0 


-12.0 


-15.7 


-15.6 


-18.2 


-21.4 


24 


16.2 


22.4 


29.6 


30.2 


35.3 


43.1 


-14.6 


-20.1 


-25.5 


-26.2 


-29.6 


-34.6 
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TABLE XVI. The percentage cross section uncertainties due to the energy scale response cor- 
rection that correspond to a given percentage confidence level for 0.1 < \ r]\ < 0.7. 



Upper Lower 
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40% 


68.3% 
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95% 
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1 O 

lz 
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4.Z 


O.O 


0. 1 


o.u 


7 9 
( .z 


-z.u 


-4.U 


-D. 1 


-o.u 


-0.4 


-u.y 


13 


2.3 


4.4 


5.5 


5.5 


6.0 


7.7 


-2.4 


-4.3 


-5.5 


-5.5 


-6.1 


-7.6 


14 


2.6 


4.7 


5.8 


5.9 


6.4 


8.1 


-2.4 


-4.5 


-5.6 


-5.6 


-6.2 


-7.7 


15 


2.9 


5.0 


6.2 


6.3 


7.2 


8.8 


-2.8 


-4.8 


-6.0 


-6.1 


-6.7 


-8.2 


16 


3.2 


5.3 


6.6 


6.8 


7.7 


9.2 


-3.0 


-5.1 


-6.3 


-6.5 


-7.4 


-8.8 


17 


3.6 


5.7 


7.0 


7.3 


8.3 


9.8 


-3.6 


-5.5 


-6.9 


-7.2 


-8.2 


-9.6 


18 


4.0 


6.0 


7.5 


7.9 


9.3 


10.9 


-4.0 


-5.9 


-7.3 


-7.6 


-8.8 


-10.3 


19 


4.7 


6.7 


8.6 


8.9 


10.4 


12.2 


-4.5 


-6.5 


-8.1 


-8.4 


-9.8 


-11.4 


20 


5.4 


7.6 


9.9 


10.1 


12.1 


14.2 


-5.2 


-7.3 


-9.4 


-9.6 


-11.3 


-13.2 


21 


6.5 


8.8 


11.8 


11.9 


14.2 


16.9 


-6.2 


-8.2 


-10.8 


-10.9 


-12.8 


-15.0 


22 


7.8 


10.5 


14.1 


13.9 


16.5 


19.8 


-7.4 


-9.7 


-12.8 


-12.6 


-14.9 


-17.8 


23 


9.9 


13.4 


18.2 


18.0 


21.4 


26.5 


-9.7 


-13.0 


-17.0 


-17.1 


-19.8 


-23.5 


24 


17.2 


24.4 


32.2 


33.0 


38.3 


47.3 


-15.9 


-22.1 


-28.0 


-28.8 


-32.2 


-37.1 



TABLE XVII. The systematic error correlations for the inclusive jet cross section for 
\if ct \ < 0.5, and 0.1 < \rf ot \ < 0.7. The correlation values above the diagonal are the corre- 
lations for | if ct | < 0.5 and the correlations below the diagonal correspond to 0.1 < | if ct \ < 0.7. In 



both cases the correlation matrices are symmetric. 





1 




2 




3 


4 


5 




6 


7 


8 




9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


21 


22 


23 


24 






1.00 


1 


00 





.99 


0.98 


0.97 





.96 


0.95 


0.93 





.92 





.91 


0.90 


0.89 


0.88 





.86 





.86 


0.84 


0.83 





.82 


0.80 


0. 


.78 





.75 


0.71 


0.66 


0.53 
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00 


1 
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0.99 
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0.96 


0.95 





.94 





.93 


0.92 


0.91 
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.89 





.88 


0.86 
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.84 


0.82 


0. 
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.77 


0.73 


0.68 


0.55 


2 


1 


1.00 






1 


.00 


0.99 


0.99 





.98 


0.97 


0.96 





.95 





.94 


0.93 


0.92 


0.92 
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.90 


0.88 


0.88 





.86 


0.84 





.82 





.79 


0.75 


0.70 


0.57 


3 


2 
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1 
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1.00 


1.00 





.99 


0.99 


0.98 
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.96 


0.96 


0.95 


0.94 
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.92 


0.91 


0.90 





.89 


0.86 


0. 


.84 





.81 


0.77 


0.72 


0.59 


4 


3 


0.99 


1 


.00 


1 


.00 




1.00 


1 


.00 


1.00 


0.99 





.98 





.97 


0.97 


0.96 


0.95 





94 





.94 


0.92 


0.92 





.90 


0.88 


0. 


.86 





.83 


0.79 


0.74 


0.61 


5 


4 


0.98 





.99 





.99 


1.00 




1 


.00 


1.00 


0.99 





.99 





.98 


0.98 


0.97 


0.96 


0. 


96 





.95 


0.94 


0.93 





.92 


0.90 





.88 





.85 


0.81 


0.76 


0.63 


6 


5 


0.97 





98 





.99 


1.00 


1.00 






1.00 


0.99 





.99 





.99 


0.98 


0.98 


0.97 


0. 


97 





.96 


0.95 


0.94 





.93 


0.91 


0. 


.89 





.87 


0.83 


0.78 


0.65 


7 


6 


0.96 


0. 


97 





.98 


1.00 


1.00 


1 


.00 




1.00 


1 


.00 


1 


.00 


0.99 


0.99 


0.99 


0. 


.98 





.97 


0.97 


0.96 





.95 


0.93 


0. 


91 





.89 


0.85 


0.80 


0.68 


8 


7 


0.95 





.96 





.98 


0.99 


1.00 


1 


.00 


1.00 




1 


.00 


1 


.00 


1.00 


0.99 


0.99 





99 





.98 


0.97 


0.97 





.96 


0.94 





.93 





.90 


0.87 


0.82 


0.70 


9 


8 


0.93 





.95 





.96 


0.98 


0.99 





.99 


0.99 


1.00 






1 


.00 


1.00 


1.00 


1.00 





99 





.99 


0.98 


0.98 





.97 


0.95 


0. 


.94 





.92 


0.88 


0.84 


0.72 


10 


9 


0.92 





,94 





.95 


0.97 


0.98 





.99 


0.99 


1.00 


1 


.00 






1.00 


1.00 


1.00 


1. 


.00 





.99 


0.99 


0.98 





.98 


0.96 


0. 


95 





.93 


0.90 


0.85 


0.74 


11 


10 


0.91 





93 





.95 


0.97 


0.98 





.98 


0.99 


1.00 


1 


.00 


1 


.00 




1.00 


1.00 


1 


.00 


1 


.00 


0.99 


0.99 





.98 


0.97 


0. 


96 





.94 


0.91 


0.87 


0.76 


12 


11 


0.90 





.92 





.94 


0.96 


0.97 





.98 


0.98 


0.99 


1 


.00 


1 


.00 


1.00 




1.00 


1 


.00 


1 


.00 


1.00 


0.99 





.99 
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TABLE XVIII. x 2 comparisons between jetrad and \r} \ < 0.5 and 0.1 < \ r) \ < 0.7 data for 
fi = O.bE™ 11 *, TZ scp = 1.3, and various PDFs. There are 24 degrees of freedom. 



PDF 


A* 




\rj\ < 0.5 


0.1 < 


| Tj | < 0.7 






, 2 

X 


Frob. 


, 2 

X 


rrob. 


CTEQ3M 


O^OE™ 3 * 


25.3 


0.39 


32.7 


0.11 


CTEQ4M 


0.50£™ ax 


20.1 


0.69 


26.8 


0.31 


CTEQ4HJ 


0.50£™ ax 


16.8 


0.86 


22.4 


0.56 


MRS (A') 


0.50Sf? ax 


20.4 


0.67 


28.5 


0.24 


MRST 


0.50^^ 


25.3 


0.39 


29.6 


0.20 


MRST( 5 |) 


0.50£f? ax 


21.6 


0.60 


30.1 


0.18 


MRST( 5 j) 


0.50£^ ax 


47.5 


0.003 


47.9 


0.003 


CTEQ3M 


0.25S?? ax 


21.4 


0.61 


28.1 


0.26 


CTEQ3M 


o^osf? 3 * 


25.3 


0.39 


32.7 


0.11 


CTEQ3M 


0.75^^ 


25.8 


0.37 


32.5 


0.11 


CTEQ3M 


l.OOE^ 13 * 


24.8 


0.42 


31.7 


0.14 



eks predictions for 0.1 < \r] \ < 0.7 all give \ 2 values with probabilities below 10% for the 
choices examined. 



I. Comparison with Previously Published Results 



The top panel in Fig. |55| shows (D — T)/T for our data in the 0.1 < \rj\ < 0.7 region 
relative to a jetrad calculation using the CTEQ4HJ PDF, fx = 0.5E™ ax , and 1Z sep = 1.3. 
Also shown are the previously published CDF data from the 1992-1993 Fermilab Tevatron 
running period relative to the same jetrad calculation. For this rapidity region, we 
have carried out a x 2 comparison between our data and the nominal curve describing the 
central values of the data of Ref. ]/]]. Comparing our data to the nominal curve, as though 
it were theory, we obtain a x 2 °f 56.5 for 24 degrees of freedom (probability of 0.02%). 
Thus our data cannot be described with this parameterization. As illustrated in the middle 



panel of Fig. [55], our data and the curve differ at low and high Et; such differences cannot 
be accommodated by the highly correlated uncertainties of our data. If we include the 
systematic uncertainties of the data of Ref. in the covariance matrix, the \ 2 is reduced 
to 30.8 (probability of 16%), representing acceptable agreement. 



J. Rapidity Dependence of the Inclusive Jet Cross Section 

D0 has subsequently extended the measurement of the inclusive jet cross section as a 
function of Et to \ t)\ < 3 in several bins of pseudorapidity |26[]. In this analysis the details 



of the jet energy scale corrections, single jet resolutions, and vertex selection were updated 
to minimize uncertainties for jets at large pseudorapidity (\rj\ > 1.5). These cross sections 
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TABLE XIX. x 2 comparisons between EKS and the data for \rj\ < 0.5 and 0.1 < \ < 0.7 
with n = DEf 3 * or DE^\ TZ scp = 1.3, and various PDFs. There are 24 degrees of freedom. 



PDF 


D where 




hi 


< 0.5 






0.1 < \r)\ 


< 0.7 






H = DE$ XX 






^max 




Ef 








x 2 


Prob. 


x 2 


Prob. 


x 2 


Prob. 


x 2 


Prob. 


CTEQ3M 


0.25 


21.1 


0.63 


17.9 


0.81 


32.3 


0.12 






CTEQ3M 


0.50 


20.7 


0.66 


19.3 


0.74 


33.7 


0.09 


33.3 


0.10 


CTEQ3M 


0.75 


20.4 


0.67 


19.4 


0.73 


33.3 


0.10 


33.0 


0.10 


CTEQ3M 


1.00 


20.2 


0.68 


19.5 


0.73 


32.9 


0.11 


32.7 


0.11 


CTEQ3M 


1.25 


20.4 


0.68 


19.8 


0.71 


32.8 


0.11 


32.8 


0.11 


CTEQ3M 


1.50 


20.8 


0.65 


20.3 


0.68 


33.1 


0.10 


33.1 


0.10 


CTEQ3M 


1.75 


21.5 


0.61 


21.2 


0.63 


33.5 


0.09 


33.6 


0.09 


CTEQ3M 


2.00 


22.4 


0.55 


22.1 


0.57 


34.2 


0.08 


34.3 


0.08 


CTEQ4M 


0.50 


19.4 


0.73 


18.2 


0.80 


33.8 


0.09 


34.3 


0.08 


CTEQ4HJ 


0.50 






23.3 


0.50 










CTEQ4A1 


0.50 






18.4 


0.78 










CTEQ4A2 


0.50 






18.3 


0.79 










CTEQ4A4 


0.50 






18.4 


0.78 










CTEQ4A5 


0.50 






19.2 


0.74 










MRS(A') 


0.50 






19.3 


0.74 






36.8 


0.05 
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FIG. 55. Top: Normalized comparisons of our data and of the CDF data in Ref. [7] to the 
jetrad prediction using CTEQ4HJ with \x = O.SE 1 ™^. Middle: Difference between our data and 
smoothed results of CDF normalized to the latter. The shaded region represents the ±lcr systematic 
uncertainties about the D0 data. The dashed curves show the ±la systematic uncertainties 
about the smoothed CDF data. Bottom: A comparison of the systematic uncertainties of the D0 
measurement and the CDF measurement. 
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are compared with jetrad predictions generated with \i 



0.5£^ ax , 



scp 



1.3, and similar 



choices of PDF given in Table [XVII]| . The data and theory were also compared using the 
same \ 2 test as used in this paper (Section |XI H| ) . The data indicate an preference for the 



CTEQ4HJ, MRST(gT), and CTEQ4M PDFs 26 



K. Conclusions 

We have made the most precise measurement to date of the inclusive jet cross section 
for Ex > 60 GeV at yfs = 1800 GeV. No excess production of high-.&r jets is observed. 
QCD predictions are in good agreement with the observed cross section for most standard 
parton distribution functions and different renormalization scales (fi = 0.25-2.00-Er where 
E T = E^ x or E^). 
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XII. RATIO OF INCLUSIVE JET CROSS SECTIONS AT ,/s 

GeV 



1800 AND 630 



A. Inclusive Jet Cross Section at y/s = 630 GeV 



The inclusive jet cross section for \rj \ < 0.5 at y/s = 630 GeV consists of data collected 
with three triggers: Jet_12, Jet_2_12, and Jet_30. To form the inclusive jet cross section, 
an £t region of each trigger is selected to maximize statistical power while maintaining full 
trigger efficiency. Any given cross section bin receives contributions from one and only one 
trigger. The luminosity in any given bin is the luminosity exposure for that trigger (given 
in Table 0). 

The inclusive jet cross section at y/s = 1800 GeV was determined prior to the 630 GeV 
analysis. To facilitate the ratio calculation as a function of xt = 2E^/y/s, the bin boundaries 
for the 630 GeV analysis were selected such that 



7-1630 

Ulrp 



630 
1800 



• E, 



1800 



;i2.i) 



i.e., such that the bin edges match in x^-space. Most of the resulting bins are 3.5 GeV 
wide, but some bins have a width of 7.0 GeV, 10.5 GeV, or more. 

Figure |56] displays the observed cross section at y/s = 630 GeV. The three different 
symbols indicate the E^ region for each jet trigger. Vertical lines (mostly hidden by the 
symbols) indicate the statistical uncertainty on each point. 

The cross section is corrected for the effects of jet resolution using the same method 
as used for the y/s = 1800 GeV cross section (Section [XI E|) . The single jet resolutions at 
y/s = 630 GeV are given in Section |XD|. The resulting ansatz fit parameters are given in 



Table XX and the unsmearing correction is plotted in Fig. 57. 



TABLE XX. Unsmearing ansatz function parameters (see Eq. 11 .3| ) for the inclusive jet cross 
section for \r] \ < 0.5 (in nb) at y/s = 630 and 1800 GeV and their uncertainties. 



yfs (GeV) 


Parameter 


Value 






Error Matrix 










A 


23.43 


4.48825 x 10" 


-3 


-1.15352 x 10" 


-3 


-6.70083 x 10" 


-3 


1800 


a 


-5.04 


-1.15352 x 10" 


-3 


2.98882 x 10" 


-4 


1.79044 x 10" 


-3 







8.23 


-6.70083 x 10" 


-3 


1.79044 x 10" 


-3 


1.21005 x 10" 


-2 




A 


22.7 


2.28649 x 10" 


-2 


-8.62822 x 10" 


-3 


-5.76500 x 10" 


-2 


630 


a 


-5.33 


-8.62822 x 10" 


-3 


3.28592 x 10" 


-3 


2.24650 x 10" 


-2 




(3 


6.58 


-5.76500 x 10" 


-2 


2.24650 x 10" 


-2 


0.16449 x 10" 


-1 



The resulting inclusive jet cross section at y/s = 630 GeV is given in Table |XXI| and is 
plotted in Fig. |58|. The uncertainties in the cross section are given in Table [XXII] and are 
also plotted in Fig. |59l. The bin-to-bin correlations of the uncertainties are shown in Fig. BO 



and are given in Table |XXII I 



The magnitude of the energy scale uncertainties are larger for the cross section at y/s = 
630 GeV than at y/s = 1800 GeV (Table |XI]| ). This is caused by several factors. The cross 
section at 630 GeV begins with jet E T > 20 GeV compared with 60 GeV at y/s = 1800 GeV. 
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FIG. 56. The observed inclusive jet cross section for y/s = 630 GeV. Symbols indicate the three 
jet triggers (shaded triangles: Jet_12; hollow diamonds: Jet_2_12; shaded circles: Jet_30). 



l 




C3 0.6 - 



05 1 1 ' 1 ' 1 ' 1 ■ 1 

0.1 0.2 0.3 0.4 0.5 

FIG. 57. The nominal unsmearing correction at ^/s = 630 GeV is given by the central line. The 
outer curves depict the uncertainty in the unsmearing due to the uncertainties in the measurement 
of the resolution of the jet energy. 
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The uncertainty in the energy scale offset correction (which is additive) has a much larger 
effect on 20 GeV jets than on 60 GeV jets. For E T > 60 GeV the cross section at y/s = 
630 GeV is much steeper than the cross section at 1800 GeV, hence the same uncertainty 
in the energy scale will lead to a larger uncertainty in the cross section. 

Figures |61] and ^ show the fractional differences between the data and several jetrad 
predictions using different choices of renormalization scale and PDF. These NLO QCD 
predictions are in reasonable agreement with the data. The data and predictions are com- 
pared quantitatively using a x 2 test (Section |XI H|) . The resulting x 2 values are given in 
Table X.X1V ; each comparison has 20 degrees of freedom. All but two of the jetrad predic- 
tions adequately describe the cross section at y/s = 630 GeV. For these, the probabilities for 
X 2 to exceed the listed values are between 6.4% and 78%. The prediction using MRST(g'T) 
and /i = 0.5i?™ ax produces the highest probability. The predictions using MRST(g j) with 
/i = 0.5.E™ ax , and CTEQ3M with /i = 2i?™ ax , thus are inconsistent with our measurements 
with probabilities < 0.4%. 
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FIG. 58. The \ rj\ < 0.5 inclusive cross section at y/s = 630 GeV. Statistical uncertainties are 
not visible on this scale (except for the last point). The histogram represents the jetrad prediction 
and the shaded band represents the ±l<r systematic uncertainty band about the prediction. 



B. The Ratio of Jet Cross Sections 

The dimensionless inclusive jet cross section (Section [IV A 2 ) is given by 

T 



El d 2 a 



a 



(12.2) 



90 





TABLE XXI 


Inclusive jet 


cross section for 


\rj \ < 0.5 at ys = 630 Gey. 




Bin Range 


Plotted 


Plotted 


Cross Section 


Systematic 








XT 


± statistical error 


Uncertainty 


(GeV) 


(GeV) 




(nb) 


(%) 


91 n 


— 94 ^ 


22.6 


0.07 


(2.56 ±0.03) x 10 2 


21.7 


94 ^ 


— 98 n 


26.1 


0.08 


(1.07 ±0.02) x 10 2 


17.2 


98 n 


Ol.J 


29.6 


0.09 


(5.14 ±0.16) x 10 1 


14.6 


31 \ 


— 3^ n 

oo.u 


33.1 


0.11 


(2.67 ±0.05) x 10 1 


13.0 


oo.u 


— 38 ^ 
oo.o 


36.7 


0.12 


(1.37 ±0.04) x 10 1 


12.1 


38 ^ 
oo.o 


— 49 n 


40.2 


0.13 


(7.96 ±0.27) x 10° 


11.5 


49 n 


— 4^ ^ 


43.7 


0.14 


(4.24 ±0.20) x 10° 


11.2 


4^ ^ 

4:0.0 


— 4Q n 


47.2 


0.15 


(2.83 ±0.16) x 10° 


11.0 


4Q n 


— ^9 ^ 
oz.o 


50.7 


0.16 


(1.81 ±0.13) x 10° 


10.9 


^9 ^ 
oz.o 


OO.U 


54.2 


0.17 


(1.14 ±0.03) x 10° 


10.9 


56 n 


- 59 5 


57.7 


0.18 


(7.35 ±0.21) x lO" 1 


11.0 


59.5 


63.0 


61.2 


0.19 


(5.07 ±0.17) x lO" 1 


11.1 


63.0 


66.5 


64.7 


0.21 


(3.29 ±0.14) x 10- 1 


11.3 


66.5 


70.0 


68.2 


0.22 


(2.42 ±0.12) x lO" 1 


11.5 


70.0 


73.5 


71.7 


0.23 


(1.64 ±0.10) x 10" 1 


11.8 


73.5 


77.0 


75.2 


0.24 


(1.18 ±0.08) x lO" 1 


12.1 


77.0 


80.5 


78.7 


0.25 


(8.79 ± 0.72) x 10~ 2 


12.4 


80.5 


94.5 


85.2 


0.27 


(3.69 ± 0.23) x 10~ 2 


13.6 


94.5 


- 112.0 


100.5 


0.32 


(1.05 ±0.11) x 10~ 2 


16.2 


112.0 


- 196.0 


136.2 


0.43 


(5.81 ± 1.19) x 10~ 4 


20.4 



TABLE XXII. Percentage cross section uncertainties for \ rj\ < 0.5 at ^/s = 630 GeV. The last 
row gives the nature of the bin-to-bin xt correlations: signifies uncorrelated uncertainties; 1 



correlated; and p partially correlated. 



X T 


Statistical 


Jet Selection 


Luminosity 


Trigger 


Unsmcaring 


Energy Scale 


Total 


0.07 


1.3 


0.2 


4.4 


2.4 


8.4 


19.4 


21.7 


0.08 


2.1 


0.2 


4.4 


0.9 


5.5 


15.6 


17.2 


0.09 


3.2 


0.2 


4.4 


0.3 


3.8 


13.3 


14.6 


0.11 


1.8 


0.2 


4.4 


0.6 


2.8 


11.9 


13.0 


0.12 


2.5 


0.2 


4.4 


0.2 


2.2 


11.0 


12.1 


0.13 


3.4 


0.2 


4.4 


0.1 


1.8 


10.5 


11.5 


0.14 


4.7 


0.1 


4.4 


0.0 


1.6 


10.2 


11.2 


0.15 


5.8 


0.1 


4.4 


0.0 


1.4 


9.9 


11.0 


0.16 


7.3 


0.1 


4.4 


0.0 


1.3 


9.8 


10.9 


0.17 


2.2 


0.1 


4.4 


1.1 


1.3 


9.9 


10.9 


0.18 


2.8 


0.1 


4.4 


0.4 


1.2 


10.0 


11.0 


0.19 


3.4 


0.1 


4.4 


0.2 


1.2 


10.1 


11.1 


0.21 


4.2 


0.2 


4.4 


0.1 


1.2 


10.3 


11.3 


0.22 


4.9 


0.2 


4.4 


0.0 


1.2 


10.6 


11.5 


0.23 


6.0 


0.2 


4.4 


0.0 


1.2 


10.8 


11.8 


0.24 


7.0 


0.2 


4.4 


0.0 


1.2 


11.2 


12.1 


0.25 


8.2 


0.2 


4.4 


0.0 


1.2 


11.6 


12.4 


0.27 


6.3 


0.2 


4.4 


0.0 


1.2 


12.8 


13.6 


0.32 


10.6 


0.3 


4.4 


0.0 


1.3 


15.5 


16.2 


0.43 


20.4 


0.3 


4.4 


0.0 


1.5 


19.8 


20.4 


Correlation 








1 





1 


P 


P 



91 



0.3 




Total Systematic Uncertainty 






★ 


Energy Scale 








• 


Luminosity 










Resolution 








T 


Trigger Efficiency 






0.2 


*\ 

A 






^* 


0.1 


*>~^_^^r*^ 
********* 








* 











* 

T T ' S '' !! ' ,!!,, 5 1 '5'**##*##* # * 
T 'TUT 'TTTTTT T T 

I.I.I. 




• 

* 

T 




0.1 


0.2 0.3 


0.4 





Jet x T 

FIG. 59. Contributions to the | rj\ < 0.5 at = 630 GeV cross section uncertainty plotted by 
component. 




FIG. 60. The correlations for the total systematic uncertainty for the inclusive jet cross section 
for \ r]\ < 0.5 at y/s = 630 GeV (Table |XXIII| ). 
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TABLE XXIII. The correlations for the total systematic uncertainty for the inclusive jet cross 

section for \rj\ < 0.5 at \fs = 630 GeV and the ratio of cross sections. The correlation values 

above the diagonal are for the cross section at y/s = 630 GeV and the values below the diagonal 
correspond to the ratio of cross sections. In both cases the correlation matrices are symmetric. 

1.00 0.99 0.96 0.91 0.86 0.82 0.79 0.74 0.69 0.66 0.62 0.60 0.58 0.56 0.54 0.52 0.50 

1.00 0.99 0.96 0.93 0.90 0.87 0.83 0.79 0.76 0.73 0.70 0.69 0.67 0.65 0.63 0.61 

1 1.00 1.00 0.99 0.97 0.95 0.93 0.90 0.87 0.84 0.82 0.79 0.78 0.77 0.75 0.73 0.72 

2 0.98 1.00 1.00 0.99 0.98 0.97 0.94 0.92 0.90 0.88 0.86 0.85 0.84 0.83 0.81 0.80 

3 0.95 0.99 1.00 1.00 1.00 0.99 0.97 0.96 0.94 0.93 0.91 0.90 0.89 0.88 0.87 0.85 

4 0.90 0.95 0.98 1.00 1.00 1.00 0.99 0.97 0.96 0.95 0.94 0.93 0.92 0.91 0.90 0.89 

5 0.84 0.91 0.95 0.99 1.00 1.00 0.99 0.99 0.98 0.97 0.96 0.95 0.94 0.94 0.93 0.92 

6 0.79 0.87 0.93 0.97 0.99 1.00 1.00 1.00 0.99 0.98 0.98 0.97 0.97 0.96 0.95 0.94 

7 0.74 0.83 0.89 0.95 0.98 0.99 1.00 1.00 1.00 0.99 0.99 0.98 0.98 0.97 0.97 0.96 

8 0.68 0.77 0.84 0.91 0.94 0.96 0.97 1.00 1.00 1.00 0.99 0.99 0.99 0.98 0.98 0.97 

9 0.62 0.72 0.80 0.87 0.91 0.94 0.96 0.99 1.00 1.00 1.00 1.00 0.99 0.99 0.99 0.98 

10 0.57 0.67 0.76 0.84 0.89 0.92 0.95 0.98 0.99 1.00 1.00 1.00 1.00 0.99 0.99 0.99 

11 0.52 0.63 0.72 0.80 0.86 0.90 0.93 0.97 0.98 0.99 1.00 1.00 1.00 1.00 0.99 0.99 

12 0.47 0.58 0.68 0.76 0.83 0.87 0.90 0.95 0.97 0.98 0.99 1.00 1.00 1.00 1.00 0.99 

13 0.43 0.54 0.64 0.73 0.80 0.84 0.88 0.93 0.95 0.97 0.99 0.99 1.00 1.00 1.00 1.00 

14 0.39 0.50 0.60 0.70 0.76 0.82 0.85 0.91 0.94 0.96 0.98 0.99 0.99 1.00 1.00 1.00 

15 0.34 0.46 0.56 0.66 0.73 0.78 0.83 0.88 0.91 0.94 0.96 0.98 0.99 0.99 1.00 1.00 

16 0.30 0.41 0.52 0.62 0.69 0.75 0.79 0.85 0.89 0.92 0.95 0.97 0.98 0.99 0.99 1.00 

17 0.25 0.37 0.47 0.57 0.65 0.71 0.76 0.82 0.86 0.89 0.93 0.95 0.96 0.98 0.99 0.99 1.00 

18 0.14 0.25 0.36 0.46 0.55 0.61 0.66 0.73 0.78 0.82 0.86 0.90 0.92 0.94 0.96 0.97 0.98 1.00 

19 -0.01 0.09 0.20 0.30 0.38 0.45 0.50 0.58 0.63 0.69 0.74 0.79 0.82 0.85 0.88 0.91 0.93 0.97 

20 -0.13 -0.04 0.06 0.15 0.24 0.30 0.36 0.44 0.50 0.56 0.62 0.68 0.72 0.75 0.79 0.83 0.86 0.93 



0.44 


0.39 


0.36 


1 


0.56 


0.51 


0.48 


2 


0.67 


0.62 


0.59 


3 


0.76 


0.71 


0.68 


4 


0.82 


0.77 


0.74 


5 


0.85 


0.81 


0.78 


6 


0.88 


0.85 


0.82 


7 


0.91 


0.88 


0.85 


8 


0.94 


0.91 


0.88 


9 


0.95 


0.93 


0.90 


10 


0.97 


0.94 


0.92 


11 


0.97 


0.95 


0.93 


12 


0.98 


0.96 


0.94 


13 


0.98 


0.97 


0.95 


14 


0.99 


0.97 


0.95 


15 


0.99 


0.98 


0.96 


16 


0.99 


0.98 


0.97 


17 


1.00 


0.99 


0.98 


18 




1.00 


0.99 


19 






1.00 


20 


1.00 








0.98 


1.00 







TABLE XXIV. x 2 comparisons of the inclusive jet cross section for 1 77 1 < 0.5 at \fs = 630 GeV 
with several theoretical predictions (20 degrees of freedom). 



PDF 


A* 


x 2 


prob. 






40.5 


0.4% 


CTEQ3M 


H/rp 


25.9 


17% 




0.5£f? ax 


30.4 


6.4% 




0.25£f? ax 


27.5 


12% 


CTEQ4M 


0.5Kp ax 


24.2 


24% 


CTEQ4HJ 




19.0 


53% 


MRST 


0.5£^ ax 


22.6 


31% 


MRST( 5 |) 


0.5£f ax 


14.9 


78% 


MRST( 5 j) 


0.5£f? ax 


51.8 


0.01% 



93 



o 



o 
H 



Q 



0.4 



0.2 - 



-0.2 - 



-0.4 



★ D0Data 

JETRAD: CTEQ3M, n = 0.5E™ x R sep = 1.3 








- 


—s K 






•ill 1 ' 1 ' ' 
:"f.J..........................——"" 






|l = 0.25E™ X 

- \l= l.OOE™" 

-■- H = 2.00E™ X 

i.i.i 



50 100 150 

E T (GeV) 

FIG. 61. The difference between the data and the prediction (jetrad), divided by the predic- 
tion for 1 7] | < 0.5 at yfs = 630 GeV. The solid stars represent the comparison to the calculation 
using CTEQ3M with /i = 0.5E^ iax . The shaded region represents the ±lcx systematic uncertainties 
about the prediction. The effects of changing the renormalization scale are also shown (each curve 
shows the difference between the alternative prediction and the standard prediction). 
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FIG. 62. The difference between the data and the prediction (jetrad), divided by the predic- 
tion for 1 77 1 < 0.5. The solid stars represent the comparison to the calculation using fi = 0.5-E™ ax 
and the PDFs CTEQ4M, CTEQ4HJ, MRST, MRST( 5 |), and MRST( 5 |). The shaded region 
represents the ±lcr systematic uncertainties about the prediction. 
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where d 2 a /dE^dr] is given by Eq. |1 1 . 1| . The ratio of inclusive jet cross sections for \rj \ < 0.5 
is calculated in bins of identical Xt- 

R(x T ) = a ^ T \ . (12.3) 

O"l800 [XT) 



C. Uncertainties in the Ratio of Jet Cross Sections 

Most of the systematic uncertainties in the inclusive jet cross section are highly corre- 
lated as a function of Et and center-of-mass energy, and cancel when the ratio of the two 
cross sections is calculated. To determine the uncertainty in the ratio, all uncertainties are 
separated into three categories, depending on the correlation (p) as a function of E T and 
CM energies. In most cases, complete correlation in E^ at one CM energy implies complete 
correlation between CM energies, but exceptions exist and are highlighted in the following 
sections. 



1. Luminosity Uncertainties 

The luminosity calculation at ^/s = 630 GeV shares many common uncertainties with 
the calculation at 1800 GeV (Section |V1]). The uncertainty from the fit to the world average 
(WA) pp total cross section determines the uncertainty in the luminosity at y/s = 630 GeV 
(Fig. fUj) . A lcr shift in the mean value of the cross section at y/s = 1800 GeV directly 
impacts the central value of the cross section at y/s = 630 GeV, resulting in a shift of 
unequal magnitude but like direction. The magnitude of the shift at 630 GeV, subtracted 
in quadrature from the interpolation uncertainty, defines two uncertainty components: the 
shift, which is completely correlated with the 1800 GeV cross section uncertainty, and the 
remainder, which is added in quadrature with the other independent luminosity uncertain- 
ties. The uncertainty components in the WA elastic and single-diffractive pp cross sections 



are handled with the same procedure. Table |XXV| lists the systematic uncertainties due to 
the luminosity for the ratio. 

TABLE XXV. The uncertainties in the ratio of cross sections due to the luminosity calculation. 
*Includes trigger matching uncertainty. 

Source Uncertainty (Percent) 
World Average pp cross section ±3.2 
Hardware efficiency ±3.6 
Geometric acceptance ±0.8 
Uncorrelated* ±2.6 



All sources* 



±4.2 
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2. Jet and Event Selection Uncertainties 

At 1800 GeV, the total uncertainty for jet cut efficiencies, the $ T cut efficiency, and 
the vertex cut efficiency is 1%. An independent study at 630 GeV determined cut uncer- 
tainties that were smaller (Table |XXVI|) . Despite some similarities in methodology, these 
uncertainties are all considered to be independent of one another in the ratio. 



TABLE XXVI. Uncertainty from jet and event selection. 





Uncertainty Source 


Uncertainty 


1800 GeV 


all selection 


1% below 350 GeV 




cut efficiencies 


2% above 350 GeV 




Jet cuts 


0.12 to 0.53% 


630 GeV 




0.03% 




vertex cut 


0.006% 



3. Resolution and Unsmearing Uncertainties 

Uncertainty in the unsmearing correction is dominated by the uncertainty in the jet res- 
olution measurement. In the case of ^/s = 1800 GeV, the systematic uncertainty dominates; 
for a/s = 630 GeV, poor statistics result in a fit uncertainty that is larger in magnitude than 
the systematic uncertainty. The systematic uncertainties in the unsmearing correction are 
assumed to be uncorrelated between the CM energies, as are the fitting uncertainties. The 
magnitudes of the resolution and unsmearing uncertainties are illustrated in Fig. 



4- Energy Scale Uncertainties 

The uncertainty in the inclusive jet cross section and in the ratio of cross sections is 
calculated using a Monte Carlo Simulation. The event generator performs several steps for 
each y/s and each cross section bin in xt- First, it generates a sample of jets with an xt 
spectrum which matches that observed in data. Second, it closely imitates true running 
conditions by simulating luminosity, vertexing, and smearing effects; thus the energy scale 
corrections of each Monte Carlo jet will closely match the corrections in real data. Third, the 
uncertainties from the energy scale corrections are calculated. Finally, the weighted average 
uncertainties and correlations in each bin are combined to form a covariance matrix. 

The jet E T distribution must be identical to the observed (smeared) jet cross section in 
data. The Monte Carlo simulation: 

1. Randomly generates the initial parton momenta X\ and x<i- 

2. Generates the corresponding and other kinematic quantities for both of the final- 
state partons (which result in jets). 
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FIG. 63. The uncertainty components in the ratio of inclusive jet cross sections as a function 
of xt plotted by component. 



3. Smears the jets according to the known resolution functions and then selects one jet 
at random. 

4. Checks that the selected jet falls within the desired xt bin and has | rf Gt \ < 0.5 (or 
starts over). 

5. Generates a weight for the jet, to reproduce the steeply falling spectrum of the inclusive 
jet cross section, using either a theoretical weight based on CTEQ4M and the scale of 
the collision, or an experimental weight based on the ansatz from unsmearing. 

Because the generated jet distribution already represents the energy-scale-corrected jet 
.Bp, and because the response correlation is given in terms of the energy before the response 
correction, the energy scale algorithm must be run "in reverse" to find the uncertainties and 
their correlations as a function of jet and CM energy. 

The ratio of inclusive jet cross sections is given in Eq. |12.2| . The elements of the covariance 
matrix are 

Cy = (Pij SRi SRj) , (12.4) 

where p expresses the correlation between the xt bins i and j, and the uncertainties in the 
ratio 5R may be expressed as 



(12.5) 
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where the two partial derivatives possess opposite signs: 

9Ri - f racial 800 - ^ 



da} 800 (a 4 1800 ) 2 ^ 1800 ' 1 ' 

Defining x = xt, the dependence of 6a on jet energy is given by: 

da a 2 da a 

Sat = -^Sxf = - sin Oi ■ 5 El (12.7) 

oxi a oxi 

The cross section uncertainty is now expressed in terms of jet energy, the jet angle, the CM 
energy (a), and the slope of the dimensionless cross section. The final expression for the 
covariance matrix elements becomes 



2 a 2 ^ a R k dot Rt da\ 



b 



C ij = Q~ sin ®k t sin 6i — 



0,6 k,i a b a k dx k °i dxi 

pt 5E a k 6E$) , (12.8) 

where a and b indicate CM energies; is the correlation between the uncertainties of the 
two jets whose energies fall in bins k and I, originating from the data sets at ^/s = a and 
b; and q is a factor that accounts for the negative sign in Eq. |12.6|: q = 1 when a = b, and 



q — — 1 otherwise. The bracket notation indicates the average. The summations indicate 
the four relevant correlations, visually described in Fig. R3. 



1800 




Ratio in 
Bin i 

FIG. 64. Correlations between two ratio bins i and j. Arrows indicate the four possible corre- 



lation (p) terms. The uppermost arrow is p\ 



630-630 

ij i 



while the 



arrow is p 



630-1800 

ij 



As mentioned previously, interpolation of the correlation matrix determines the values of 
Pm for the response uncertainty. For the completely correlated uncertainties, all p's take the 
value of unity; for the uncorrelated uncertainties, all p's are zero. The major contribution 
originates from the partially correlated response uncertainty. 



5. Combined Uncertainty in the Ratio 



The individual uncertainties of the earlier sections fall into several classifications, summa- 
rized in Table [XXV11| . Complete cancellation of uncertainties occurs when the uncertainties 
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are completely correlated between CM energies. The components of the systematic uncer- 



tainties for the ratio of cross sections are plotted in Fig. ^ and given in Table |XXVIII| . The 



uncertainty in the energy scale correction dominates at each end of the spectrum; resolution 
and contributions from other sources (primarily the luminosity uncertainty) become impor- 
tant at intermediate values of xt- Figure |65| plots the point-to-point uncertainty correlations 
between data points. 



TABLE XXVII. Uncertainty correlations in the ratio of cross sections. "0" indicates no corre- 



lation, "1" indicates complete correlation. 



Uncertainty Source 


Correlation 


in 

Jet E T 


Comments 


Luminosity 


partial 


1 




Filter Match 





1 


1800 GeV Only 


Event Cuts 










Jet Cuts 










Resolution 








Fits 


partial 


1 




Closure 


1 


1 




Unsmearing Fits 





1 




Energy Scale 








Offset 


partial 


1 




Response fit 


1 


partial 




Response at 630 GeV 





1 




Showering 


1 


1 





D. Results and Comparison to Theoretical Predictions 



The ratio between the inclusive jet cross sections at yfs = 630 and 1800 GeV is given 
in Table |XXIX| . Figures ^6] and B7| show the ratios of cross sections compared with jetrad 
predictions using different choices of PDF and renormalization scale. The measured ratios 
lie approximately 10% below the theoretical predictions, which have an uncertainty of ap- 
proximately 10% (Section [LV A 2 ). Table XXX lists the x 2 distributions for the ratio of cross 
sections compared to selected theoretical predictions. The \ 2 values lie in the range 15.1-24 
for 20 degrees of freedom (corresponding to probabilities in the range 28% to 77%). The 
best agreement occurs for extreme choices of renormalization scales: fj, = (0.25, 2.00)i?™ ax . 
As expected, there is very little dependence on the choice of PDF. 

Different renormalization scales can be selected for the different CM energies since there 
is no explicit theoretical need for identical scales at = 630 and 1800 GeV. Figure |68| 
depicts a comparison between the ratio and theoretical predictions where the renormalization 
scales at the two CM energies are not equivalent. The resulting x 2 indicate good agreement 
between the data and the predictions (Table [XXXlp . 
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FIG. 65. The correlation matrix for the ratio of cross sections. Axes indicate the bin numbers. 



TABLE XXVIII. Percentage uncertainties in the ratio of inclusive jet cross sections at y/s = 630 
and 1800 GeV for \ n\ < 0.5. 



X T 


Statistical 


Jet Selection 


Trigger Match 


Luminosity 


Trigger Efficiency 


Unsmearing 


Energy Scale 


Total 


0.07 


1.5 


1.1 


1.7 


4.1 


2.4 


6.6 


9.5 


12.7 


0.08 


2.4 


1.1 


1.7 


4.1 


0.9 


4.2 


7.4 


9.7 


0.09 


3.5 


1.1 


1.7 


4.1 


0.3 


2.7 


6.0 


8.0 


0.11 


1.9 


1.1 


1.1 


4.1 


0.6 


1.9 


5.0 


7.0 


0.12 


2.7 


1.1 


1.1 


4.1 


0.2 


1.3 


4.4 


6.3 


0.13 


3.5 


1.1 


1.1 


4.1 


0.1 


1.0 


4.0 


6.0 


0.14 


4.9 


1.1 


1.1 


4.1 


0.0 


0.8 


3.7 


5.8 


0.15 


5.8 


1.1 


0.0 


4.1 


0.0 


0.7 


3.4 


5.5 


0.16 


7.3 


1.1 


0.0 


4.1 


0.0 


0.7 


3.2 


5.4 


0.17 


2.4 


1.1 


0.0 


4.1 


1.1 


0.7 


3.2 


5.4 


0.18 


3.0 


1.1 


0.0 


4.1 


0.4 


0.7 


3.2 


5.4 


0.19 


3.5 


1.1 


0.0 


4.1 


0.2 


0.7 


3.3 


5.4 


0.21 


4.4 


1.1 


0.0 


4.1 


0.1 


0.8 


3.4 


5.5 


0.22 


5.1 


1.1 


0.0 


4.1 


0.0 


0.8 


3.5 


5.5 


0.23 


6.2 


1.1 


0.0 


4.1 


0.0 


0.8 


3.7 


5.6 


0.24 


7.3 


1.1 


0.0 


4.1 


0.0 


0.8 


3.9 


5.8 


0.25 


8.5 


1.1 


0.0 


4.1 


0.0 


0.9 


4.1 


5.9 


0.27 


6.6 


1.1 


0.0 


4.1 


0.0 


1.0 


4.7 


6.4 


0.32 


11.0 


1.1 


0.0 


4.1 


0.0 


1.1 


6.3 


7.7 


0.43 


20.5 


1.1 


0.0 


4.1 


0.0 


1.3 


8.7 


9.7 
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1ABLE XXIX. 


lire ratio of the inclusive 


jet cross sections for \ i] < 0.5 at ^ 


Js = 630 and 1800 


GeV. 










M 


Ratio of Cross Sections 


Systematic 


Bin Range 


Plotted 


± statistical error 


Uncertainty (%) 


0.067 - 0.078 


0.072 


1.72 ±0.03 


12.7 


0.078 - 0.089 


0.083 


1.64 ±0.04 


9.7 


0.089 - 0.100 


0.094 


1.62 ±0.06 


8.0 


0.100 - 0.111 


0.105 


1.67 ±0.03 


7.0 


0.111 - 0.122 


0.116 


1.57 ±0.04 


6.3 


0.122 - 0.133 


0.127 


1.59 ±0.06 


6.0 


0.133 - 0.144 


0.139 


1.48 ±0.07 


5.8 


0.144 - 0.156 


0.150 


1.63 ±0.09 


5.5 


0.156 - 0.167 


0.161 


1.64 ±0.12 


5.4 


0.167 - 0.178 


0.172 


1.64 ±0.04 


5.4 


0.178 - 0.189 


0.183 


1.62 ±0.05 


5.4 


0.189 - 0.200 


0.194 


1.67 ±0.06 


5.4 


0.200 - 0.211 


0.205 


1.60 ±0.07 


5.5 


0.211 - 0.222 


0.216 


1.74 ±0.09 


5.5 


0.222 - 0.233 


0.228 


1.69 ±0.10 


5.6 


0.233 - 0.244 


0.239 


1.78 ±0.13 


5.8 


0.244 - 0.256 


0.250 


1.81 ±0.15 


5.9 


0.256 - 0.300 


0.271 


1.74 ±0.11 


6.4 


0.300 - 0.356 


0.319 


1.85 ±0.20 


7.7 


0.356 - 0.622 


0.432 


1.83 ± 0.38 


9.7 



TABLE XXX. The calculated x 2 for the ratio of cross sections (20 degrees of freedom). 



PDF Renormalization Scale X 2 Prob. 

2£™ ax 17^9 60% 

£3p ax 21.6 36% 

CTEQ3M 0.75^ ax 23.1 28% 

0.5Ef ax 20.5 43% 

0.25£^ ax 15T 77% 

CTEQ4M 0.5££ iax 22A 32% 

CTEQ4HJ 0.5£f? ax 21.0 40% 

MRST 0.5^ ax 22.2 33% 

MRST(c/T) 0.5Ef ax 19.5 49% 

MRST(^|) 0.5£f? ax 24.1 24% 
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FIG. 66. The ratio of dimensionless cross sections for 1 77 1 < 0.5 compared with jetrad predic- 
tions with fx = 0.5££ iax and the CTEQ3M, CTEQ4M, CTEQ4HJ, and MRST PDFs. The shaded 
band represents the ±la systematic uncertainty band about the prediction. 



TABLE XXXI. x 2 comparisons for the ratio of cross sections for \rj\ < 0.5 where the renor- 
malization scale is mismatched between CM energies. 



PDF 


Renormalization 
630 GeV 


Scale 
1800 GeV 


x 2 


Prob. 




2^ max 


0.5£^ ax 


14.9 


78% 


CTEQ3M 


£max 


0.5Ef ax 


17.2 


64% 




0.25£^ ax 




23.1 


28% 
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FIG. 67. The ratio of dimensionless cross sections for \rj\ < 0.5 compared with jetrad pre- 
dictions with various values of fi and the CTEQ3M PDF. The shaded band represents the ±la 
systematic uncertainty band about the prediction. 
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FIG. 68. The ratio of dimensionless cross sections for \rj\ < 0.5 compared with jetrad pre- 
dictions with n = 0.5£f? ax at yfe = 1800 GeV, /z = (0.25, 1.0, 2.0)Ef ax at yfs = 630 GeV, and 
the CTEQ3M PDF. The shaded band represents the ±lcr systematic uncertainty band about the 
prediction. 
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An additional analysis was carried out to measure the significance of the normalization 
difference between the data and the theoretical predictions. The data are reduced to a 
single value by fitting them to a constant (horizontal line), resulting in a value of 1.60 ± 
0.08. The uncertainty in this value is given by the uncertainty of the fit. Each of the 
theoretical predictions is also reduced to a single value. Each x T point of the prediction 
is assigned a weight given by the statistical uncertainty of the corresponding point in the 
data (Table |XXX1J| ). The uncertainty in the value representing the theoretical prediction 
is assumed to be zero. The resulting x 2 values are given in Table pCXXIIj , and lie in the 
range 1.4-13.2 (corresponding to probabilities of 0.03% to 23%). In every case, discarding 
the shape information in favor of a comparison of normalization results in poorer agreement 
between data and the theory. 



TABLE XXXII. Normalization-only predictions for the ratio of cross sections and the x 2 com- 
parison with the data (1.60 ± 0.08 for one degree of freedom). 



PDF 


Renormalization Scale 


Theory Normalization 


x 2 


Prob. 






1.75 


3.3 


6.8% 




pmas 


1.82 


7.1 


0.8% 


CTEQ3M 


Q.lhE^ 


1.87 


10.7 


0.1% 




0.5 Ef™ 


1.85 


9.5 


0.2% 




0.25£f? ax 


1.70 


1.5 


22.9% 


CTEQ4M 


0.5£5p ax 


1.90 


13.2 


0.03% 


CTEQ4HJ 


O.SEf™ 


1.87 


10.7 


0.1% 


MRST 


OXEf™ 


1.89 


12.6 


0.04% 


MRST(#T) 


0.5 E^ K 


1.87 


11.1 


0.09% 


MRST( 5 j) 


0.5£f? ax 


1.90 


12.9 


0.03% 



630 GeV 



1800 GeV 



CTEQ3M 



Hirp 

0.25£f? ax 



0.5£f ax 
O.bE^ 



1.46 
1.71 
1.44 



2.7 
1.8 
3.7 



10.2% 
18.1% 
5.4% 



E. Conclusions 

We have made the most precise measurement to date of the ratio of the inclusive jet cross 
sections at ^fs = 630 and 1800 GeV. This measurement is nearly insensitive to the choice 
of parton distribution functions. The ratio of cross sections is therefore a more stringent 
test of QCD matrix elements. The NLO QCD predictions yield satisfactory agreement with 
the observed data for standard choices of renormalization scale or PDF. In terms of the 
normalization however, the absolute values of the standard predictions lie consistently and 
significantly higher than the data. 
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XIII. DUET ANGULAR DISTRIBUTION 

The dijet angular distribution is given by: 

1 d 3 a 1 N 



13.1) 



E cr dM dx dr/ hoost £ Ni AMAxAr) hoost 
where the invariant mass is calculated assuming massless jets: 

M 2 = 2E$ tl E$ t2 [cosh (At?) - cos (A0)] , (13.2) 

the pseudorapidity of the center-of-mass of the dijet system is given by 77* = ~(Ar]); the 
pseudorapidity boost is given by r] hoost = ±(m + 7/a); X = exp (| At] |) = exp (2|t?*|); E^ 1 , r) X , 
and 4>i refer to the values associated with the jet with the largest Ey in an event; E^ t2 , 772, 
and 02 refer to the values associated with the jet with the second largest E T in an event; 
Arj = 1 771 — 772 I ; Acf) = <pi — <p 2 , and N is the number of events in a given x an d mass bin. 
If the individual jet masses are taken into account, the change in the dijet invariant mass 
is less than 1% for jets used in this analysis. Since the bins of Arj^ oost are constant and we 
plot the angular distribution for a given mass bin, AM, we choose to measure da/dx which 
is uniform for Rutherford scattering) 

1 da 1 N4 , 

(13.3) 



A. Data Selection 

The selected data are events with two or more jets which satisfy the set of inclusive jet 
triggers and pass the standard jet and event quality requirements (Section |VII1|) . Events are 
removed unless both of the leading two jets pass the jet quality requirements. The vertex 
of the event must be within 50 cm of z = 0. The x distributions were corrected for the 
efficiencies of the standard jet quality cuts and the $ T cut. 

To ensure that the jet triggers did not introduce a bias, the trigger requirement was 
verified by comparing the x distribution of a lower trigger threshold to the x distribution of 
the desired trigger threshold. It is known that the lower threshold trigger is 100% efficient 
in the desired region and thus a comparison would show an inefficiency in the desired trigger 
sample. No differences were seen. The E^s of all second jets are well within the region 
of 100% jet reconstruction efficiency, so an additional Et requirement on the second jet 
was not necessary. The final energy-scale-corrected Et requirement placed on each trigger 



sample is summarized in Table XXXIU. 



B. Acceptance: Limits on Mass and x 



Event acceptance is calculated using the kinematic relationships between mass, x> an d Et 
shown in Fig. |69|. Since an Et requirement is placed only on the leading jet, the maximum 
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TABLE XXXIII. The cut on the E T of the leading jet to ensure that the trigger is 100% 
efficient. 



Trigger Corrected Ei Limit on Leading Jet(GeV) 
Jet_30 55^0 
Jet_50 90.0 
Jet_85 120.0 

Jet_115 175.0 



X with 100% acceptance is determined from the Et requirement placed on the leading jet 
and the desired mass bin using the following formula: 



M 



2E 2 



Tl 



[cosh (ln(x)) + 1]. 



(13.4) 



In this formula the Et's of the two leading jets are assumed to be identical. Four mass bins 
were chosen in order to maximize the number of events per \ bin, and to attain a maximum 
X of 20 (corresponding to rj* = 1.5). These mass bins are listed in Table [XXXI V. 



TABLE XXXIV. 


The average mass, 


maximum \ measured, 


and the number 


of events after 


applying all kinematic cuts. 








Trigger Et 


Mass 


Average 


Xmax 


Number 


Threshold 


Range 


Mass 




of Events 


(GeV) 


(GeV/c 2 ) 


(GeV/c 2 ) 






55 


260-425 


302 


20 


4621 


90 


425-475 


447 


20 


1573 


120 


475-635 


524 


13 


8789 


175 


>635 


700 


11 


1074 



Once the x limit is known, a limit on r] hoost can be calculated. The r] boost parameter is 
used to restrict the x distribution to the physical limits of the detector (Fig. |7Dp . The //boost 
limit is calculated using 

l^boostl = ||?7*| - \VmaxW (13.5) 

= ||1.5| - 1 3 .0 1 1 = 1.5, 

where |?7 m ax |= 3.0 is the maximum rj used for this analysis. The boost cut is chosen to be 
//boost < 1-5. For M > 475 GeV/c 2 , | 7] boost | is kinematically restricted to a value less than 
1.5. These mass bins are listed together with the average dijet invariant mass, the maximum 
X measured, and the number of events for each of four mass ranges in Table |XXX1 V . 
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E T > 175 GeV 
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475 < M <-'635 
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E T > 55 GeV 



260 < M < 425 



200 




FIG. 69. In the mass versus x plane, the curves shown are contours of constant Et- The 
simplest form of uniform acceptance in this plane is a rectangle. For a chosen mass region, the 
limit on \ corresponds to the intersection of the lower mass limit and the Et contour. The shaded 
regions shown are the mass bins chosen for this analysis. 
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acceptance 



FIG. 70. The 1 77 1 's of the two leading jets were required to be less than 3.0. For a maximum 
77* of 1.5, the |r?boost I was chosen to less than 1.5 to restrict the measurement to a region of 100% 
acceptance. 

C. Systematic Studies 

In order to study the systematic effects of the jet and event selection requirements, 
and the various corrections that are applied to the data, a series of systematic studies was 
performed. For each requirement or correction, we measured the effect on the angular 
distribution by varying the requirement or correction by an appropriate amount. 

To determine the systematic uncertainties on the shape of the angular distribution, each 
distribution is fit with a function: F = A + B\ + C/x + D/x 2 + E/x 3 - The effect of varying 
each of the selection criteria or corrections is measured by taking the ratio of the distribution 
with the nominal selection criteria, and with the adjusted criteria (Fig. |7l|) , giving the size 
of the systematic uncertainty. 

The largest source of uncertainty involves the r\ dependence of the jet energy scale. Small 
uncertainties in the relative response as a function of 77 have large effects on the angular 
distribution. The uncertainties in the jet energy scale are less than 2% up to an \r\ \ of 2.0 
and become large near ~ 3.0. The uncertainty in the showering correction is less than 
2% for 1 77 1 < 2.0 and becomes large at high \ rj\. The effect of the //-dependent energy scale 
uncertainties are given in Fig. |72|(a). 

The resolution of our measurement of the jet energy can also affect the angular distribu- 
tion. This was determined by measuring the difference between the smeared and unsmeared 
theory calculations. Since we are not unsmearing the data for the effects of rj and Ej- 
smearing, we apply this as an uncertainty in the measurement (Fig. 172(b)). 
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FIG. 71. Technique to determine changes in angular distribution due to systematic uncertain- 
ties. Top: Normalized angular distribution for the mass range 475-635 GeV/c 2 compared to a fit 
to the data. Middle: Normalized angular distribution for the mass range 475-635 GeV/c 2 after 
removing the jet quality cuts and a fit to the data. Bottom: The (ratio— 1) of the two fits shows 
the effects of the jet quality cuts on the shape of the angular distribution. 



The effect on the angular distribution due to the rj bias in the jet reconstruction algorithm 
(Section |III G| ) was studied by applying a correction for the bias. The difference between 
the corrected and uncorrected distributions was 1% on average (Fig. |72|(c)). 

The overall energy scale does not affect the shape of the distribution, because a shift 
in the overall energy scale shifts the entire distribution in mass. The angular distribution 
changes very slowly with mass, so a small shift would not cause a significant change in the 
shape. 

For the Jet_30 and Jet_50 triggers, an online mitool (see Section [VII] ) requirement was 
used in the trigger for part of the run. To determine if the mitool requirement biased the 
angular distribution, runs with no MITOOL requirement were compared to runs with the 
requirement. A small shape difference was seen and an uncertainty equal to the difference 
between the two measurements was assigned. 

The effects of multiple interactions on the distributions were studied. A secondary in- 
teraction adds approximately 0.6 GeV of per unit At] x A<p (Fig. |20|). Since the angular 
distribution is measured in regions in which the E^s of the two leading jets are in excess 
of 50 GeV and are often above 100 GeV, the effect of this additional energy on the two 
leading jets is minimal. It is possible that a second interaction may produce a vertex which 
is incorrectly used as in the primary vertex for the leading two jets. This would cause an 
error in the measured r\ positions of the jets as well as the measured Et of the jets. We 
studied the effect of not selecting the primary vertex by minimizing the St in the event 
(Section [VI]|) . This has a negligible effect on the angular distribution. 
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It is possible that the vertex produced by a second interaction is the only vertex found 
in the event. This would also cause an error in the measured 77 and E? values of the jets. 
We studied the possibility of multiple interactions affecting the angular distribution in this 
manner by the following method. For a determined percentage of events, we switched the 
vertex to a randomly chosen vertex. The new vertex was based on the measured vertex dis- 
tribution, which has an approximate mean of z = and a a ~ 30 cm. We then recalculated 
the 77 and E? of the two leading jets in the event and measured the angular distribution. 
The percentage of events with a new vertex was determined based on the efficiency of vertex 
reconstruction for events with large Et jets (~ 70%), and the percentage of multiple inter- 
actions in the data used for this analysis (~ 60%). The number of vertices switched was 
20%, which is an estimate of the number of times that the vertex reconstruction is incorrect. 
The size of the effect is less than 2% and is dependent on the value of x ( see Fig- |72](d)). 

The jet quality requirements and their corresponding efficiency corrections are necessary 
to remove noise from the event sample. Their effect on the shape of the angular distribution 
is minimal. 

The D0 jet algorithm allows for the splitting and merging of jets. This can cause a shift 
in the 77 of the jet, and therefore affect the angular distribution. The effect on the shape of 
the distribution of removing those events in which either of the leading two jets were split 
or merged is minimal. Since the theoretical predictions are expected to properly address 
merging and splitting, no uncertainty was assigned due to this effect. 
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0.05 
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X 

FIG. 72. Ratios of parameterized curves showing the effects of uncertainties on the shape of the 
angular distribution. Shown are the four largest uncertainties in the mass range 475-635 GeV/c 2 . 
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D. Results and Comparison to Theory 



The measurement of the dijet angular cross section is given in Table |XXXV| . The leading- 
order and next-to-leading-order theory predictions were obtained using the jetrad parton- 
level event generator |3| with CTEQ3M and \i = 0.5£^ ax . Four mass ranges are compared 
to the LO and NLO predictions of QCD in Fig. [73|. The band at the bottom of each plot 
represents the ± 1 a systematic uncertainties. They are obtained by adding in quadrature 
all of the parameterized curves describing the shape uncertainties discussed earlier. 



TABLE XXXV. Dijet angular cross section (100/A r )(dA r /dx) ± statistical ± systematic un- 
certainties for the four mass bins (GeV/c 2 ). 



X 


260< M <425 




425< M <475 




475 < M <635 




M >635 




value ± 


stat.i 


sys. 


value ± 


stat.rt 


sys. 


value ± 


stat.i 


sys. 


value ± 


stat.± 


sys. 


1.5 


5.95 ± 


0.35± 


0.58 


7.58 ± 


0.66± 


2.08 


10.08 ± 


0.33± 


0.63 


11.98 ± 


0.99± 


0.49 


2.5 


5.50 ± 


0.33± 


0.54 


4.26 ± 


0.50± 


0.75 


7.56 ± 


0.28± 


0.36 


12.49 ± 


l.Olrb 


0.78 


3.5 


4.59 ± 


0.30± 


0.31 


4.96 ± 


0.53± 


0.67 


7.83 ± 


0.29± 


0.33 


9.11 ± 


0.86± 


0.61 


4.5 


4.57 ± 


0.30± 


0.28 


5.54 ± 


0.56± 


1.04 


7.71 ± 


0.28± 


0.25 


9.79 ± 


0.89± 


0.23 


5.5 


4.56 ± 


0.30± 


0.25 


5.29 ± 


0.55± 


0.86 


7.87 ± 


0.29± 


0.17 


10.06 ± 


0.91± 


0.26 


6.5 


5.10 ± 


0.32± 


0.23 


6.26 ± 


0.60± 


0.73 


8.17 ± 


0.29± 


0.16 


9.58 ± 


0.88± 


0.51 


7.5 


5.10 ± 


0.32± 


0.19 


4.83 ± 


0.53± 


0.33 


8.70 ± 


0.30± 


0.20 


9.30 ± 


0.87± 


0.57 


8.5 


5.61 ± 


0.34± 


0.15 


4.40 ± 


0.50± 


0.16 


7.91 ± 


0.29± 


0.21 


8.08 ± 


0.81± 


0.42 


9.5 


4.93 ± 


0.32± 


0.09 


5.60 ± 


0.57± 


0.25 


8.46 ± 


0.30± 


0.24 


8.96 ± 


0.85± 


0.30 


10.5 


6.04 ± 


0.35± 


0.06 


5.22 ± 


0.55± 


0.37 


8.62 ± 


0.30± 


0.27 


10.65 ± 


0.93± 


0.41 


11.5 


5.40 ± 


0.33± 


0.04 


4.30 ± 


0.50± 


0.40 


8.38 ± 


0.30± 


0.29 








12.5 


5.33 ± 


0.33± 


0.08 


4.75 ± 


0.52± 


0.52 


8.69 ± 


0.30± 


0.36 








13.5 


5.41 ± 


0.33± 


0.14 


5.43 ± 


0.56± 


0.65 














14.5 


5.40 ± 


0.33± 


0.20 


5.69 ± 


0.57± 


0.70 














15.5 


5.60 ± 


0.34± 


0.28 


6.18 ± 


0.60± 


0.76 














16.5 


4.81 ± 


0.31± 


0.30 


4.70 ± 


0.52± 


0.57 














17.5 


4.95 ± 


0.32± 


0.38 


4.83 ± 


0.53± 


0.56 














18.5 


5.78 ± 


0.34± 


0.53 


5.01 ± 


0.54± 


0.55 














19.5 


5.37 ± 


0.33± 


0.57 


5.17 ± 


0.55± 


0.55 















Also shown in Fig. [73 



are comparisons to NLO theory predictions calculated using 
/i = _E™ ax . With the large angular reach measured, the angular distribution is sensitive to 
the choice of renormalization scale. The QCD theoretical predictions are in good agreement 
with the measured angular distributions. 



E. Compositeness Limits 

A comparison to theory is made to test for quark compositeness (Section |IV B|) . Predic- 
tions of the theory of compositeness are available at LO. In order to simulate NLO prediction 
with compositeness, we generated LO curves at various values of A. We measured the frac- 
tional differences between the LO angular distribution with A = oo and those with finite A 
values. We then multiplied the NLO prediction of the angular distribution by these fractional 



differences. The results are shown in Fig. 74 for the mass bin with M > 635 GeV/c 2 . 

To remove the point-to-point correlated uncertainties, the distribution can be character- 
ized by a single number: R x = N(x < X)/N(X < x < Xmax), the ratio of the number of 
events with x < X to the number of events between X < x < Xmax (where Xmax is given 
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FIG. 73. Dijet angular distributions for D0 data (points) compared to JETRAD for LO (dashed 
line) and NLO predictions with renormalization/factorization scale fi = O^I?™^ (dotted line). The 
data are also compared to jetrad NLO predictions with /x = E™ x (solid line). The errors on 
the data points are statistical only. The band at the bottom represents the ± 1 a systematic 
uncertainty. 
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FIG. 74. Data compared to theory for different compositeness scales. See text for an expla- 
nation of the compositeness calculation. The errors on the points are statistical and the band 
represents the ± 1 a systematic uncertainty. 



in Table |XXXIV|) . The choice of the value of X in the definition of R x is arbitrary. To 
optimize the choice of X, the following study was performed. The percentage change in the 
largest mass bin between NLO QCD and NLO QCD with a contact term of A^ L = 2.0 TeV 
was measured as a function of the definition of R x . The change due to compositeness in- 
creases as one chooses smaller values of X. However, the measurement error also increases 
for smaller values of X. Only the statistical error was used to optimize the choice of X so 
as not to bias the optimization with the data. The ratio of percentage change to percentage 
statistical error peaks at X = 4; hence we chose R x = N(x < 4)/iV(4 < \ < Xmax)- 

To determine the errors on R x , the nominal value was compared to the value after 
each systematic uncertainty was varied within error. Table |XXXVI| shows the size of the 
systematic uncertainties for the smallest and largest mass ranges. 

Table XXXVI l| shows the experimental ratio R x for the different mass ranges with their 
statistical and their systematic uncertainties, which are fully correlated in mass. Figure [F5 
shows R x as a function of M for two different renormalization scales, along with the theo- 
retical predictions for different compositeness scales. The effects of compositeness should be 
greatest at the highest masses. Note that the two largest dijet invariant mass bins have a 
lower Xmax value (Table |XXX1V| ), and thus a higher value of R x is expected independent of 
compositeness assumptions. Also shown in Fig. ^ are the \ 2 values for the four degrees of 
freedom for different values of the compositeness scale. 

The method chosen to obtain a compositeness limit uses Bayesian statistics |5rj |. The 
compositeness limit is determined using a Gaussian likelihood function for R x as a function 
of dijet mass. The likelihood function is defined as 



115 



TABLE XXXVI. The systematic uncertainties on the measurement of R x for the smallest and 
largest mass ranges. 





260-425 GeV/c 2 


Mass Range 


> 635 GeV/c 2 


Misvertexing 


0.24 




0.001 


MITOOL 


0.0076 




0.000 


Jet quality cuts 


0.002 




0.010 


rj bias 


0.007 




0.009 


Energy scale 


0.01 




0.023 


Resolution 


0.004 




0.010 



TABLE XXXVII. The dijet angular ratio R x and its statistical and systematic uncertainty. 
Also listed are the jetrad predictions with TZ scp = 1.3, the CTEQ3M PDF, and fi = 0.5Ef ax and 
Ef ax . 



Mass Range 




Theory 




GeV/c 2 


R x ± Stat. ± Sys. 


H = 0.5E^ ax 


H = Ef ax 


260-425 


0.191 ±0.0077 ±0.015 


0.198 


0.180 


425-475 


0.202 ±0.0136 ±0.010 


0.206 


0.185 


475-635 


0.342 ± 0.0085 ± 0.018 


0.342 


0.344 


>635 


0.506 ±0.0324 ±0.028 


0.506 


0.458 



116 



0.7 
0.6 

0.5 
0.4 

0.3 
0.2 



0.1 



|l=0.5E™ x 
•D0 Data 



A 
O 



A = 1.8 TeV 

LL 

A + = 2.0 TeV 

LL 

NLO 



=17.8 
=10.9 
=6.3 

L_ 



H=E! 



A 

O 



A = 1.8 TeV 

LL 



10.2 



X 

L 

A + = 2.0TeV y 2 =5.7 

LL 

NLO Z 2 =3-9 
I 



250 



500 



750 250 

M (GeV/c 2 ) 



500 



750 



FIG. 75. R x as a function of dijet invariant mass for two different renormalization scales. The 
inner error bars are the statistical uncertainties and the outer error bars include the statistical and 
systematic uncertainties added in quadrature. The \ 2 values for the four degrees of freedom are 
shown for the different values of the compositeness scale. The data are plotted at the average mass 
for each mass range. The NLO points are offset in mass to allow the data points to be seen. 



\S\2u 



exp 



2 L 



m\ T v- i [d-f{i)\\p{i) 



;i3.6) 



where d is a 4 component vector of data points for the different mass bins, /(£) is a 4 
component vector of theory points for the different mass bins for different values of £ where 
£ is related to A (see below) , V" 1 is the inverse of the covariance matrix, and P(£) is the 
prior probability distribution, P(£). The covariance matrix is defined so that the element 

is 



i,j of the covariance matrix, Vi 



i j ■ 



V i4 



AcTj • Act 



;i3.7) 



where Act is the sum of the systematic and statistical uncertainties added in quadrature 
if i = j, and the systematic uncertainty only if i ^ j. The systematic uncertainties are 
assumed to be 100% correlated as a function of mass. 

The compositeness limit depends on the choice of the prior probability distribution, 
P(£). Motivated by the form of the Lagrangian, P(£) is assumed to be flat in £ = 1/A 2 . 
Since the dijet angular distribution at NLO is sensitive to the renormalization scale, each 
renormalization scale is treated as a different theory. To determine the 95% confidence level 
(C.L.) limit in A, a limit in £ is first calculated by requiring that Q(£) = Jq L(R x \£')d£' = 
0.95Q(oo). The limit in £ is then transformed into a limit in A. Table [XXXVH"I| shows the 
95% C.L. limit for the compositeness scale obtained for different choices of models. These 



results supersede those reported in |23[ following the correction of an error in the program 
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used to calculate the effects of compositeness in that paper. The resulting limits are reduced 
by approximately 150 GeV. If the prior distribution is assumed to be flat in 1/A 4 , the limits 



are slightly reduced, as shown in Table XXXIX 



TABLE XXXVIII. The 95% confidence level limits for the left-handed contact compositeness 
scale for different models. The prior probability distribution is assumed to be flat in 1/A 2 . 



Compositeness scale 


H = £™ ax 


H = 0.5Ef ax 


A+ 

^LL{ud) 


2.0 TeV 
2.0 TeV 
1.8 TeV 
1.8 TeV 


2.1 TeV 

2.2 TeV 
1.8 TeV 
2.0 TeV 


TABLE XXXIX. 
scale for different mc 


The 95% confidence level limits for the left-handed contact compositeness 
dels. The prior probability distribution is assumed to be flat in 1/A 4 . 


Compositeness scale 


H = Ef ax 


fj, = 0.5Ef ax 


All 

All 
A+ 

^LL(ud) 


2.0 TeV 
1.9 TeV 
1.7 TeV 
1.7 TeV 


2.0 TeV 

2.1 TeV 
1.7 TeV 
1.9 TeV 


Recently published results from CDF |19| on dijet angular distributions 
model in which all quarks are composite, yielding 95% confidence limits A£ L 


compare to the 
> 1.8 TeV and 



h~ LL > 1.6 TeV. 



F. Coloron Limits 



Predictions of the dijet angular distribution with colorons are available at LO (Sec- 
tion IV C ). To simulate NLO predictions, coloron LO predictions are generated for sev- 



eral values of M c j cot 9. The fractional differences between the angular distribution with 
M c j cot 9 = oo and the distributions with finite values of M c j cot 9 are measured. The col- 
oron NLO predictions are then obtained by multiplying the NLO QCD prediction obtained 
using jetrad by the LO fractional differences obtained above. The results are shown in 



Fig. 76 



Limits on the coloron mass are calculated using the same method as in the previous 
section. For a renormalization scale of \x = i?™ ax , the 95% CL. limit on the coloron mass 
is M c /cot9 > 759 GeV/c 2 . If fi = 0.5£^ ax , the 95% CL. limit is MJ cot 9 > 786 GeV/c 2 . 



The resulting limits are shown in Fig. Tn\. The shaded region shows the 95% CL. exclu- 



sion region for the DO dijet angular distribution measurement (M c / cot# > 759 GeV/c 2 ). 
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FIG. 76. D0 data compared to theory for different values of M c j cot 6 (see text for details of the 
coloron distribution calculation). The errors on the points are statistical and the band represents 
the correlated ±lcr systematic uncertainty. 



The horizontally-hatched region at large cot 9 is excluded by the model [|38| , f40|1 . The 
diagonally-hatched region is excluded by the value of the weak-interaction p parameter 
(M c / cot# > 450 GeV/c 2 ) JEJ. The cross-hatched region is excluded by the CDF search for 



new particles decaying to dijets |f41 |. These limits are then converted into more general limits 
on color-octet vector current-current interactions: A^ g > 2.1 TeV assuming a s (Mz) = 0.12 
(Section IIVBD . 



G. Conclusions 



We have measured the dijet angular distribution over a large angular range. The data 
distributions are in good agreement with NLO QCD predictions. The compositeness limits 
depend on the choice of the renormalization/factorization scale, the model of compositeness, 
and the choice of the prior probability function. Models of quark compositeness with a 
contact interaction scale of less than 2 TeV are ruled out at the 95% C.L. 
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FIG. 77. Limits on the coloron parameter space: coloron mass M c vs. mixing parameter cot 9. 
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XIV. THE INCLUSIVE DUET MASS SPECTRUM 



The dijet mass spectrum is calculated using the relation: 



a 



dMdrjidri2 A^AMA^A^ 



14.1) 



where iVj is the number of events in mass bin i; Ci is the unsmearing correction; Li is the 
integrated luminosity; t; L is the efficiency of the trigger, vertex selection, and the jet quality 
cuts; AM is the width of the mass bin; and Az/i^ are the widths of the pseudorapidity bins 
for jets 1 and 2. The dijet mass is calculated assuming massless jets using Eq. [L3.2| . If we 
define the mass using four vectors, m 2 = {E x + E 2 ) 2 — (pi + P2) 2 , the cross section changes 
by less than 2%. 

A. Data Selection 

The selected data are events with two or more jets which satisfy the set of inclusive jet 
triggers and pass the standard jet and event quality requirements (see Section |V111| ). Events 
are removed unless both of the two leading jets pass the jet quality requirements. The vertex 
of the event must be within 50 cm of z = 0. The efficiency for each event is then given by the 
product of the efficiencies (e^) of the jet quality cuts, the efficiency (e mct ) of the cut on $ T , 
the efficiency (e tr ig g er) for an event to pass the trigger, and the efficiency (e verte x) for passing 
the vertex cut. The reciprocal of the resulting efficiencies (the event weights) is plotted as 
a function of M in Fig. [78|. The efficiency of the vertex requirement is 90 ± 1%. The data 
are used to select a sample where both jets have pseudorapidity \rf ct \ < 1.0. To examine 
the inclusive dijet cross section more closely, two sub-samples are created where both jets 
satisfy either | rf et | < 0.5 or 0.5 < | rf ct | < 1.0. 
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FIG. 78. The reciprocal of the event efficiencies €j for each mass bin. 



To determine the mass at which a given trigger (Jet_XX) becomes fully efficient, the 
event efficiencies are plotted for each of the triggers in the chosen mass ranges (Table KL|) 
in Figs. [79| and |8(| This plot shows that the triggers are > 95% efficient for most of the 
data. 
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The mass spectra obtained from the triggers Jet_30, Jet_50 and Jet_85 are then scaled 
to match the Jet_115 mass spectrum by the scale factor Sput which is given by 

**, = x (14.2) 

\ *-Jet_115 / e vertexi 

The values of S-put used to scale the data in this analysis are given in Table |XL| . These scales 
(and the event weights) are then applied to the data to produce the mass spectra (two such 
spectra are depicted in Fig. RTJ) . The error plotted for each point is given by the statistical 
errors for that bin. 



TABLE XL. The mass ranges and scale factors for the triggers used in this analysis. 





Data Sets satisfying \if et \ < 1.0, and 0.5 < \ rf et 


< 1.0. 


Trigger Name 


Mass Range 


Scaling Factor 




(GeV/c 2 ) 


Spilt 


Jet_30 


200-270 


289.3 ± 14.4 


Jet_50 


270-350 


21.7 ± 1.1 


Jet_85 


350-550 


1.845 ± 0.005 


Jet_115 


550-1400 


1.095 ±0.009 



Data Set satisfying \ rf et \ < 0.5. 



Trigger Name 


Mass Range 


Scaling Factor 




(GeV/c 2 ) 


Stilt 


Jet_30 


150-200 


289.3 ± 14.4 


Jet_50 


200-300 


21.7 ±1.1 


Jet_85 


300-390 


1.845 ± 0.005 


Jet_115 


390-1400 


1.095 ± 0.009 



B. Vertex Selection Biases 



The vertex selection procedure chooses the vertex with the smallest value of St (Sec- 
tion | VI I]) . This selection criterion may be biased for events where both of the two leading 



jets have the same absolute rapidity. In this case the vertex chosen would be the one that 
minimizes the for both of the leading jets and not necessarily the correct one. This 



bias was studied using the pythia |57[ MC event generator to generate events with multiple 
vertices at the same rate as the Jet_85 and Jet_115 triggers. For dijet events with | rf et | < 1.0 
the number of incorrectly chosen vertices is 5%. 

The effect on the dijet mass cross section is measured by calculating the ratio of the 
mass spectrum produced using the selected vertex to that of the correct vertex. The result 



of this calculation is given in Fig. |82| for | rf et \ < 1.0 and shows that the effect is of the order 



of 1% and that it is reasonably uniform as a function of mass. A 2% uncertainty in the cross 
section, uncorrelated as a function of mass, was assumed. 
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FIG. 79. The trigger efficiencies of the events included in the dijet mass spectrum. Note that 
most events have an efficiency greater than 99%. 
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FIG. 80. The average trigger efficiencies for each trigger as a function of mass. 
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FIG. 81. The scaled dijet mass spectrum for \ if ct \ < 1.0 and \ rf et \ < 0.5. 
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FIG. 82. The effect on the mass spectrum of incorrectly identified vertices for \rf ct | < 1.0. 
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C. Dijet Mass Resolution 

The dijet mass resolutions were calculated using the measured single jet resolutions 
(Section |X|). The dijet mass resolutions depend on the Et and r\ values of the two leading 
Et jets in each event. Hence the mass resolutions are determined by using a Monte Carlo 
(MC) event generator to convolute the measured single jet resolutions (Table |V|). For each 
MC event generated, the individual particle jets are smeared by the measured single jet 
resolutions. The unsmeared and smeared dijet masses are calculated and used to determine 
the mass smearing. The values of the mass smearing are plotted in discrete mass bins and 



fitted to a Gaussian ansatz (see Fig. |83|for an example). The distribution is well-represented 
by a Gaussian with only a small fraction (<C 1%) of events forming a tail (due to events 
where the jets are reordered after smearing). The resolution at each of these masses is given 
by the width of the Gaussian. The results obtained for \ 7f et \ < 1.0 using the pythia MC 
are plotted in Fig. R3. 




0.8 1 1.2 

Mass Resolution 

FIG. 83. The distribution of M smeared /M unsmcared for PYTHiA-generated events with 
\rf et \ < 1.0 and 490 < M unsmeared < 510 GeV/c 2 . 
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FIG. 84. The mass resolutions for | r/ jet | < 1.0 generated using the pythia MC. The solid curve 
and the MC data points show the resolutions determined using the nominal jet energy resolutions; 
the dashed lines show the resolutions determined with the ±lcr jet energy resolution uncertainties. 
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The mass resolutions are then fitted using the functional form: 



a{M)/M(%) = y/A + BM + CM 2 + DM 3 . 



14.3) 



The results are depicted in Fig. IS3. The resulting fit parameters for all r] regions considered 



in this analysis are given in Table |XLI 



TABLE XLI. 

PYTHIA MC. 


Parameterizations 


of the mass resolut 


ions (in percent) 


generated using the 


Data Set 


A 


B 


C 


D 


|r/j ct | < 1.0 
|r?j ct | < 0.5 
0.5 < \ rji ct \ < 1.0 


3.40 ± 1.01 
3.78 ± 0.94 
5.24 ± 0.83 


0.761 ± 0.045 
0.701 ± 0.041 
0.709 ± 0.058 


0.0302 ± 0.0032 
0.0231 ± 0.0025 
0.0389 ± 0.0022 


0.0002 ± 0.0005 
0.0 ± 0.0003 
0.0 ± 0.0004 



The mass resolution dependence on the MC generator used to convolute the single jet 



resolutions has been estimated by using the herwig and jetrad event generators. The 
jetrad program is used at LO with renormalization scales of [i = (0.25, 0.5, 1.0)-E™ ax and 
H = (0.25, 0.5, 1.0) VI with CTEQ3M To ensure that the choice of PDF does not affect 
the resolutions, jetrad was run with the CTEQ3L and MRS(A') pi PDFs. 




200 400 600 800 1000 

M (GeV/c 2 ) 

FIG. 85. A comparison of the mass resolutions for | rf et | < 1.0 obtained by running pythia 
(solid line), herwig (dashed line), and jetrad (dotted line) with /i = O.SEj^ and the CTEQ3M 
PDF. The resolutions are divided by the nominal pythia resolutions. The upper and lower curves 
show the effect of the ±1<t uncertainties on the measured jet resolutions. 



126 



D. Data Unsmearing 



The jet energy scale corrects only the average response to a jet. The steeply falling dijet 
mass spectrum is distorted by the jet energy resolution and, to a negligible extent, by the rj 
resolution. The observed mass spectrum is corrected for resolution smearing by assuming a 
trial unsmeared spectrum 



F(M') 



N . M '~ a (l-¥L 



which is convoluted with the measured mass resolutions 



f(M) = J F(M')o (M' - M, M') dAf, 



14.4) 



;i4.s) 



such that the number of events in any mass bin i is given by integrating / over that bin. The 



data were fitted using a binned maximum likelihood method and the minuit package |58 
to determine the values of N, a, and /3. The smearing correction is given by 



JF(M)dM 
JJdM 



(14.6) 



To account for any uncertainties in the choice of trial function the data are fitted with two 
additional functions: 



F(M) = N ■ M~ a 
F(M) = N ■ M~ a exp 



-P 




Viooy 



Viooy 



;i4.7) 
;i4.s) 



The nominal smearing correction is given by the fit to the data using the trial function 
given in Eq. |14.4| and the obtained mass resolutions (Table |XL1| ). The resulting fit for 



njet 



< 1.0 has a 



10.3 for 13 degrees of freedom and is given in Fig. 



The 



magnitude of the correction is approximately 5% at 209 GeV/c 2 and drops to approximately 
2% at 500 GeV/c 2 , and then rises to 8% at 1 TeV/c 2 . The uncertainty in the smearing 
correction is obtained by fitting the data with each of the trial functions and all of the 
mass resolutions generated with the different MC generators. The error is given by the 
maximum and minimum corrections obtained for each mass bin and is approximately 2%. 



The resulting smearing corrections are shown in Fig. 37. 



E. Energy Scale Corrections 

The uncertainties in the dijet mass cross section due to the energy scale have also been 
determined using a Monte Carlo program. The MC program generates two initial state 
partons each with a uniform distribution in x (the fraction of the proton momentum carried 
by the parton). The kinematic quantities of the two jets that result from this interaction 
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FIG. 86. Top: The fit to the data for | if** \ < 1.0 using Eq. |14-4 Bottom: The residuals of the 
fit are plotted, (Data — Fit) /Fit. 
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FIG. 87. The smearing correction factor to be applied to the data for \ rf ,et \ < 1.0. The upper 
and lower curves represent the ±ler uncertainties of the smearing correction. 
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are determined by generating a random value of x from a uniform distribution. The jet E^ 
and rj are uniquely determined for the event; it is assumed that the two jets are back-to- 
back in <f). The event is accepted if it satisfies the requirement l^ 6 *! < 1.0. Each event is 
weighted by M~ 4 ' 75 x Pr(a;i) x Pr(a;2) where Pr(xi 5 2) is the probability of finding a parton 
with momentum fraction x in the CTEQ4M PDF |TJj]. The exponent, -4.75, was chosen to 
obtain a dijet mass spectrum with similar normalization and shape as the data. Finally, each 
of the resulting jets has its smeared by the measured single jet resolutions. Figure 
shows a comparison of the mass spectrum produced by the Monte Carlo and the data; the 
two are in reasonable agreement. The effect of changing the weight applied in the Monte 
Carlo has been studied. If the weight is changed to M~ 4 ' 5 or M~ 5 ' the resulting energy 
scale error changes by less than 1%. 
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FIG. 88. A comparison of the data and the mass spectrum produced by the Monte Carlo to 
study the energy scale uncertainties. The solid curves show the ±ler energy scale uncertainties. 
The dashed curves show the MC predictions where the weights are set to M -4 ' 5 (upper) and M~ 5,0 
(lower) . 

The energy scale uncertainties are calculated by generating a sample of MC jet events 
(in which the jets are fully corrected) and applying the inverse of jet energy scale correc- 
tion. This sample of uncorrected jets then have the nominal, high (nominal+la), and low 
(nominal— la) energy scale corrections applied. The error due to the energy scale is split 
into components: the uncorrelated error, the fully correlated error, a partially correlated 
error, and the error due to the showering correction. The resulting errors are plotted in 
Fig. along with errors obtained from fitting the data with the high and low energy scale 
corrections applied. The errors obtained by the two methods are in agreement. 

We calculated the correlation matrix for the partially correlated component of the energy 
scale uncertainty. The correlations have been calculated as a function of the jet energy 
(Section [IX Dp ; hence the relationship between the jet energy and the dijet mass spectrum 
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FIG. 89. The energy scale errors obtained from the Monte Carlo. The full circles show the 
total energy scale uncertainties, the open squares show the correlated error, the open triangles 
show the uncorrelated error, and the stars show the partially correlated error. The curve shows 
the energy scale uncertainties obtained from fitting the data with the high and low energy scale 
corrections applied. 



needs to be determined. The correlation matrix for the dijet mass spectrum can be written 

as 



dMi dMj 



(14.9) 



where 5X represents the shift in variable X due to a systematic error parameter a; i and j 
denote mass bins and a is the cross section. In the limit of massless jets, the dijet effective 
mass can be approximated by: 



Mi = ^2E iX E l2 {\ -cose il2 ) (14.10) 
where E a and E i2 are the energies, and 9n 2 is the angle between the jets for event i. Hence: 



dE n ^ El1 + dE i2 ^ El2 ^ dcos 6 il2 



5cos9 il2 ; 



(14.11) 



as we are only concerned with the energy scale, the angle error is ignored. Therefore: 

(14.12) 



dM t E l2 

— (1 - cos il2 ), 



dE n Mi 



SMi = —(E a 5E i2 + Ei 2 5Ea){l - cos0 a2 ) 



(14.13) 
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and 



MiMj 



cos 6^2) (1 — cos 6 



jl2) 



EiiEj 1 {5E i2 Ej 2 ) + E i2 Ej 1 (5E il Ej 2 ) + 

EiiEj 2 {SE i2 E jl ) + E i2 E j2 {5E il E jl )\ (14.14) 

Using this relationship, the correlations between jets due to the uncertainties in the jet 
energy scale can be translated into correlations between mass bins for the dijet mass cross 
section using the Monte Carlo. The resulting correlations are plotted for a selection of mass 



bins in Fig. 90. 
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FIG. 90. The correlations between mass bins for |r/ jet | < 1.0. The four plots show the mass 
correlations relative to four different mass bins: 200-220 GeV/c 2 (209 GeV/c 2 weighted average), 
300-320 GeV/c 2 (309 GeV/c 2 ), 470-510 GeV/c 2 (488 GeV/c 2 ), and 600-700 GeV/c 2 (639 GeV/c 2 ). 



F. Summary of Systematic Uncertainties 

In addition to the uncertainties in the luminosity, smearing correction, and the energy 
scale, there are uncertainties associated with the selection of the events that contribute to 
the data sample. These uncertainties are due to the jet quality cuts (see Eqs. [O, K% O , 
and H2I) as well as the procedure used to add hot cells back into jets (Section [VIII A| ). These 
uncertainties contribute a 1% uncorrelated uncertainty to the cross section. In addition, the 
uncertainty due to the unsmearing is assumed to be fully correlated as a function of the 
dijet mass in each event. 
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A complete description of the systematic uncertainties in the dijet mass cross section is 
given in Table [XLII| . The total systematic error in each mass bin is given by the sum in 
quadrature of these errors. The uncertainties are combined appropriately to obtain an overall 
correlation matrix for the bin-to-bin systematic uncertainties in the dijet mass spectrum. 



TABLE XLII. Common systematic errors on the cross section. 



Source 


Percentage Error 


Comment 


Jet Selection 


1. 


Statistical 
Uncorrelated 


Vertex Selection 


2. 


Systematic 
Uncorrelated 


Luminosity Scale 


5.8 


Systematic 
Fully Correlated 


Luminosity Match 




Systematic 


Jet_30 


4.9 


Statistical 


Jet_50 




Correlated for 
Jet_30 and 50 


Unsmearing Correction 


0.5-3.0 


Systematic 
Fully Correlated 


Energy Scale 


7.0-30.0 


Systematic 
Mixture 



G. Cross Section 



The dijet mass cross section is calculated using Eq. |14.1| for the pseudorapidity range 



rf et \ < 1.0 , in mass ranges starting at 200, 270, 350, and 550 GeV/c 2 , corresponding to 
the jet E T thresholds of 30, 50, 85, and 115 GeV. 

The cross section for the mass spectrum is plotted in Fig. 91, and given in Ta- 
ble |XLIII| . The data are plotted at the mass-weighted average of the fit function for each bin 
(/ MFdM I J FdM). The systematic uncertainties are dominated by the uncertainties in the 
jet energy scale, which are 7% (30%) for the 209 (873) GeV/c 2 mass bins. The bin-to-bin 
correlations of the uncertainties are shown in Fig. [92| and are given in Table [XL IV] [55 |. 
The dijet mass cross section measurement was then repeated for | rf ct | < 0.5, and 0.5 < 



if 



jet 



< 1.0. The resulting cross sections are given in Tables [XLV| and |XLVI . 
Most of the systematic uncertainties in the measurement of the inclusive dijet mass 
spectrum are highly correlated as a function of dijet mass and rj and to a good approximation, 
cancel when a ratio of two cross sections is made. For this reason the cross section ratio for 
the rapidity ranges |?]' ct | < 0.5 and 0.5 < \ rf et \ < 1.0 will be calculated: 



14.15) 





| 7yj e t 


< 0.5) 




k(0.5 < 




< 1.0) 
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TABLE XLIII. Dijet mass cross section for |r/j et | < 1.0. High (low) systematic uncertainties 
are the sum in quadrature of the uncertainties from the ±1<7 variations in the energy calibration, 
the unsmearing, the vertex corrections, luminosity matching, jet selection, and the uncertainty in 
the luminosity. Also included is the jetrad prediction with fj, = 0.5£vp ax , lZ scp = 1.3, and the 
CTEQ3M PDF. 





Mass Bin (GeV/c 2 ) 




Cross Section 


Systematic Error 


Theoretical 


Min. Max. 


Weighted 




± Statistical Error 


Low 


High 


Prediction 






Center 




(nb/(GeV/c 2 )/(Ar/) 2 ) 


(%) 


(%) 


(nb/(GeV/c 2 )/(A?7) 2 ) 


200 


220 


209.1 


918 


(3.66 ±0.12) x 10~ 2 


11.2 


11.7 


3.57 x 10~ 2 


220 


240 


229.2 


507 


(2.03 ±0.09) x 10~ 2 


11.1 


11.4 


2.12 x 10~ 2 


240 


270 


253.3 


419 


(1.13 ±0.06) x 10~ 2 


11.3 


11.6 


1.17 x 10~ 2 


270 


300 


283.4 


2944 


(5.98 ±0.11) x 10~ 3 


11.4 


11.8 


6.00 x 10~ 3 


300 


320 


309.3 


1123 


(3.43 ±0.10) x 10~ 3 


11.4 


12.0 


3.53 x 10- 3 


320 


350 


333.6 


1006 


(2.06 ±0.06) x 10~ 3 


11.8 


12.1 


2.17 x 10~ 3 


350 


390 


367.6 


8749 


(1.14 ±0.01) x 10~ 3 


10.9 


11.5 


1.15 x 10 -3 


390 


430 


407.8 


4323 


(5.66 ±0.09) x 10~ 4 


11.4 


12.0 


5.67 x 10~ 4 


430 


470 


447.9 


2137 


(2.80 ±0.06) x 10" 4 


11.8 


12.7 


2.92 x 10" 4 


470 


510 


488.0 


1210 


(1.59 ±0.05) x 10~ 4 


12.3 


13.4 


1.54 x 10" 4 


510 


550 


528.0 


646 


(8.47 ±0.33) x 10~ 5 


12.7 


14.2 


8.36 x 10~ 5 


550 


600 


572.0 


699 


(4.35 ±0.16) x 10~ 5 


13.3 


15.2 


4.31 x 10~ 5 


600 


700 


638.9 


542 


(1.68 ±0.07) x 10~ 5 


14.8 


17.1 


1.55 x 10~ 5 


700 


800 


739.2 


144 


(4.43 ± 0.37) x 10~ 6 


17.5 


20.7 


3.75 x 10~ 6 


800 


1400 


873.2 


46 


(2.32 ±0.34) x 10~ 7 


23.1 


28.9 


1.95 x 10~ 7 
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FIG. 91. Dijet mass cross section (Pa/dMdrndm for \rf et \ < 1.0. The D0 data are shown 
by the solid circles, with error bars representing the ±1<t statistical and systematic uncertainties 
added in quadrature (in most cases smaller than the symbol). The histogram represents the jetrad 
prediction. 




FIG. 92. The correlations between systematic uncertainties in bins of dijet mass (see Ta- 
bles |XTffl] and pCLIVI ) for \ rf ct \ < 1.0. The correlations are calculated using the average systematic 
uncertainty. The discontinuities arise from the uncorrelated errors (adjacent to correlations of 1.0) 
and to luminosity matching. 
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TABLE XLIV. The systematic error correlations for the dijet cross section for | rf ct \ < 1.0, and 
the ratio «;(|r/j ct | < 0.5)/k(0.5 < \ri^ ct \ < 1.0). The correlation values above the diagonal are the 
correlations corresponding to the cross section and the correlations below the diagonal correspond 
to the ratio. In both cases the correlation matrices are symmetric. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 






1.00 


0.88 


0.88 


0.87 


0.87 


0.87 


0.76 


0.74 


0.73 


0.72 


0.71 


0.69 


0.66 


0.61 


0.55 


1 






1.00 


0.89 


0.89 


0.89 


0.89 


0.78 


0.77 


0.76 


0.74 


0.73 


0.71 


0.68 


0.63 


0.56 


2 


1 


1.00 




1.00 


0.90 


0.90 


0.90 


0.79 


0.78 


0.77 


0.76 


0.75 


0.73 


0.70 


0.65 


0.59 


3 


2 


0.58 


1.00 




1.00 


0.90 


0.91 


0.80 


0.79 


0.79 


0.78 


0.77 


0.75 


0.72 


0.68 


0.61 


4 


3 


0.61 


0.60 


1.00 




1.00 


0.91 


0.81 


0.80 


0.80 


0.79 


0.78 


0.77 


0.74 


0.70 


0.64 


5 


4 


0.61 


0.59 


0.63 


1.00 




1.00 


0.82 


0.82 


0.81 


0.81 


0.80 


0.79 


0.76 


0.73 


0.67 


6 


5 


0.52 


0.51 


0.54 


0.54 


1.00 




1.00 


0.89 


0.90 


0.89 


0.89 


0.87 


0.85 


0.81 


0.76 


7 


6 


0.56 


0.54 


0.57 


0.58 


0.87 


1.00 




1.00 


0.90 


0.90 


0.89 


0.88 


0.87 


0.84 


0.78 


8 


7 


0.60 


0.58 


0.62 


0.63 


0.54 


0.58 


1.00 




1.00 


0.91 


0.91 


0.90 


0.89 


0.86 


0.82 


9 


8 


0.60 


0.58 


0.61 


0.62 


0.53 


0.57 


0.63 


1.00 




1.00 


0.91 


0.91 


0.90 


0.88 


0.84 


10 


9 


0.59 


0.58 


0.61 


0.62 


0.53 


0.57 


0.62 


0.63 


1.00 




1.00 


0.92 


0.91 


0.90 


0.86 


11 


10 


0.58 


0.57 


0.60 


0.61 


0.53 


0.56 


0.61 


0.62 


0.62 


1.00 




1.00 


0.92 


0.91 


0.89 


12 


11 


0.59 


0.58 


0.61 


0.63 


0.53 


0.58 


0.62 


0.63 


0.63 


0.62 


1.00 




1.00 


0.93 


0.92 


13 


12 


0.60 


0.58 


0.62 


0.63 


0.54 


0.58 


0.63 


0.62 


0.62 


0.61 


0.63 


1.00 




1.00 


0.95 


14 


13 


0.58 


0.56 


0.59 


0.61 


0.53 


0.56 


0.61 


0.60 


0.60 


0.59 


0.61 


0.62 


1.00 




1.00 


15 


14 


0.53 


0.52 


0.55 


0.57 


0.50 


0.53 


0.56 


0.56 


0.56 


0.55 


0.57 


0.58 


0.58 


1.00 






15 


0.51 


0.50 


0.53 


0.55 


0.48 


0.52 


0.54 


0.54 


0.54 


0.53 


0.55 


0.57 


0.57 


0.55 


1.00 
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TABLE XLV. Dijet mass cross section for | if ct \ < 0.5. High (low) systematic uncertainties 
are the sum in quadrature of the uncertainties from the ±1<7 variations in the energy calibration, 
the unsmearing, the vertex corrections, luminosity matching, jet selection, and the uncertainty in 
the luminosity. Also included is the jetrad prediction with [i = 0.5Ej} ax , TZ scp = 1-3, and the 
CTEQ3M PDF. 





Mass Bin (GeV/c 2 ) 


AT 

Ni 


Cross Section 


Systematic Error 


1 neoretical 


Min. Max. 


Weighted 




± statistical Juror 


Low 


High 


Prediction 






Center 




(nb/ (LreV 1 c )j \l\r\ ) ) 


(%) 


(%) 


(nb/ (GeV / c ) / \l\rj ) ) 


150 


160 


154.7 


467 


(1.46 ±0.07) x lO" 1 


10.8 


11.0 


1.46 x 10" 1 


160 


180 


168.9 


552 


(8.69 ±0.37) x 10- 2 


10.5 


10.5 


9.03 x 10~ 2 


180 


200 


189.0 


315 


(4.99 ±0.28) x 10~ 2 


10.5 


10.9 


4.91 x 10~ 2 


200 


220 


209.1 


2243 


(2.69 ±0.06) x 10~ 2 


10.3 


10.5 


2.81 x 10~ 2 


220 


240 


229.2 


1390 


(1.67 ±0.04) x 10~ 2 


10.2 


10.7 


1.68 x 10~ 2 


240 


270 


253.3 


1055 


(8.52 ±0.26) x 10~ 3 


10.5 


10.5 


9.35 x 10~ 3 


270 


300 


283.4 


550 


(4.47 ±0.19) x 10~ 3 


10.7 


10.8 


4.82 x 10~ 3 


300 


320 


309.3 


2671 


(2.78 ±0.05) x 10~ 3 


9.1 


10.0 


2.86 x 10~ 3 


320 


350 


333.6 


2434 


(1.69 ±0.03) x 10~ 3 


9.9 


9.9 


1.78 x 10~ 3 


350 


390 


367.7 


1823 


(9.50 ±0.22) x 10~ 4 


9.9 


10.3 


9.49 x 10~ 4 


390 


430 


407.8 


1459 


(4.52 ±0.12) x 10~ 4 


10.3 


10.8 


4.78 x 10" 4 


430 


470 


448.0 


831 


(2.58 ±0.09) x 10~ 4 


10.7 


11.4 


2.50 x 10" 4 


470 


510 


488.0 


480 


(1.49 ±0.07) x 10~ 4 


11.1 


12.1 


1.35 x 10~ 4 


510 


550 


528.1 


231 


(7.17 ±0.47) x 10~ 5 


11.7 


12.8 


7.43 x 10~ 5 


550 


600 


572.2 


156 


(3.87 ±0.31) x 10~ 5 


12.4 


13.7 


3.91 x 10~ 5 


600 


700 


639.4 


125 


(1.55 ±0.14) x 10~ 5 


13.8 


15.5 


1.46 x 10~ 5 


700 


800 


739.8 


30 


(3.71 ±0.68) x 10~ 6 


16.5 


19.3 


3.72 x 10~ 6 


800 


1400 


878.1 


14 


(2.86 ±0.77) x 10~ 7 


22.0 


27.8 


2.08 x 10~ 7 
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TABLE XLVI. Dijet mass cross section for 0.5 < \r]^ ct \ < 1.0. High (low) systematic un- 
certainties are the sum in quadrature of the uncertainties from the ±lcr variations in the energy 
calibration, the unsmearing, the vertex corrections, luminosity matching, jet selection, and the 
uncertainty in the luminosity. Also included is the jetrad prediction with jjl = O^E™*, lZ sep = 



1.3, and the CTEQ3M PDF. 





Mass Bin (GeV/c 2 ) 




Cross Section 


Systematic Error 


Theoretical 


Min. Max. 


Weighted 




± Statistical Error 


Low 


High 


Prediction 






Center 




(nb/(GeV/c 2 )/(Ar/) 2 ) 


(%) 


(%) 


(nb/(GeV/c 2 )/(A?7) 2 ) 


200 


220 


209.1 


275 


(4.39 ±0.26) x 10~ 2 


12.3 


12.7 


4.56 x 10~ 2 


220 


240 


229.1 


170 


(2.73 ±0.21) x 10~ 2 


11.7 


12.1 


2.70 x 10~ 2 


240 


270 


253.2 


139 


(1.49 ±0.13) x 10~ 2 


12.0 


12.4 


1.49 x 10~ 2 


270 


300 


283.4 


964 


(7.87 ±0.25) x 10~ 3 


11.9 


12.7 


7.60 x 10~ 3 


300 


320 


309.3 


371 


(4.55 ±0.24) x 10~ 3 


12.1 


12.3 


4.46 x 10~ 3 


320 


350 


333.6 


292 


(2.40 ±0.14) x 10~ 3 


12.5 


13.0 


2.75 x 10~ 3 


350 


390 


367.6 


2682 


(1.41 ±0.03) x 10~ 3 


11.6 


12.4 


1.44 x 10~ 3 


390 


430 


407.8 


1445 


(7.62 ±0.20) x 10~ 4 


12.0 


13.0 


7.16 x 10" 4 


430 


470 


447.9 


689 


(3.64 ±0.14) x 10" 4 


12.6 


13.6 


3.70 x 10" 4 


470 


510 


488.0 


408 


(2.16 ±0.11) x 10~ 4 


13.1 


14.1 


1.97 x 10~ 4 


510 


550 


528.1 


219 


(1.16 ±0.08) x 10~ 4 


13.4 


15.0 


1.07 x 10~ 4 


550 


600 


572.2 


244 


(6.11 ±0.39) x 10~ 5 


13.6 


16.0 


5.59 x 10~ 5 


600 


700 


639.4 


192 


(2.40 ±0.17) x 10~ 5 


14.7 


17.5 


2.05 x 10~ 5 


700 


800 


739.8 


49 


(6.10 ±0.87) x 10~ 6 


17.1 


20.4 


5.19 x 10~ 6 


800 


1400 


878.8 


20 


(4.06 ±0.91) x 10~ 7 


22.4 


28.0 


2.92 x 10~ 7 
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The uncertainty in the theoretical prediction of this ratio due to the choice of PDF is 
less than 3%, and 6% from the choice of renormalization and factorization scale (excluding 
fi = 0.25.E™ ax ). The luminosity matching error only contributes to those bins where the data 
from triggers Jet_30 and Jet_50 overlaps with the data from triggers Jet_85 and Jet_115 (i.e. 
for masses between 300 and 350 GeV/c 2 ). The errors from the vertex selection cancel when 
the data in a bin come from the same trigger for each of the cross sections. The errors due 
to the unsmearing and the (partially) correlated part of the energy scale are assumed to be 
correlated for the two cross sections and mostly cancel out leaving small errors (<^C 1%). In 
addition the uncertainty due to the hot cell restoration is assumed to be correlated between 
the two cross-sections. All other errors are assumed to be uncorrelated between the two 
measurements. For the purposes of calculating a covariance matrix, the correlated energy 
scale and unsmearing errors are assumed to be fully correlated as a function of mass. 

The resulting cross section ratios are given in Table pCLVIIj and plotted in Fig. |93|. Taking 
the ratio of the cross sections reduces the systematic uncertainties to less than 10%. The 



correlations of the systematic uncertainties are given in Table XLIV [55 



TABLE XLVII. The ratio K(\rf et \ < 0.5)/«(0.5 < |r? jct | < 1.0). The systematic uncertainties 
are the sum in quadrature of the uncertainties from the ±1<7 variations in the energy calibration, 
the unsmearing, the vertex corrections, luminosity matching, jet selection, and the uncertainty 
in the luminosity. Also shown is the jetrad prediction with [i = 0.5-Ej? ax , lZ S ep = 1-3, and the 
CTEQ3M PDF. 

Mass Bin Ratio of Mass Spectra Theoretical 



(GeV/c 2 ) k(\ rf ct | < 0.5)//s(0.5 <\rf ct \ < 1.0) Prediction 



Min. 


Max. 


(± stat. error ± syst. errror) 




200 


220 


0.613 ± 0.039± 0.037 


0.616 


220 


240 


0.614 ± 0.050± 0.030 


0.621 


240 


270 


0.570 ± 0.051± 0.029 


0.627 


270 


300 


0.568 ± 0.030± 0.027 


0.635 


300 


320 


0.610 ± 0.034± 0.050 


0.642 


320 


350 


0.705 ± 0.044± 0.058 


0.648 


350 


390 


0.672 ± 0.020± 0.032 


0.657 


390 


430 


0.593 ± 0.022± 0.030 


0.667 


430 


470 


0.708 ± 0.036± 0.037 


0.676 


470 


510 


0.690 ± 0.046± 0.036 


0.685 


510 


550 


0.620 ± 0.058± 0.033 


0.693 


550 


600 


0.634 ± 0.065± 0.033 


0.701 


600 


700 


0.647 ± 0.074± 0.034 


0.710 


700 


800 


0.608 ± 0.141± 0.035 


0.718 


800 


1400 


0.705 ± 0.246± 0.046 


0.711 
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FIG. 93. The ratio of cross sections for | rf ct | < 0.5 and 0.5 < | rf ct | < 1.0 for data (solid circles) 
and theory (various lines). The error bars show the statistical uncertainties. The shaded region 
represents the ±1<t systematic uncertainties about the prediction. The effects on the prediction of 
changing the renormalization scale are also shown. 

H. Comparison of Data with Theory 

Figure |54] shows the ratio (Data — Theory)/Theory for \r]^ Qt \ < 1.0 and the jetrad 
prediction using CTEQ3M with \i = 0.5-E™ ax . The effect of varying the renormalization 
scale in the prediction is also shown; all are in good agreement except for fi = 0.25-E™ 



max 



which lies approximately 30% below the data. Figure 95 shows (Data — Theory)/Theory for 
jetrad predictions with different choices of PDFs. Given the experimental and theoretical 
uncertainties, the predictions can be regarded as being in good agreement with the data. 



Figure [SB] shows that the data and jetrad predictions are in agreement for | rf^ | < 0.5 and 



0.5 < \rf ct \ < 1.0. The data are also in agreement, within the uncertainties, with the cross 
section measured by CDF [jLs|| . 

A x 2 can be calculated for each of the comparisons between the data (cross sections and 
ratio of cross sections) and the predictions. The x 2 is given by: 

x 2 -J2 s iViJ% ( 14 - 16 ) 

where Si is the difference between the data and theory for mass bin i, and Vij is element i,j 
of the covariance matrix: 

Vij = p ir Aai- Aaj, (14.17) 

where Act is the sum of the systematic error and the statistical error added in quadrature 
if i = j and the systematic error if i ^ j, and p,j is the correlation between the systematic 
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FIG. 94. The difference between the data and the prediction (jetrad) divided by the prediction 
for |r/ jet | < 1.0. The solid circles represent the comparison to the calculation using CTEQ3M with 
fi = 0.5.Ep ax . The shaded region represents the ±lu systematic uncertainties about the prediction. 
The effects of changing the renormalization scale are also shown (each curve shows the difference 
between the alternative prediction and the prediction using CTEQ3M with /i = 0.5i?™ ax ). 



uncertainties of mass bins i and j as given in Table [XL1VI The systematic uncertainty is 
given by the fractional uncertainty times the theoretical prediction. The resulting y 2 values 
are given in Table |XLV11J| for all of the theoretical choices described above. The choice of 
PDF and renormalization scale is varied; all choices are in good agreement with the data, 
except for fi = 0.25i?™ ax which is excluded by the data. 



I. Compositeness Limits 

The ratio of the mass spectra can be used to place limits on quark compositeness (Sec- 
tion [IV B| ). Currently there are no NLO compositeness calculations available; therefore a 
LO event generator (pythia) is used to simulate the effect of compositeness. The ratio of 
these LO predictions with compositeness, to the LO with no compositeness, is used to scale 



the jetrad NLO prediction, shown in Fig. Tl . 



The data show no evidence of compositeness and are used to set 95% confidence level 
limits on A^ L . This was done using the same method that was used to extract compositeness 
limits from the dijet angular distribution (Section [XIII E|) . Figure |)8] shows the probability 



distribution for a theoretical prediction obtained using jetrad with the CTEQ3M PDF 
and a renormalization scale of \i = i?™ ax . The 95% CL. limit on the compositeness scale is 
At L > 2.7 TeV. 

Limits were also set for several different theoretical choices of PDF and renormalization 
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FIG. 95. The difference between the data and the prediction (jetrad) divided by the prediction 
for |r/ Jct | < 1.0. The solid circles represent the comparison to the calculation using [i = 0.5E'™ ax 
and the PDFs CTEQ4M, CTEQ4HJ, MRS(A'), and MRST. The shaded region represents the 
zbla systematic uncertainties about the prediction. 
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FIG. 96. The difference between the data and the prediction (jetrad) divided by the prediction 
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The shaded region represents the ±la systematic uncertainties 
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TABLE XLVIII. The calculated x 2 for K(\r? et \ < 1.0) (15 degrees of freedom), 
n(\7]i ct \ < 0.5) (18 d.o.f.), and k(0.5 < \rf et \ < 1.0) (15 d.o.f.) and for the ratio 
K,(\if et | < 0.5)/k(0.5 < |r/j ct | < 1.0). The probability of obtaining a larger x 2 is also given. 



PDF 




K(\rf ct \ 


< 1.0) 


K(\rf ct \ 


< 0.5) 


«(0.5 < \ rf 


et | < 1.0) 


Ratio 






x 2 


Prob. 


x 2 


Prob. 


x 2 


Prob. 


x 2 


Prob. 


CTEQ3M 


0.25^ ax 


24.7 


0.05 


26.4 


0.09 


38.3 


0.001 


29.1 


0.02 


CTEQ3M 


0.50£f ax 


5.7 


0.98 


11.2 


0.89 


8.9 


0.88 


14.1 


0.52 


CTEQ3M 


0.75£f ax 


6.1 


0.98 


11.2 


0.89 


9.1 


0.87 


13.6 


0.56 


CTEQ3M 


l.OOEf 3 * 


6.3 


0.97 


12.1 


0.84 


9.2 


0.87 


13.3 


0.58 


CTEQ3M 


2ME!p ax 


6.0 


0.98 


12.5 


0.82 


11.5 


0.71 


13.0 


0.60 


CTEQ3M 


0.25^/xiX2S 


12.7 


0.63 


28.7 


0.05 


10.2 


0.81 


14.9 


0.46 


CTEQ3M 


0.50-^X1X25 


6.1 


0.98 


14.5 


0.70 


8.8 


0.89 


13.8 


0.54 


CTEQ3M 


1.00^/XlX2S 


7.7 


0.93 


13.4 


0.77 


13.3 


0.58 


14.3 


0.51 


CTEQ4M 
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8.3 


0.91 
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0.51 
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13.9 


0.54 
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0.19 
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0.19 
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FIG. 97. The ratio of cross sections for | rf ct \ < 0.5 and 0.5 < | rf \ < 1.0 for data (solid circles) 
and theoretical predictions for compositeness models with various values of A^ L (various lines, see 
Section IV B for model details). The error bars show the statistical uncertainties. The shaded 
region represents the ±1<t systematic uncertainties about the jetrad prediction. 
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FIG. 98. The probability distribution (solid curve) P(a\t;')P(£')/Q(oo) for the theoretical pre- 
diction jetrad with [i = £vp ax . The dashed curve shows the integral of the probability distribution 
and the dotted line shows the 95% C.L. limit on the compositeness scale, 2.73 TeV. The most prob- 
able value for the compositeness scale is A~^ L = oo. 
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scales for both the NLO jetrad and LO compositeness predictions. The limits on the 



compositeness scale are summarized in Table |XLLX| . The dijet mass spectrum rules out 
quark compositeness models at the 95% confidence level where Aj L is below 2.7 TeV and 
tC LL is below 2.4 TeV. 

TABLE XLIX. The 95% confidence level limits in TeV for the left-handed contact composite- 
ness scales for different models. 
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Limits on models with color-singlet (octet) vector or axial contact interactions were also 
set using an analytic LO calculation instead of the pythia event generator. The result- 
ing limits are given in Table |]. The limits on the scale of Ay g can be converted into limits 
on a flavor-universal coloron [|5{|, resulting in a 95% CL. limit of M c / cot# > 837 GeV/c 2 
(see Section |1V C| for a description of the theory). 

The robustness of the confidence limits are tested in two ways. The first assumes that 
the systematic uncertainties are completely uncorrelated as a function of mass, which results 
in a degradation of the limit by 10 GeV (negligible compared to the scale of the limit). The 
second doubles the size of the systematic uncertainty, which results in a degradation of the 
limit by 20 GeV. 



J. Conclusions 



We have measured the inclusive dijet mass spectrum for a pseudorapidity range of | r]i ct \ < 
1.0 and 200 < M < 1400 GeV at ^/s = 1.8 TeV to an accuracy of 10% to 30% as a function of 
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TABLE L. 95% confidence level limits in TeV for different contact compositeness scale for 
different models calculated using an analytic LO prediction |38| (see Section IV B for a description 
of the models). 
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mass. QCD NLO predictions, using several PDFs, show good agreement with the observed 
inclusive dijet mass spectrum. 

The ratio of the inclusive dijet mass cross sections for \if et \ < 0.5 and 0.5 < \rf et \ < 
1.0 has also been measured with a systematic uncertainty that is less than 10%. The 
data distributions are in good agreement with NLO QCD predictions. Models of quark 
compositeness with a contact interaction scale of less than 2.2 TeV are excluded at the 95% 
confidence level. 
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XV. CONCLUSIONS 

We have presented a series of measurements of high energy jets at the Fermilab Tevatron 
which are sensitive to the various components of QCD predictions: the parton distributions, 
the matrix elements, and the scales. Measurements of the cross section as a function of 
jet E^, and dijet invariant mass have been presented. By taking the ratio of the inclusive 
cross sections at two energies, both the experimental errors and the sensitivity to the parton 
distributions were reduced, providing a stringent test of the E T dependence of the QCD 
matrix element at next-to-leading order. By looking at both the dijet angular distribution 
at fixed mass and the ratio of dijet invariant mass distributions in two different rapidity 
ranges, we have again minimized the experimental uncertainties and tested the angular 
dependence of the matrix element calculation. 

We have made the most precise measurement to date of the inclusive jet cross section 
for E? > 60 GeV at \fs = 1800 GeV. No excess production of high-E T jets is observed. 
QCD predictions are in good agreement with the observed cross section for standard parton 
distribution functions and different renormalization scales (/i = 0.25-2.00-Er where = 
E™ ax and E 3 ^ ). We have also made the most precise measurement to date of the ratio of 
the inclusive jet cross sections at \/s = 630 and 1800 GeV. The NLO QCD predictions 
yield satisfactory agreement with the observed data for standard choices of renormalization 
scale or PDF. In terms of the normalization however, the absolute values of the standard 
predictions lie consistently and significantly higher than the data. 

We have measured the dijet angular distribution over a large angular range and the inclu- 
sive dijet mass spectrum for a pseudorapidity range of \i]i ot | < 1.0. QCD NLO predictions, 
using several PDF's, show good agreement with the observed inclusive dijet mass spectrum. 
Since we found good agreement, the data have permitted us to provide sensitive limits on 
the existence of possible non-standard model phenomena. 
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APPENDIX: x 2 STUDIES 

In this paper we have made quantitative x 2 comparisons between theoretical predictions 
and data to determine which predictions provide better agreement. The systematic uncer- 
tainties in the inclusive jet cross section (Sections [X| and |XII|) and the dijet mass spectrum 



(Section |XIV|) are highly correlated. An inappropriate definition of the uncertainties in 
X 2 analyses may result in theoretical predictions that have an average normalization be- 
low the data yielding a better fit (Peelle's Pertinent Puzzle ||60|| ). The first section of this 
Appendix describes alternative methods for calculating the x 2 an d our choice of an appro- 
priate method. The second section describes studies of the probability distributions for the 
analyses presented in this paper. 



1. Definition of x 2 

The x 2 is given by 

x 2 = E^7^ ( A1 ) 



'•j 



where 5i is the difference between the data and the expected cross section for bin i, and 
is element i, j of the covariance matrix, with each element given by: 

Va = Pij (Aaf at Aaf % + Aaf s Aaf s ) (A2) 

where 5ij is the Kronecker delta function, py = 1 for i = j, and is the correlation of the 
systematic uncertainties between cross section bins i and j. 

The analyses presented in this paper is based on using the fractional systematic uncer- 
tainties in each bin, but there are several ways of calculating the impact of the absolute 
systematic uncertainty on the x 2 values. We can use: 

1. Fractional uncertainty multiplied by the observed cross section. 



2. The fractional uncertainty multiplied by a smooth fit to the observed cross section |6T 
(which is normalized to the observed integrated cross section). 

3. The fractional uncertainty multiplied by a theoretical prediction. 

This Appendix discusses these choices. In previous publications of the inclusive jet cross 
section M and the dijet mass spectrum |25j the x 2 values were calculated using the first 



option. 

The choice of calculation for the absolute systematic uncertainty used in the x 2 is i n_ 
vestigated using the measurement of the dijet mass spectrum (Section |XIV|) . A theoretical 



prediction, called the Ansatz (A), based on a fit to the observed cross section (Fig. ^) is 
obtained by normalizing the fit to the observed integrated cross section (cf. Option 2). We 
also define a Floating Ansatz (FA) through a multiplicative factor X that is used to change 
the normalization of the Ansatz (FA = XA). A comparison between the Ansatz and the 
data is given in Fig. 
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FIG. 99. The difference between the dijet mass cross section for 1 771,2] < 1-0 and the Ansatz 
(see text). The dashed line shows the best fit obtained by using the standard x 2 an d absolute 
systematic uncertainties obtained using the product of the fractional systematic uncertainties and 
the measured cross section in each bin (Option 1). 
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If the systematic uncertainty is given by the product of the fractional uncertainty and the 
observed cross section in each bin (Option 1), the minimum value of the x 2 °f the Floating 
Ansatz is obtained for a normalization of X = 0.965 (the dashed line in Fig. |99|). This 
is clearly not the best visual fit to the observed cross section. When this test is repeated 
using Option 2, the preferred normalization is X = 1.0 (Fig. |100| ). Using several different 
predictions from jetrad (Option 3) also yields X = 1.0 as a best fit (not shown). 



N * 14 




0.95 1 1.05 

X (Normalization) 

FIG. 100. The x 2 for the Ansatz as a function of the floating normalization X for Option 2 
(see text). The minimum x 2 is obtained for a normalization of X = 1.0. The short vertical line 
indicates a normalization of X = 0.965, illustrating the bias of Option 1. 



Calculating systematic uncertainties using the observed cross section per bin introduces 
a statistical component to the systematic uncertainty; i.e., when the cross section fluctuates 
to a small value in a given bin the absolute systematic uncertainty also fluctuates to a smaller 
value. The smaller values of cross section therefore appear to be more precise relative to 
any given theory. This bias has been called Peelle's Pertinent Puzzle pU|. 

We choose to rely on Options 2 and 3 for determining systematic uncertainties. The 
choice depends on the question that is posed. In our work we wish to 

"Determine the probability that the theoretical prediction could have produced 
the observed number of events." 

This requires that we determine the systematic uncertainties using the theoretical predictions 
(Option 3). For example, if we underestimated the integrated luminosity the number of 
predicted events would also be underestimated. 
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This choice of x 2 definition means that the current results differ from those published 
previously for the inclusive jet cross section || and the dijet mass spectrum |25[. Table [C| 
summarizes the differences in y 2 values for the dijet mass analysis. The x 2 values in Table [LI 
are calculated using the same luminosity definition as given in Ref. [p5fl , and differ from 
those given in Table [XLV111| . The \ 2 changes most for theoretical predictions with the 
largest normalization differences with the data. If the theoretical prediction has a smaller 
normalization than the data then the size of the systematic uncertainties are reduced, hence 
increasing the value of the x 2 ■ 



TABLE LI. The \ 2 f° r the cross section in dijet mass for \ rf et \ < 1.0 (15 degrees of freedom). 
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2. Probabilities 



The probability that a given theoretical prediction agrees with the data for a given x is 
calculated assuming that the x 2 is given by the standard distribution fl62j: 



/ (z; n) 



x W 2 -i) eX p (-x/2) 
2™/ 2 r (ra/2) ' 



(A3) 



where n is the number of degrees of freedom (d.o.f.) in the data. The probability of getting 
a value of x 2 larger than the one obtained is then given by 
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P(x 2 ;n) 



X 



f [x] n) dx. 



(A4) 



Hence, for the probabilities quoted in Sections [Xl] and |XIV1 to be reliable, the x 2 distribution 
for comparisons between theoretical predictions and the data must follow Eq. [A3| . 

The distribution of x 2 f° r comparisons with the dijet mass spectrum was tested by 
developing a Monte Carlo program that generates many trial predictions based on the Ansatz 
(with a total of 15 bins, or 15 degrees of freedom). The first step is to generate trials based 
on statistical fluctuations. The true number of events per bin is given by the Ansatz. The 
trial spectra are then generated for each bin by sampling a Poisson distribution with a mean 
defined by the true number of events. The x 2 f° r each of these trials is calculated using the 
difference between the true and the generated values. Figure |101| shows the x 2 distribution 



for all of the generated trials. The distribution is fitted to Eq. |A3|, with the best fit obtained 
for n = 15.08 ± 0.20, which is consistent with the expected value of n = 15 for a normalized 
distribution of bins. 
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FIG. 101. The x 2 distribution generated by sampling the Ansatz cross section using only sta- 
tistical fluctuations. The histogram shows the generated distribution, and the curve is a fit to the 



histogram using Eq. A3 . The fitted number of degrees of freedom is 15.08. 



The final step is to assume that the uncertainties are the same as the uncertainties in the 
measurement of the dijet cross section. Trial spectra are generated using these uncertainties 
in order to obtain a x 2 distribution (see the dotted curve in Fig. |102|) . It is clear that the 
X 2 distribution is very similar to that predicted by Eq. |A3|; hence any probability generated 
using Eq. [A4] should be approximately right. The resulting x 2 distribution was fitted using 
Eq. [A3] (Fig. [103]) and yielded n = 14.6 ± 0.2. 
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FIG. 102. The x 2 distribution generated by sampling the Ansatz cross section using only sta- 
tistical fluctuations (solid curve), and fluctuations based on the uncertainties in the dijet cross 
section as a function of dijet mass (dotted curve). 
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FIG. 103. The x 2 distribution generated by sampling the Ansatz cross section using all of the 
systematic uncertainties of the dijet cross section. The histogram shows the generated distribution 



and the curve is a fit to the histogram using Eq. A3. The fitted number of degrees of freedom is 
14.6 ±0.2. 
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The study of the x 2 distribution was repeated for the measurement of the inclusive jet 
cross section, which has 24 bins (d.o.f.). Figure 104 shows the resulting distributions for 
statistical fluctuations (solid curve) and the systematic uncertainties in the inclusive jet cross 
section (dotted curve). The two distributions agree for \ 2 values below approximately 15, 
and then begin to diverge. The distribution based on the cross section uncertainties has a 
longer tail than the statistical x 2 distribution. This implies that all the probabilities quoted 
in Section [XI] are slightly underestimated. 
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FIG. 104. The x 2 distribution generated by sampling the inclusive jet cross section Ansatz 
using only statistical fluctuations (solid curve), and fluctuations based on the uncertainties in the 
inclusive jet cross section (dotted curve). 



Finally, the ratio of inclusive jet cross sections (Section [XI I|) , was also examined with 

The resulting distribution is similar to the one 



the results of the study given in Fig. |105 



obtained for the inclusive jet cross section, with the distribution based on the uncertainties 
having a larger tail than the standard x 2 distribution. The maximum deviation between the 
probability obtained assuming the standard distribution and the measured distribution is 
2.9%, and probabilities quoted in Section [XI 1| will therefore be slight underestimates of the 
correct probabilities. 

The studies presented describe the x 2 comparisons made between the observed data and 
the theoretical predictions. We have demonstrated that they give an accurate representation 
of the probability of agreement between a given theoretical prediction and the data. 
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FIG. 105. The x 2 distribution generated by sampling the ratio of inclusive jet cross sections 
Ansatz using only statistical fluctuations (solid curve) , and fluctuations based on the uncertainties 
in the inclusive jet cross section (dotted curve). 
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